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1.  INTRODUCTION 


Materials  that  conduct  electrical  currents  very  poorly  even  in  the 
presence  of  strong  electric  fields  are  known  as  insulators  or  dielectrics. 
Such  materials  find  wide  use  in  cable  insulation,  high  voltage  standoffs, 
capacitors,  and  other  components.  When  a dielectric  material  is  sub- 
jected to  sufficiently  strong  electric  fields,  it  undergoes  rapid  struc- 
tural changes  and  its  conductance  changes  markedly.  The  phenomena 
associated  with  these  changes  are  referred  to  as  dielectric  breakdown. 
Breakdown  properties  are  frequently  studied  by  constructing  capacitors 
with  two  parallel  planar  electrodes  separated  by  a slab  of  dielectric 
material  or  with  one  planar  electrode  and  a second  nonplanar  electrode 
separated  by  a dielectric  slab. 

The  mechanisms  responsible  for  dielectric  breakdown  have  received 
attention  for  over  100  years  (1,2).  It  is  useful  to  distinguish  two 
principal  mechanisms,  thermal  and  electrical.  Thermal  breakdown  gener- 
ally occurs  in  those  situations  where  leakage  currents  are  large  and  the 
field  is  less  than  10^  V/cra.  In  thermal  breakdown,  leakage  current 
causes  a temperature  rise  which  enhances  the  electrical  conductivity, 
causing  still  further  temperature  rise,  the  process  continuing  until  the 
dielectric  bums,  fractures,  blisters,  or  otherwise  loses  its  integrity. 
Thermal  breakdown  is  generally  a relatively  slow  process,  with  a time 
scale  of  many  seconds.  Destruction  to  the  dielectric  involves  a size- 
able portion  of  the  insulator,  and  the  associated  light  emission  is  ther- 
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mal  in  character.  No  new  conduction  mechansim  has  to  be  hypothesized  in 
thermal  breakdown;  high  conductance  follows  from  the  positive  temperature 
dependence  of  the  conductivity  with  temperature  and  damage  occurs  through 
Joule  heating,  by  oxidation,  chemical  dissociation,  thermal  stress,  and 
mechanical  stress. 

Elecrrical  breakdo™  occurs  at  field  strengths  greater  than 
5 X 10^  V/cm.  It  is  accompanied  by  formation  of  a localized  gaseous 
channel  that  grows  within  a few  tens  of  nanoseconds  from  inception  and 
light  is  emitted  characteristic  of  the  gas  in  the  channel.  The  two 
mechanisms,  thermal  and  electric,  compete  with  each  other  in  a given 
material  under  given  environmental  conditions  (2) . 

Historically,  another  term  has  also  been  employed  in  describing 
breakdown  mechanisms,  "intrinsic  breakdown"  (2).  Intrinsic  breakdown 

refers  to  the  (unique)  ultimate  breakdown  strength  of  a material 

implying  that  each  material  has  a unique  ultimate  strength.  However, 
attempts  to  measure  a unique  strength  have  not  been  successful  (2). 

Thus,  the  most  carefully  annealed  single  crystals  do  not  give  higher 
breaKdown  strengths  than  crystals  with  strains.  Deliberately  introduced 
impurities  do  not  necessarally  reduce  breakdown  strength. 

The  present  study  is  concerned  with  the  theory  of  breakdown  in 
solid  insulators.  In  Chapter  II,  a short  review  of  earlier  experimental 
observations  and  theoretical  work  wi31  be  given  as  background  information. 
This  review  will  include  a discussion  of  a phenomenological  model  of 
Budenstein  (4)  which  will  provide  the  basis  for  the  present  study. 

Chapter  III  considers  one  aspect  of  this  model  as  it  might  pertain  to 
alkali  halides  and  similar  materials.  Chapter  IV  contains  the  major 
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effort  of  this  thesis,  a cal«,vilation  of  the  ground  state  energy  of  the 
doubly  negative  hydrogen  molecular  ion.  Chapter  V considers  the  impli- 
cations of  the  results  of  Chapter  IV  to  the  general  breakdown  model 
proposed  by  Budenstein,  and  recommendations  are  given  for  further  theo- 
retical work  which  needs  to  be  done  in  evaluating  this  model. 


II.  BACKGROUND  ON  DIELECTRIC  BREAKDOWN 


In  order  to  put  the  present  study  into  proper  perspective,  it  is 
necessary  to  Indicate  the  results  of  previous  work.  There  exists  a 
copious  literature  on  breakdown  in  solids,  including  several  cogent 
review  articles  and  monographs  (1,7, 8, 9).  Our  review  here  will,  of 
necessity,  emphasize  those  points  which  relate  most  directly  to  the  model 
we  wish  to  investigate,  and  without  which  a proper  understanding  of  the 
differences  between  this  model  and  others  proposed  previously  cannot 
he  achieved.  This  review  will  first  describe  a broad  range  of  breakdown 
phenomena  obtained  from  studies  of  thin  film  and  bulk  samples,  and  will 
then  examine  the  existing  theories  of  electric  breakdown. 

A.  Experimental  Background 

fl 

Dielectric  breakdown  occurs  at  field  strengths  of  about  10°  V/cm 
in  homogeneous  dielectrics.  A statistical  nature  rf  breakdown  is 
indicated  by  the  wide  variation  in  the  breakdown  field  strength  found 
in  a given  substance.  Indeed,  in  a thin  film  capacitor  that  breaks 
down  with  isolated  single  hole  channels,  the  breakdown  voltage  may 
easily  vary  by  100%  (10).  Similar  variations  have  been  observed  in 
carefully  grown,  annealed  single  crystals  of  various  alkali  halides(l]). 
Overlapping  of  the  distributions  of  the  breakdown  field  strengths  can 
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be  expected  for  most  homogeneous  insulators,  even  though  the  materials 
may  be  structurally  dissimilar. 

^•Jhen  dielectric  breakdown  occurs,  the  dielectric  changes  abruptly 
from  a nonconducting  state  to  one  in  x^hlch  a localized  current  passes 
through  the  material,  the  local  resistivity  decreasing  by  about  ten 
orders  of  magnitude  (12).  The  conduction  process  is  always  accompanied 
by  local  destruction  of  the  dielectric.  If  the  device  is  a thin  film 
capacitor,  the  electrode  material  may  be  evaporated  or  blown  back  from 
the  breakdown  channel,  and  the  capacitor  left  non-shorted.  Thus 
multiple  breakdown  studies  may  be  performed  on  a single  sample.  In  a 
macroscopic  insulator,  the  localized  region  of  destruction  of  a breakdown 
becomes  a high  conductance  gaseous  path  through  the  sample.  Local 
removal  of  the  electrode  material  does  not  isolate  this  path  because 
surface  flashover  from  one  electrode  to  the  gaseous  channel  occurs  upon 
re-application  of  the  field  before  the  solid  dielectric  can  break  down 
in  another  place. 

The  destruction  within  the  dielectric  causes  the  formation  of  a 
channel  through  the  material,  extending  completely  through  the  dielec- 
tric. The  channel  may  contain  frequent  branching,  giving  a "tree"  or  a 
"star"  pattern.  This  pattern  is  especially  prominent  for  nonuniform 
(point-plane)  applied  field  geometries.  In  other  cases,  such  as 
parallel  plate  geometry,  a straight  channel  with  essentially  no 
deviations  may  occur,  or  tree  channels  formed,  normally  from  the  cathode 
side.  The  formation  of  one  or  the  other  type  of  channel  may  be  polarity 
dependent . 


The  following  illustrations  are  the  most  graphic  manner  in  which 
the  evidence  of  the  breakdown  channel  may  be  presented.  The  orientation 
of  the  view  is  along  the  direction  of  the  electric  field  lines  in  each 
case  for  Figs.  1-4.  Figures  1,  3,  and  4 are  from  thin  film  capacitor 
studies,  while  Fig.  2 is  in  a bulk  sample. 

In  Fig.  1,  the  typical  single-channel  unbranched  breakdown  is 
illustrated.  This  type  of  breakdown  is  normal  in  thin  film  capacitors, 
where  the  voltage  is  applied  relatively  slowly  to  the  capacitor.  The 
central  region  is  the  main  channel,  the  next  region  the  dielectric  where 
the  upper  electrode  has  been  removed  by  the  breakdown,  and  the  outermost 
region  the  upper  electrode.  Two  points  of  particular  note  appear  in  this 
photograph:  Firstly,  the  appearance  of  the  electrode  is  similar  to  that 
of  an  explosion  cavity,  with  much  peeling  back  of  the  electrode  from  the 
channel,  and  lines  of  metal  fragments  extending  quite  some  distance  from 
the  center  of  the  channel.  Secondly,  the  dielectric  appears  to  have 
suffered  structural  damage  in  the  regions  visible  in  the  photograph. 

The  "star"  pattern  of  breakdown  channel  occurs  frequently  in 
point-plane  geometries,  but  can  occur  also  in  other  electrode  configur- 
ations. The  breakdown  illustrated  in  Fig.  2 makes  the  origin  of  the 
name  self-evident.  Although  the  scales  of  the  two  illustrations  are 
different,  it  will  be  noted  that  the  breakdown  event  in  Fig.  1 is 
significantly  different  from  that  in  Fig.  2.  In  the  former,  the  main 
channel  was  the  dominant  feature,  whereas  the  side  channels  are  much 
more  prominent  in  the  latter  case.  The  "tree"  pattern  of  breakdown  is 
similar  to  the  "star",  but  shows  more  extensive  side-branching  and  a 


Fig.  1.  Normal  single  channel  breakdown  in  a thin  film 
showing  damage  to  dielectric  and  electrodes. 
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much  lower  degree  of  overall  symmetry  in  the  positioning  of  the  side 
channels.  The  "tree"  pattern  is  not  illustrated. 

The  next  two  illustrations  show  explicitly  that  breakdown  channels 
may  serve  as  nucleation  sites  for  additional  breakdowns  in  thin  film 
capacitors.  Both  of  these  figures  are  somewhat  atypical;  the  sort  of 
clustering  shown  in  Fig.  3 apparently  occurs  only  in  CaF2.  Clustering 
of  breakdown  channels  can  also  occur  in  other  materials;  that  of  Fig.  4 
is  in  SiO,  and  is  the  more  common  sort  of  clustering. 

The  process  of  dielectric  breakdown  in  liquids  seems  to  require 
production  of  a gaseous  channel  (13) . These  gaseous  channels  frequently 
exhibit  a structure  similar  to  the  beading  which  may  occur  in  channels 
in  solids.  A gaseous  channel  appears  to  be  an  essential  part  of  the 
breakdown  in  solids.  This  observation  is  certainly  compatible  with  the 
experimental  observation  that  there  exists,  in  thin  film  studies,  some 
minimum  voltage  at  which  the  breakdown  event  terminates.  For  a given 
capacitor,  this  voltage  minimum  is  quite  sharp,  and  seems  to  be  the  same 
for  similar  capacitors,  regardless  of  thickness.  This  is  explicable  in 
terms  of  a gaseous  discharge,  which  exhibits  the  same  effect. 

Examination  of  the  breakdown  region  in  both  thin  film  and  bulk 
systems  provides  evidence  for  the  existence  of  a high  pressure  gas 
internally  during  the  breakdown  process  (14).  Material  is  thrown  away 
from  the  breakdown  site  (as  in  Fig.  1),  fracture  lines  emanate  from  the 
breakdown  region,  and  electrodes  are  deformed.  Inspection  of  the  region 
of  destruction  in  thin-filra  capacitors  for  a given  breakdown  event 
reveals  a central  region  of  complete  destruction  of  dielectric  and  both 
electrodes.  The  outer  (evaporated)  electrode  is  further  eroded  or  blown 


Fig.  3.  Multiple  breakdown  pattern  observed  in  a thin-film 
sample  of  CaF2.  This  type  of  multiple  breakdown  seems  to  occur 
only  in  this  material. 
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back  from  this  central  region,  giving  a cratered  appearance,  The  entire 
configuration  suggests  that  a high  pressure,  high  temperature  gas 
existed  in  the  central  region  during  breakdown  (10) . 

The  thickness  dependence  of  the  breakdown  strength  has  been  studied 
for  both  thin  films  and  for  thicker  specimens.  For  thin  films  the 
field  strength  increases  as  sample  thickness  decreases  (1).  A similar 
but  less  pronounced  variation  was  observed  in  thicker  samples  (15) . 

Another  parameter  which  has  been  studied  is  the  temperature  depen- 
dence of  breakdown.  It  is  generally  found  in  the  case  of  electric 
breakdowns  that  the  small  variation  in  breakdown  field  strength  with 
temperature  acts  to  lower  the  field  strength  with  increasing  tempera- 
ture (11) . 

Breakdown  strength  increases  slightly  with  the  rate  of  application 
of  voltage  to  the  specimen  (15,16).  This  effect  seems  to  be  found  in 
all  materials,  and  persists  in  both  thin  film  and  bulk  samples. 

In  a study  by  Budenstein  and  Hayes  (3,  17),  several  observations 
were  made  which  relate  to  the  most  likely  sites  for  breakdown.  One 
observation  was  that,  in  Al-SiO-Al  capacitors,  imperfections  with  dia- 
meters of  ''5-1  pm  were  frequently  the  sites  where  individual  breakdowns 
would  occur.  The  identity  of  these  defects  was  not  discernable, 
however.  If  the  beam  of  an  electron  microscope  were  focussed  on  a defect 
and  the  voltage  across  the  capacitor  gradually  increased,  breakdown 
would  always  occur  at  the  irradiated  defect.  However,  if  none  of  the 
above  defects  were  within  the  illuminated  area,  breakdown  would  occur  at 
some  other  site,  out  of  the  field  of  illumination.  The  capacitors 

O O 

studied  varied  in  thickness  from  about  2000  A to  around  10000  A. 
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In  another  study  CIO) , it  was  found  that  the  material  of  the  elec- 
trode made  a significant  difference  in  breakdown  strength;  it  has  also 
been  noted  that  he  change  from  metallic  to  electrolytic  electrodes  makes 
a considerable  difference  in  the  pre-breakdown  behaviour  of  the  capaci- 
tor (18) . This  may  be  attributable  to  changes  in  the  work  function  at 
the  electrode-dielectric  Interface.  It  has  also  been  determined  that 
the  polarity  of  the  applied  voltage  influences  somewhat  the  character 
of  the  breakdown  (17).  The  effect  of  polarity  decreased  for  increasing 
dielectric  thickness. 

Studies  have  also  been  carried  out  which  analyze  the  light  emission 
produced  during  breakdown  (10).  In  the  case  of  thin  film  capacitors, 
Smith  found  that  the  spectral  character  of  the  light  emitted  was 
entirely  a gaseous  bright-line  spectra  of  the  electrode  and  dielectric 
materials.  This  study  was  carried  out  using  Al-SlO-Al  capacitors  pri- 
marally,  but  similar  results  were  obtained  using  Al-SiO-Cu  capacitors, 
and  capacitors  with  MgF2  as  the  dielectric  material.  Neutral  state  and 
second-ionization  state  radiation  was  observed  in  these  spectra,  but  the 
emission  spectra  of  the  singly  ionized  species  predominated.  The  light 
emission  spectra  was  observed  within  a time  of  about  10  nsec  from  the 
onset  of  breakdown  conduction  as  indicated  by  the  voltage  waveform.  A 
search  to  detect  any  thermal  component  in  the  radiation  indicated  that 
the  thermal  component  was  less  than  one  percent  of  the  total  emission, 
if  it  was  present  at  all.  Another  point  of  some  significance  is  that 
when  the  electrodes  were  of  different  materials,  the  spectra  of  both 
were  present  from  the  first,  in  addition  to  that  of  the  components  of 


the  dielectric. 
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Cooper  and  his  covorkera  at  Manchester  have  also  carried  out  a 

number  of  experiments  which  provide  information  on  the  breakdown 

transition  (19,20).  In  the  majority  of  their  studies,  the  specimens 

used  were  alkali  halide  single  crystals  about  0.5  mm  thick.  PrebreakJown 

—8 

currents  up  to  about  10  amperes  occurred  in  samples  for  field 
strengths  of  about  60  percent  of  the  breakdown  value.  The  transit  time 
for  electrons  in  the  samples  they  employed  were  determined  to  be  about 
5 nsec  (19).  In  these  experiments,  they  were  able  to  investigate  the 
prebreakdown  light  emission  to  within  about  2 nsec  of  voltage  breakdown. 
The  initial  Intensity  of  the  emitted  light  was  so  small  that  an  image 
intensifier  was  required  to  determine  the  light  output.  This  technique 
did  not,  unfortunately,  allow  them  to  determine  any  spectral  differen- 
tiation in  the  light.  Light  emission  was  first  detected  about  20  nsec 
before  the  voltage  drop  characteristic  of  breakdown.  At  19.5-20  nsec, 
a faint  glow  was  detected  near  the  cathode,  extending  for  some  distance 
along  the  cathode,  but  not  for  any  great  distance  into  the  dielectric. 
This  diffuse  light  gradually  extended  itself  across  the  dielectric, 
growing  somewhat  in  intensity.  This  core  completely  crossed  the 
dielectric  about  10  nsec  before  the  onset  of  breakdown.  The  light  emis- 
sion grew  rapidly  in  intensity  from  this  point  on.  It  was  also  obseirved 
that  the  light  from  the  core  region  showed  a headed  appearance,  typical 
of  many  breakdown  patterns.  This  core  region  did  indeed  become  the 
breakdown  channel. 

From  the  inability  of  doing  spectral  analysis  of  the  light,  Cooper 
was  not  able  to  say  what  physical  process  was  producing  the  light.  From 
investigation  of  temperature  variations  he  felt  that  it  was  not  arising 
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from  therraoluminescence  (22) , but  was  from  the  same  process  as  the 
breakdown,  or  from  a process  which  was  Intimately  connected  with  the 
formation  of  the  breakdown  channel. 

Cooper  ^ ^ (20)  also  studied  the  orientation  of  breakdown  channels 
in  single  crystal  samples  which  had  been  carefully  annealed.  The 
breakdown  channel  would,  for  a particular  alkali  halide,  normally  align 
in  a simple  crystallographic  direction  ( usually  [100]  but  [110]  was 
occasionally  observed) , even  when  the  applied  field  was  in  some  other 
direction.  A photograph  is  shown  of  a breakdown  channel  in  KCl  when  the 
applied  field  was  in  the  [110]  direction.  The  channel  started  out  in 
the  [100]  direction,  shifted  to  a [010]  direction,  and  then  changed  back 
to  its  original  direction.  Unannealed  and  polycrystalline  samples  did 
not  exhibit  this  behaviour,  but  would  break  down  in  the  direction  of 
the  applied  field. 

Grinberg  and  de  Grinberg  (21)  report  directional  effects  in 
breakdown.  In  this  case,  the  breakdown  channel  lies  predominantly  in  a 
[110]  direction,  with  side  branches  in  a [100]  direction.  No  treeing  or 
dendritic  growth  to  the  channel  was  noted;  each  side  branch  was  single. 
Their  crystals  were,  however,  prepared  by  cleaving  after  annealing,  and 
thus  are  not  comparable  with  those  of  Cooper's  work,  where  it  was  demon- 
strated that  the  presence  of  the  strains  caused  by  cleavage  could  alter 
the  breakdown  direction.  Grinberg  and  de  Grinberg  also  find  the  pres-r 
ence  of  free  atomic  alkali  metal  in  amorphous  form  about  the  breakdown 
channel . 

Another  point  determined  by  the  study  of  Al-SiO-Al  capacitors  by 
Budensteln  and  Hayes  (17)  is  of  some  interest  in  this  connection. 
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Selected  area  transmission  electron  diffraction  studies  of  breakdowns 
determined  that  "balls”  or  crystalline  silicon  of  about  1 pm  diameter 
were  always  found  In  the  breakdown  region.  X-ray  diffraction  studies  of 
these  capacitors  before  and  after  breakdown  showed  that  the  silicon  was 
a result  of  the  breakdown  process,  and  did  not  exist  within  the  specimen 
prior  to  breakdown. 

Similarly,  in  electron  beam  interaction  studies  of  Insulating  mater- 
ials, a chemical  reduction  of  cations  has  frequently  been  observed  at 
relatively  high  beam  currents  (23,24).  In  this  case,  the  cations  nor- 
mally form  into  amorphous  clusters  of  essentially  atomically  pure  mater- 
ial in  the  region  where  the  beam  current  is  relatively  large,  frequently 
globular  in  shape.  Also,  it  sometimes  occurs  that  breakdowns  occur  in 
these  samples;  these  breakdowns  have  channels  similar  to  those  of  the 
breakdowns  which  occur  in  capacitor  studies. 

In  another  study  of  thin  film  SiO  capacitors,  Tsuchida  and  Ueda 

report  the  "breakdown  strength"  (for  SiO  films  of  thickness  greater  than 
0 

about  200  A)  "increases  rapidly  with  decreasing  thickness  and  increased 
linearly  with  increasing  work  function  of  cathode  metal,  and  temperature 
had  a little  effect  on  the  breakdown  strength"  in  the  range  of  thickness 
investigated  (25).  For  extremely  thin  films,  they  found  the  breakdown 
strength  to  increase  somewhat  for  decreased  thickness.  For  such  very 
thin  films,  however,  the  current-voltage  characteristics  in  the 
prebreakdown  region  are  those  of  a Zener  diode,  and  the  conduction  was 
characteristic  of  Fowler-Nordheim  emission  at  the  cathode.  They  conclu- 
ded that  breakdown  in  such  very  thin  films  was  correctly  described  by 
Zener's  model. 
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B,  Theoretical  Background 

The  literature  dealing  vith  the  theoretical  aspects  of  dielectric 
breakdovm  is  extensive,  with  a number  of  review  articles  available  which 
cover  most  of  the  major  work  (1,11,26).  The  salient  points  of  most 
major  earlier  theoretical  work  will  be  Indicated  In  this  section, 
grouped  by  basic  mechanism. 

Avalanche  Theories 

Most  current  theories  of  breakdown  in  solid  dielectrics  begin  with 
the  assumption  of  electron  multiplication  by  an  avalanche  process.  The 
basic  argument  is  that  there  exist  conduction  electrons  wl.thln  the 
dielectric  capable  of  gaining  sufficient  energy  from  the  electric  field 
to  enable  them  to  ionize  valence  electrons  from  the  atoms  in  their  path 
by  a scattering  event.  Then  the  "new”  electrons,  as  well  as  the  incident 
electrons,  will  repeat  the  process,  and,  for  a sufficiently  strong  elec- 
tric field,  lead  to  a geometrical  progression  in  the  number  of  conduction 
electrons.  It  is  assumed  that  the  addition  of  electrons  to  the  conduc- 
tion band  is  the  direct  cause  of  the  high  currents  characteristic  of 
breakdown,  and  that  disruption  of  the  dielectric  follows  the  high 
current  by  Joule  heating  of  the  material.  We  shall  examine  briefly  two 
types  of  avalanche  theories,  one  in  which  the  applied  field  is  the  domi- 
nating factor,  the  other  in  which  space  charge  plays  an  important  role. 
Other  theoretical  treatments  will  also  be  considered  following  these. 

Probably  the  earliest  of  the  field-critical  theories  was  suggested 
by  von  Hlppel,  and  developed  analytically  by  Callen  (27,28).  For  a small 
applied  field,  conduction  electrons  do  not,  on  the  average,  gain  energy 
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from  the  field;  they  remain  near  the  bottom  of  the  conduction  band.  With 
increasing  field  strength,  the  average  energy  slowly  increases,  but  ener- 
gy losses  associated  with  scattering  events  still  dominate.  However,  at 
a sufficiently  high  field  strength,  the  energy  gain  between  collisions 
is  sufficient  to  lower  significantly  the  collision  probability.  Once 
beyond  this  "friction  barrier",  an  electron  gains  sufficient  energy 
to  cause  secondary  ionization  through  a collision  process  with  a bound 
electron  in  the  valence  band.  The  friction  barrier  is  in  the  range  of 

0.1-0. 2 eV  above  the  edge  of  the  conduction  band.  The  critical  field 
for  breakdown  is  taken  to  be  the  electric  field  at  which  the  average 
conduction  electron  energy  is  equal  to  the  friction  barrier  (1) . The 
principal  features  of  this  treatment  include: 

1.  The  field  strength  necessary  to  cause  breakdown  is  that  field 
capable  of  accelerating  the  average  conduction  electron  to 
ionization  energy,  regardless  of  its  Initial  energy. 

2.  The  breakdown  field  will  increase  with  an  increase  in  tempera- 
ture. This  is  a direct  result  of  the  enhanced  scattering  at  higher 
temperatures.  Smith  (15)  Indicates  an  increase  in  field  strength 
in  NaCl  of  from  about  1.1  MV/cm  at  0°K  to  1.6  MV/cm  at  500°K. 

3.  The  field  required  for  breakdown  in  independedt  of  thickness. 
Frohlich  suggested  that  it  is  not  necessary  for  the  average  energy 

of  the  conduction  electrons  to  exceed  the  friction  barrier,  but  that 
only  the  highest  energy  electrons  need  be  accelerated  to  the  lonlzatlcn 
energy  (1) . The  critical  field  for  breakdown  under  this  assumption  is 
taken  to  be  that  field  at  which  the  rate  of  energy  lost  to  the  lattice 
in  a collision  is  equal  to  the  rate  of  energy  gained  from  the  electric 
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field  for  an  electron  at  the  ionization  energy  Cl).  At  the  critical 
field  strength,  electrons  with  higher  initial  energies  will  attain 
energies  greater  than  the  ionization  energy.  Thus  the  high  energy  tall 
of  the  electron  energy  distribution  is  gaining  energy  from  the  field 
more  rapidly  than  energy  is  being  lost  in  collisions  with  the  lattice. 

We  thus  have  a supply  of  electrons  available  to  cause  the  avalanche 
multiplication  of  electrons.  Frohlich's  criterion  is  similar  to  that  of 
von  Hippel  and  Gallen  in  that  both  predict  a critical  field  strength 
that  depends  on  an  energy  balance.  The  difference  between  the  two  is 
the  emphasis  on  the  status  of  the  average  conduction  electron  by  von 
Hippel  and  Gallen  and  on  the  high  energy  tail  by  Frohlich.  As  the  field 
required  to  accelerate  the  average  electron  through  the  friction  barrier 
is  considerably  greater  than  that  required  to  accelerate  the  high-energy 
electrons  through  the  barrier  the  predicted  field  strength  of  Frohlich 
is  significantly  less  than  that  predicted  by  Gallen 's  criterion.  The 
other  basic  prediction  of  Frohlich’ s model  which  differs  from  that  of 
Gallen  is  a greater  temperature  dependence.  Frohlich  predicts  a 
breakdown  field  in  NaCl  of  0.4  MV/cm  at  O^K  and  of  0.8  MV/cm  at  bOO^K  (1) . 

The  theory  of  Forlani  and  Minna j a (1)  assumes  multiplication  of  the 
Fowler-Nordhelm  emission  current  from  che  cathode.  Through  a consid- 
eration of  current  continuity  obtained  by  equating  the  injection  and 
multiplication  currents,  they  obtain  a critical  field  for  breakdown 
which  varies  approximately  as  the  inverse  square  root  of  the  sample 
thickness.  Their  paper  does  not  treat  temperature  dependence. 

Zener  (29)  proposed  an  alternate  electron  multiplication  model 
which  was  further  developed  by  Franz  (30).  In  this  model,  explicit 
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avalanching  by  colliaions  does  not  play  any  significant  role,  but  the 
field  production  of  electrons  by  tunneling  is  taken  to  be  the  main 
source  of  conduction  electrons.  Once  an  electron  has  entered  the  conduc- 
tion band,  it  is  accelerated  by  the  field  in  the  direction  of  the  anode, 
Since  conduction  electrons  will  be  generated  throughout  the  material,  the 
greatest  current  density  will  occur  near  the  anode.  The  critical  field 
strength  increases  with  the  rate  of  application  of  the  field  and  with 
temperature.  The  mechanism  for  the  disruption  of  the  dielectric  is  again 
Joule  heating  by  the  current. 

In  each  of  the  preceding  theories,  the  basic  concept  was  that 
breakdown  was  an  Intrinsic  property  of  the  dielectric,  with  the  breakdown 
conduction  being  determined  largely  by  the  band  gap  of  the  material.  As 
experimental  evidence  relating  breakdown  in  thin  film  and  in  bulk  single 
crystal  systems  became  available,  these  ideas  became  more  difficult  to 
support.  In  the  past  decade,  an  increasing  emphasis  has  been  placed  on 
the  statistical  nature  of  breakdown  events,  reflecting  the  comparatively 
wide  range  of  breakdown  fields  observed  in  even  the  most  carefully 
prepared  specimens  (9) . The  theories  which  have  been  proposed  most  re- 
cently have  considered  the  role  of  space  charge  in  the  breakdown  process 
Since  space  charge  distributions  depend  strongly  on  the  specimen  geom- 
etry, its  past  electrical  history,  leakage  currents,  and  the  presence  of 
defects  and  mechanical  strains,  their  variability  can  well  account  for 
the  fluctuations  observed  in  breakdown  strengths,  The  following  theor- 
ies each  consider  the  effects  of  space  charge  on  the  process,  with  less 
emphasis  being  placed  on  the  occurrence  of  avalanching. 
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O'Dwyer  C31"33)  took  the  space  charge  into  account  by  requiring  the 
continuity  of  current,  assuming  that  the  current  is  carried  solely  by 
electrons.  His  work  also  included  physically  meaningful  boundary  con- 
ditions at  the  electrodes,  including  electrode  injection  of  carriers, 
and  modifications  to  the  field  due  to  the  space  charge.  This  model  still 
considers  that  colllslonal  ionization  plays  a significant  role.  The 
results  obtained  in  (31)  predict  a sharply  rising  field  at  the  cathode 
when  the  overall  field  exceeds  some  critical  value.  The  curve  obtained 
here  requires  ionization  rates  which  are  not  unreasonably  high;  the  crit- 
ical field  (expressed  in  dimensionless  form)  varies  approximately  lin- 
early with  the  dimensionless  inter-electrode  separation.  The  results 
of  (33)  are  much  more  Involved,  and  require  numerical  or  graphical  eval- 
uation. 

Watson  and  Heyes  (16)  have  considered  the  build-up  of  space  charge 
by  ionic  migration  to  the  cathode  to  predict  the  field  strength  for 
breakdown,  and  the  dependence  of  breakdown  field  on  the  rate  of  appli- 
cation of  the  field.  The  critical  field  yielded  by  this  treatment  seems 
to  exceed  those  normally  accepted  as  bulk  breakdown  field  strengths,  but 
the  critical  field  obtained  by  their  approach  is  the  field  at  the 
cathode.  This  may  correspond  to  a critical  field  for  the  build-up  of 
space  charge  near  the  cathode.  Good  agreement  between  theory  and  exper- 
iment is  claimed. 

DiStefano  and  Shatzkes  (6)  have  considered  an  impact-ionization 
model  in  which  the  effects  of  recombination  are  taken  into  account. 

Space  charge  is  created  which  enhances  the  field  in  the  cathode  region. 
Ultimately,  an  instability  develops  as  electrons  Injected  from  the 


cathode  gain  energy  from  the  field.  Impact  ionization  of  the  lattice  by 
the  electrons  enhances  the  positive  space  charge  near  the  cathode,  which 
increases  the  cathode  field,  which  enhances  the  electron  emission  at  the 
cathode,  etc.  In  the  case  of  Si02  their  model  predicts  a current  density 
for  the  instability  (breakdown)  which  increases  rapidly  for  films  with 

O 

thickness  below  200  A. 

The  impact  ionization  (avalanche)  models  seek  to  explain  the  obser- 
ved increase  in. current  accompanying  electric  breakdown.  Such  a current 
rise  is  taken  as  a necessary  condition  and  occurs  prior  to  any  damage. 
Damage  to  the  dielectric  is  assumed  to  follow  the  high  current  and  to  be 
thermal  in  nature.  The  direction  of  the  current  should  be  the  same  as 
the  field  direction.  Light  emission  should  be  either  a recombination 
radiation  in  the  solid  or  a thermal  emission  from  the  solid. 
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Critique  of  Avalanche  Theories 

Having  examined  both  the  experimental  evidence  and  the  various  av- 
alanche theories  of  electric  breakdown,  it  is  now  possible  to  summarize, 
The  avalanche  approach  has  several  virtues,  including: 

1.  Avalanche  theories  explain  the  large  increase  in  current  that 
is  characteristic  of  breakdown. 

2.  They  explain  localization  because  small  geometric  variations 
can  give  rise  to  high  local  fields  and  the  size  of  an  avalanche 
Increases  greatly  with  electric  field. 

3.  They  explain  the  insensitivity  of  breakdown  to  the  structure  of 
the  dielectric,  since  it  depends  only  on  the  gross  features  of  the 
electrical  band  structure. 

A.  They  provide  some  basis  for  understanding  the  thickness 
dependence  of  the  breakdown  field. 

Thus  there  are  a number  of  points  which  are  well  covered  by  the  avalan- 
che theories. 

However,  there  are  significant  observations  that  do  not  appear  to 
be  explicable  in  terms  of  the  avalanche  hypothesis.  The  most  notable  of 
these  include: 

1.  Failure  to  account  for  a gaseous  channel  prior  to  the  rise  in 
current  and  fall  in  voltage  indicating  breakdown. 

2.  Failure  to  account  for  the  line  spectra  observed  from  the 
beginning  of  breakdown  conduction  in  thin  film  studies. 

3.  Failure  to  account  for  the  breakdown  configuration  in 
point-plane  geometries  and  in  annealed  single  crystals:  the  channel  in 
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thc''“  situatipng  will  follow  favored  crystallographic  orientations 
rather  than  parallelling  the  electric  field, 

4.  The  nonexistence  of  current  precursors  to  breakdown  renders 
unlikely  any  mechanism  which  requires  the  formation  of  high  current 
densities  prior  to  the  disruption  of  the  dielectric  (25) . 

5.  Failure  to  account  for  the  geometry  of  electrode  damage  in 

breakdowns in  an  avalanche  model,  such  damage  should  be  primarily 

by  Joule  heating.  In  practice,  damage  is  strongly  dependent  on  the 
high  pressure  gases  formed  within  the  dielectric. 

6.  The  chemical  composition  of  the  reaction  products  from  a 
breakdown  channel  is  not  explicable  in  teirms  of  thermal  decomposi- 
tion products,  but  is  expllcab  e on  the  assumption  that  the  material 
is  broken  down  into  its  separate  atoms  and  then  these  atoms 
subsequently  chemically  recombine  (34) . 

7.  There  does  not  seem  to  be  any  direct  experimental  evidence 
showing  the  existence  of  substantial  avalanche  current  in  any  solid 
dielectric  system  so  far  investigated. 

8.  Breakdown  occurs  in  thin  film  systems  at  voltages  too  low  for 
the  process  to  exist.  For  a field  strength  at  breakdown  of  5 MV/cm, 

O 

a sample  having  a thickness  of  500  A will  break  dcrwn  with  an 
applied  voltage  of  25  Volts.  If  the  insulator  in  question  has  a 
band  gap  of  10  eV,  not  more  than  two  ionizations  can  occur  per 
incident  electron,  as  the  electron  will  lose  most  of  its  energy  in 
each  collision.  Thus  avalanche  cannot  directly  cause  the  large 
increase  in  current  required  in  the  avalanche  theories  of  breakdown. 
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It  might  be  veil  to  use  Cooper's  words  to  summarize  the  state  of 
theory  and  experiment  at  this  point;  "Seitz’s  single  avalanche  theory, 
and  also  other  theories  in  which  current  multiplication  occurs  by  elec- 
tron avalanching  and  of  which  the  theory  recently  proposed  by  Forlanl 
and  Mlnnaja  is  an  example,  are  incompatible  with  the  pattern  of  light 
emission"  exhibited  in  Cooper’s  work  (19).  The  article  this  quotation 
is  taken  from  was  published  in  1966,  but  the  views  are  those  which  Cooper 
still  holds  (22).  Others,  notably  Budensteln  (17),  also  expressed  doubts 
about  the  applicability  of  the  avalanche  approach  at  about  this  same  time. 

Gaseous  Model  of  Electric  Breakdown 
A model  of  breakdown  in  solids  has  been  proposed  by  Budensteln  (4,35) 
that  attempts  to  account  for  a broader  range  of  breakdown  phenomena  than 
seems  possible  with  the  avalanche  models.  The  basic  assumption  of  the 
model  is  that,  in  a sufficiently  high  electric  field,  breakdown  in  a 
solid  occurs  when  a highly  conducting  gaseous  channel  is  formed  that 
bridges  the  electrodes.  Thus  the  theory  of  breakdown  must  be  able  to 
explain  how  this  conducting  channel  is  formed.  In  the  formation  of  the 
channel,  the  following  sequence  of  steps  is  hypothesized: 

1.  The  basic  assumption  is  made  that  the  breakdown  starts  when  the 
density  of  excess  charge  in  a local  region  containing  many  atoms  exceeds 
a critical  value.  This  charge  density  may  be  produced  by  any  of  several 
processes,  such  as  field  emission  from  the  electrodes,  ionization  assoc- 
iated with  the  absorption  of  light,  electron  injection  from  a high  energy 
electron  beam,  or  some  other  means. 
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2.  The  effect  gf  this  excess  charge  density  is  the  destruction  of 
chemical  bonds  in  the  local  region.  Vnien  the  bonds  are  broken,  the  prod- 
ucts form  a gaseous  plasma  containing  excited  atoms,  ions,  and  free  elec- 
trons, This  plasma,  to  sustain  Itself  within  the  solid,  must  be  of  high 
density,  high  pressure,  and  high  temperature. 

3.  If  the  free  electrons  within  the  plasma  are  of  sufficient 
density,  so  that  they  produce  at  the  plasma-solid  interface  the  critical 

charge  density,  then  the  process  is  self-sustaining. 

4.  As  the  reactions  proceed,  a gaseous  cavity  is  formed  within  the 
solid.  The  pressures  generated  within  the  gas  can  lead  to  the  formation 
of  cracks  within  the  remaining  solid  material,  extending  the  reaction 

zone. 

5.  When  the  channel  completely  bridges  the  electrodes,  current 
begins,  and,  in  most  circuits,  the  voltage  between  the  electrodes  falls. 


III.  ANALYSIS  OF  BOND  DISRUPTION 


The  present  study  was  performed  to  provide  some  Insight  concerning 
the  nature  of  the  excess  charge*  density  necessary  for  disrupting  chemical 
bonds  and  of  the  energies  of  the  released  particles  after  bond  disruption, 
i.e.,  to  provide  a basis  for  the  first  step  in  the  gaseous  model  of  elec- 
tric breakdovm  described  in  Section  II. 3.  For  the  gaseous  model  to  be 
feasible,  the  excess-charge-induced  reaction  must  be  exothermic. 
(Conversion  of  a solid  to  a gas  Is  normally  an  endothermic  process. 
However,  In  the  normal  situation,  the  system  is  not  subjected  to  an 
external  electric  field  in  the  MV/cm  range.)  In  this  chapter,  the  ener- 
getics of  the  formative  stages  of  breakdown  are  analyzed  in  two  stages. 
The  first  involves  energy  storage  within  the  dielectric,  the  second  the 
solid-to-gas  exothermic  reaction.  For  definiteness,  these  will  be 
explored  in  terras  of  a particular  material,  NaCl,  an  alkali  halide,  in 
a uniform  field  geometry. 

A.  Storage  of  Energy  within  the  Dielectric 
1.  Generation  of  Defect  centers  by  Impact  Ionization 
At  low  applied  voltages,  the  small  leakage  current  through  a bulk 
dielectric  arises  from  mobile  ions  and,  by  means  of  the  hopping  conduc- 
tion mechanism,  or  by  the  Poole-Frenkel  mechanism,  from  electrons.  If 
there  is  no  injection  of  electrons  or  holes  at  the  contacts  (blocking 
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contacts) , only  a transient  current  will  flow  while  space  charge  is  built 
up  at  the  electrodes. 

As  the  applied  voltage  increases,  the  potential  barrier  at  the 
cathode-dielectric  interface  can  be  penetrated  by  Schottky  or 
Fowler-Nordheim  emission  of  electrons  into  the  dielectric,  and  a signif- 
icant injection  current  may  begin.  (This  occurs  at  field  stren'^ohs  less 
than  one-tenth  of  the  breakdown  strength  (36).  Values  for  the  work 
function  at  the  cathode-dielectric  interface  do  not  seem  to  be  available, 
but  a representative  value  is  one-half  of  the  band  gap  of  the 
dielectric  (37)  ). 

As  the  voltage  continues  to  increase,  the  injection  current 
increases,  and  some  injected  electrons  pass  through  the  friction  barrier 
and  attain  energies  sufficient  to  eject  electrons  from  their  bound  sites. 
If  an  electron  of  a Cl  atom  is  ejected,  then  a colour  center  may  be 
created  within  the  dielectric.  The  method  of  production  of  Vj^  centers 
has  been  studied  by  Dienes  (38) , and  their  properties  are  described  by 
Fowler  (39) , Townsend  and  Kelly  (40) , and  Schulman  and  Compton  (41) . 

A Vj^  center  in  an  alkali  halide  consists  of  a molecula’*  ion  occu- 
pying the  lattice  sites  of  two  normal  halide  ions,  with  the  atoms  com- 
prising the  molecule  displaced  from  the  lattice  sites  toward  a common 
center.  The  molecular  ion  is  usually  formed  from  the  ions  originally 
occupying  the  lattice  sites,  A qualitative  diagram  of  this  center  is 
given  in  Fig.  5.  A more  comprehensive  discussion  of  this  colour  center, 
along  with  the  H and  F centers,  is  given  in  Appendix  D.  According  to 
Dienes  (38) , the  center  is  normally  formed  by  ionization  of  one  of  the 
atomic  sites,  followed  by  a combination  of  the  neutral  atom  with  an 
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adjacent  ion  to  form  a molecular  ion.  The  latter  is  more  stable  in 
crystals  at  low  electric  fields  than  the  neutral" ion  configuration  with 
each  atom  remaining  at  its  original  equilibrium  site.  The  stabilization 
energy  of  the  state  in  NaCl  is  about  1.4  eV,  and  the  center  forms  an 
isolated  energy  level  in  the  solid  at  about  6.5  eV  above  the  valence 
band  (38) . Using  the  normal  criterion  for  the  incident  energy  to  be 
three-halves  the  band  gap,  this  would  require  the  incident  electron  to 
have  an  energy  of  about  11.8  eV  to  allow  the  formation  of  this  state 
(the  band  gap  of  NaCl  is  about  7.85  eV) . The  equilibrium  separation  of 
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the  nuclei  in  the  molecular  ion  is  about  2.3  A,  compared  to  2.8  A in 
the  normal  lattice. 

Paracchlnl  (36,42)  has  observed  cathodelumlnescence  at  liquid  nitro- 
gen temperatures  and  below  for  alternating  fields  greater  than  10^  V/cm 
in  several  alkali  halides,  although  not  for  NaCl.  In  the  case  of  a d.c. 
electric  field,  a space  charge  Internally  built  up  in  the  dielectric 
reaches  equilibrium  with  the  applied  field,  and  luminescence  is  not 
observed.  Paracchini  attributes  this  luminescence  to  Vj^  centers  created 
by  impact  ionization  Initiated  by  electrons  injected  from  the  cathode.  A 
field  of  10^  V/cm  is  substantially  below  that  required  for  breakdown. 

The  cathodelumlnescence  is  not  seen  at  higher  temperatures,  as  the  tran- 
sition is  no  longer  radiative  at  higher  temperatures.  Paracchini  finds 
that  the  centers  are  created  predominantly  in  the  cathode  region,  as 
indicated  by  the  light  emission  pattern.  For  each  center  produced,  there 
will  be  a contribution  to  the  space  charge  of  a positive  electronic 
charge,  as  the  center  production  process  has  involved  the  removal  of  an 


Fig.  5-  Sample  configuration  of  Vj^  center  in  alkali  halide 

crystal.  Note  that  the  two  halide  atoms  comprising  the  molecular 

ion  occupy  sites  close  to  those  of  the  original  halide  ions. 

(After  Townsend  and  Kelly) 
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Fig.  6.  The  effects  of  polarization  on  the  formation  of  Vj^ 
centers  is  illustrated  by  this  sequence  of  drawings,  (a)  Shows  a normal 
lattice,  (b)  Shows  the  shift  of  atomic  species  away  from  the  lattice 
points  (indicated  by  dots)  due  to  the  action  of  the  electric  field. 

The  distance  between  halide  ions  is  indicated  by  d.  (c)  Shows  a 
neutralized  halide  atom.  The  distance  indicated  by  d'  is  less  than  d. 
(d)  Shows  the  center  fomed  by  the  neutral  halide  atom. 


32 


electron,  leaving  an  unbalanced  positive  charge.  The  low  temperature 
luminosity  la  due  to  a radiative  recombination  of  an  electron  with  the 
center. 


2.  Properties  of  Defect  Centers 
The  production  process  of  the  centers  may  be  represented  In 
chemical  notation  by  the  expressions 


Cl"  + e"  + EO  

1 

1 
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00 

ti 
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eV 

Cl°  + Cl"  

> CI2 

+ El 

El"l. 4 

eV 

Or,  we  may  describe  the  system  by  the  single  equation 

2 Cl"  + e~  + E'  > Cl"  + 2 e"  E'=  EO-El 

To  attain  the  excitation  energy  at  which  cathodelumlnescence  is  observed 
requires  accelerating  distances  in  excess  of  10  um  in  Paracchinl's  exper- 
iment, and  decreases  to  about  0.1  pm  at  an  applied  field  of  10^  V/cm. 

The  effect  of  polarization  on  the  formation  of  centers  is  Illus- 
trated in  Fig.  6.  The  lattice  positions  of  the  crystal  with  no  applied 
field  are  shown  in  Fig.  6a,  where  the  equilibrium  separation  between  sites 
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Is  (d) , the  separation  being  about  2.6  A.  Polarization  shifts  the  ions 
away  from  their  equilibrium  positions  in  the  neutral  unstressed  crystal, 
as  shown  in  Fig.  6b.  After  impact  ionization  has  occurred  at  a Cl~  site, 
the  adjacent  positive  ions  relax  by  moving  toward  the  neutral  (Fig.  6c). 
This  latter  relaxation  lowers  the  separation  of  the  neutral  and  an  adja- 
cent halide  ion,  thus  favoring  the  formation  of  centers.  The  condi- 

tions needed  for  the  formation  of  centers  exist  in  breakdown  studies. 

The  reverse  reaction  is  also  possible; 

CI2  + e~ 


•>  2 Cl"  + E’  . 
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The  energy  given  off  in  the  reverse  reaction  may  appear  as  lattice  vib- 
rations (phonons),  go  into  the  production  of  paired  H and  F centers,  or, 
at  temperatures  below  SO^K,  may  appear  in  the  form  of  visible  light  (38) . 
The  H center  is  chemically  the  same  as  the  Vj^  center,  and  consists  of  a 
CI2  molecular  ion  which  occupies  a single  Cl  site.  The  F center  is  a 
halide  vacancy;  F centers  serve  as  electron  traps.  The  production  mech- 
anism proposed  by  Pooley  (43)  for  the  paired  H and  F centers  via  the 
recombination  of  an  electron  with  a center  requires  that  the  electron 
be  trapped  by  the  F center. 

The  CI2  state  is  frequently  referred  to  as  a "trapped-hole"  colour 
center,  as  the  now  unbalanced  charge  of  the  Na"^  ion  adjacent  to  the 
center  represents  a positive  charge  localized  in  a region.  This  unbal- 
anced charge  adds  to  the  positive  space  charge  in  the  cathode  region, 
leading  to  the  enhancement  of  the  injection  current,  provided  that  a 
compensating  negative  space  charge  does  not  build  up  also.  The  effect 
of  this  space  charge  is  to  steepen  the  voltage  gradient  in  the  cathode 
region,  while  lowering  the  field  in  the  bulk  of  the  Insulator  relative 
to  its  maximum  value.  The  effect  of  this  electric  field  gradient  is 
two-fold:  first,  the  electrons  injected  into  the  material  can  now  gain 
energy  from  the  field  in  a shorter  distance,  and,  second,  Vj^  centers  are 
generated  in  a smaller  volume  of  the  crystal.  Each  center  is  a region 
of  high  potential  energy;  work  has  been  done  by  the  field  to  supply  this 
energy.  Thus  the  localization  of  the  centers  by  the  electric  field 
gradient  implies  a high  energy  density  as  well  as  a high  space  charge 
density  in  a local  region.  Defects  may  play  a significant  role  in  the 
trapping  of  charge  and  the  localization  of  the  centers.  The  possible 
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role  of  defects  on  the  formation  of  the  initial  space  charge  density  and 
their  subsequent  effects  on  breakdown  will  be  considered  later. 

The  Vj^  center  is  a state  of  high  potential  energy  with  respect  to 
the  normal  lattice.  The  source  of  this  excess  energy  is  the  applied 
electric  field.  centers  are  created  directly  from  the  ions  of  the 
lattice  by  impact  ionization,  and  they  are  annihilated  by  electron 
capture,  with  the  energy  on  annihilation  normally  going  into  phonon  pro- 
duction. There  is,  however,  another  method  of  decay  for  the  center 
on  capture.  This  mechanism,  proposed  by  Pooley  (43),  involves  the  con- 
version of  the  Vjj  center  into  a pair  of  H and  F centers.  The  H center, 
although  chemically  similar  to  the  center,  cannot  easily  decay  back 
to  the  normal  lattice  configuration.  In  order  to  return  to  the  normal 
lattice,  it  is  necessary  for  the  H center  to  recombine  with  an  F center. 

The  mechanism  of  Pooley  requires  that  the  F center  and  the  H center 
with  which  it  is  paired  be  separated  by  several  lattice  spacings  (at 
least  four,  and  normally  more),  so  the  likelihood  of  H-F  recombination 
is  small.  The  H center  retains  the  primary  energy  characteristics  of 
the  parent  center.  Both  the  and  the  H centers  represent  sites  at 
which  energy  is  stored  within  the  lattice,  and  both  of  these  centers  are 
the  primary  result  of  the  build-up  stage  of  breakdown. 

B,  Solid-Gas  Reaction 
1,  Energetics  of  the  Reaction 

According  to  the  proposed  model,  breakdown  begins  with  the  pro- 
duction of  a high  temperature,  high  pressure  gaseous  plasma  within  the 
crystal,  The  rate  of  conversion  of  solid  to  gas  must  be  very  rapid  and 
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highly  exothermic  if  the  gag  is  not  to  condense.  In  the  present  case,  a 
reaction  sequence  that  fulfills  the  requirements  is 


a. 

solid 

E3  

^^^2^gas 

4* 

e~ 

E3=1.4  eV 

b. 

(Na  ) solid 

^^“^gas 

+ 

E4 

E4=5.95eV 

c.. 

(Na*) 

gas 

E5  — 

4* 

e" 

E5=5.9  eV 

The  values  given  for  the  energies  of  the  above  reactions  are  taken  from 
tables  given  In  (44,45),  For  two  atoms  of  Cl  and  two  of  Na  to  go  to  the 
gaseous  state  from  the  form  of  an  Isolated  (or  H)  center,  we  add 
equation  (a)  to  2 times  equation  (b)  to  equation  (c) , and  obtain 

■l-  C«l2>solld  + «l2)gas  + 0<d*)gaa+?d 

where  E^j  ® 4,45  eV  is  the  net  energy  of  the  reaction,  starting  with  the 
CI2  state  within  the  crystal.  Note  that  reaction  (d)  is  exothermic  and 
does  not  alter  the  number  of  electrons  within  the  system.  This  reaction 
sequence  presupposes  clustering  of  the  centers  within  a region  suffi- 
cient to  alter  the  bonding  characteristics  of  the  Na. 

The  reaction  sequence  of  equation  (d)  could  occur  at  any  (or  H) 
center  with  a minute  gaseous  pocket  being  formed.  This  gas  will  be 
rapidly  quenched  unless  additl'^nal  conditions  are  met  which  lead  to 
repetition  of  the  same  or  similar  reactions  in  the  immediate  vicinity 
of  the  first. 


2,  Growth  of  Reaction  Zone 

The  dielectric  is  under  high  electrical  stress  due  to  the  external 
electric  field  when  breakdown  occurs  in  a solid.  When  the  gas  is  formed 
as  described  above,  it  creates  a small  conducting  region  within  the 
dielectric.  The  field  about  the  region  is  distorted  as  shown  in  Fig,  7, 


Fig.  7.  Electric  field  distribution  in  the  neighborhood  of  the 
internal  gas  pocket.  Additional  randomly  positioned  Vj^  centers  in 
the  vicinity  of  the  gas  pocket  are  also  indicated.  Note  particularly 
the  distortion  of  the  field  due  to  the  gas  pocket. 


(t) 

Fig.  8.  Typical  oitxiations  of  Vj^  and  H centers  in  an  alkali 
halide  in  a non-uniform  electric  field  region.  In  (a)  two  Vj^  and 
one  H are  illustrated;  (b)  shows  fovir  Vi.  centers. 
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Thus  a fjLeld  gradient  exists  in  the  region  of  the  gas  pocket,  with  the 
region  of  strong  field  being  centered  about  the  gas  pocket,  Such  a field 
will  exert  a force  on  a dipole  that  is  directed  toward  the  region  of 
enhanced  field.  Suppose  that  there  is  a distribution  of  H (or  Vj^)  centers 
about  the  gaseous  region.  Typical  situations  are  illustrated  in 
Figs.  8a  and  8b.  The  field  will  polarize  both  the  Cl  and  the  Cl"'  atoms 
of  the  (or  H)  center.  The  Cl'^  is  electrostatically  more  strongly  bound 
to  the  local  position  than  the  Cl  because  of  the  surrounding  Na  ions. 

The  Cl  will  be  polarized  by  the  external  field  with  a net  force  on  it 
directed  toward  the  high  field  region.  Thus  the  neutral  changes  its 
affiliation  to  an  adjacent  Cl",  that  is,  the  V-^  (or  H)  center  moves 
toward  the  gaseous  region  under  the  action  of  the  external  field  and  the 
gradient  due  to  the  gas  pocket. 

0 

The  charge  on  a Vj^  center,  a Cl  r-Cl  pair,  can  be  expected  to  move 
from  atom  to  atom  of  the  pair.  When  the  atom  in  the  high  field  region  is 
neutral,  conditions  are  favorable  for  another  jump.  Thus  a mechanism  is 
available  for  motion  of  the  centers  toward  the  plasma  region.  The  reac- 
tion of  Eq.  (d)  can  then  take  place  and  the  plasma  will  grow  in  size, 

3,  Reactions  in  the  Gas 

In  addition  to  the  reactions  of  Eq.Cd),  further  reactions  can  occur 
in  the  gas.  These  include 
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SoSe  of  the  electrons  will  be  Injected  into  the  high  field  region  near 
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the  gas.  These  can  create  additional  Vj^  centers,  the  energy  being 
provided  by  the  field  source. 

C.  Implications  of  the  Gaseous  Model 

The  proposed  model  indicates  that  materials  with  relatively  large 
dielectric  constants,  and  hence  relatively  large  displacements  due  to 
polarization,  should  break  dovm  more  easily.  The  reason  for  this  is  the 
added  ionic  displacement  in  the  initial  step  in  the  formation  of  the  V]^ 
centers.  While  this  will  probably  not  be  the  dominant  effect,  as  changes 
in  the  bonding  structure  between  materials  will  strongly  Influence  the 
breakdown,  there  should  be  a tendency  among  structurally  similar  mater~ 
ials  toward  a reduction  in  breakdown  strength  with  increasing  dielec- 
tric constant.  Such  an  effect  is  observed  experimentally  (1).  Polar- 
ization also  contributes  to  the  stored  energy  in  the  system,  although 
the  contribution  is  normally  small. 

The  model  indicates  also  that  there  should  be  areas  In  which  gas 
production  has  occurred,  but  the  gas  has  been  quenched  before  completion 
of  a channel.  Such  quenching  could  leave  localized  amorphous  (glassy) 
or  otherwise  structurally  altered  regions  within  the  dielectric.  The 
latter  could  include  the  amorphous  clusters  of  free  metal  observed  in 
both  alkali  halides  and  other  materials  (23,24),  and  the  balls  of 
crystalline  silicon  produced  at  the  edges  of  breakdowns  in  thin  film  SiO 
capacitors  (12) . 

Neither  of  these  points  has  been  explained  by  other  theories  of 
breakdown.  Additionally,  the  properties  of  the  channel  are  substan- 
tially in  agreement  with  the  experimental  evidence.  Thus  the  existence 
of  a hot  gas  within  the  channel  will  lead  to  electrode  damage,  and  the 
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light  emiss.ipn  follpw3,n8  breakdown  would  naturally  be  expected  to  be  that 
of  a gaseous  plasma.  The  occurrence  of  treeing  is  also  explicable  in 
terms  of  the  pressures  built  up  In  the  region  of  the  channel  during  the 
formativle  stages  of  the  channel.  The  production  of  gas  in  the  region 
surrounding  the  channel  may  also  aid  in  the  formation  of  cracks. 

The  observed  Increase  in  breakdown  strength  with  the  rate  of  appli- 
cation of  voltage  is  explained  in  the  present  model  by  the  requirement 
that  sufficient  density  of  centers  be  built  up  to  initiate  the  form- 
ation of  the  gas. 

More  will  be  said  later  about  the  dependence  of  the  breakdown  on 
the  presence  of  defects  within  the  crystal,  directional  effects,  ageing, 
and  other  effects. 

The  model  presented  applies  directly  to  breakdown  in  lonically 
bonded  materials.  The  question  arises,  however,  concerning  the  appli- 
cability of  the  model  to  covalent  materials.  One  requirement  of  any 
breakdown  model  is  that  it  be  relatively  insensitive  to  bonding  types 
within  the  material,  as  it  has  been  shown  earlier  that  breakdown  does 
not  depend  strongly  on  the  nature  of  the  chemical  bonds.  It  is  thus 
necessary  to  extend  our  investigation  to  cover  the  case  of  materials 
which  are  covalently  bonded,  and  to  show  that  similar  conditions  may  be 
expected  to  hold  in  this  case  also. 

Cooper  has  recently  shown  that  light  emission  occurs  during  the 
formative  stage  (prior  to  the  breakdown  current)  in  the  cathode  region 
when  the  material  used  for  the  dielectric  was  polymethylmethacrylate 
(PMMA)  (46,47) . This  emission  is  very  low  in  intensity  and  diffuse. 
Cooper  feels  it  is  a necessary  part  of  the  sequence  of  events  leading  to 
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breakdown.  In  the  above  model,  it  is  related  to  the  process  of  building 
up  sufficient  energy  density  through  formation  of  a space  charge  to  make 
the  system  unstable  as  a solid. 

In  the  next  chapter,  forming  the  main  effort  of  this  thesis,  is  a 
calculation  to  estimate  quantitatively  the  degree  of  instability  produced 
by  the  addition  of  excess  charge  to  a molecularly  bonded  system.  The 
problem  has  been  reduced  to  the  simplest  possible  case  so  that  it  could 
be  solved  without  questionable  compromising  assumptions.  By  introducing 
great  simplification,  however,  the  problem  is  removed  from  a complex 
environment  of  the  atom  in  a solid  to  that  .of  a simple  molecular  cluster. 

The  system  to  be  investigated  is  the  ground  state  energy  of  the 
hydrogen  molecular  ion  having  two  excess  electrons.  It  is  knoim  that 
H2  and  H2  molecules  are  stable  in  their  ground  state. 

Flscher-Hjalmars  (48)  has  shown  that  the  ground  state  of  is  stable. 
However,  no  work  has  been  done  previously  on  the  ground  state  energy  of 
the  doubly  negative  hydrogen  molecular  ion.  By  making  such  a calculation 
and  comparing  its  results  with  those  of  the  H2  neutral  and  H2  ion,  an 
estimate  of  the  energy  of  formation  of  the  H2"^  ion  and  the  energy  ' 
release  upon  its  dissociation  can  be  made.  These  have  important  impli- 
cations in  the  context  of  the  gaseous  model  of  breakdown. 


IV.  CALCULATION  OF  THE  GROUND  STATE  ENERGY  FOR 

A.  Formulation  of  the  Problem 

The  ground  state  energy  for  the  double  negative  hydrogen  molecular 
ion  is  given  by  the  solutions  to  the  SchrUdinger  equation 

= Etji 

where  E is  the  energy,  \(»  is  the  wavef unction  appropriate  to  the  state, 

±.2  j 

and  H = - y ^ V. 

The  wavefunction  for  any  molecular  system  may  be  constructed 

from  a suitable  combination  of  atomic  wavefunctions . Then  the  total 

wavefunction  should  correctly  describe  the  system  in  the  limit  of 

large  internuclear  separation,  which  is  just  the  case  of  separate 

atomic  systems.  It  is  desirable  to  choose  a wavefunction  that  can 

be  treated  in  closed  form,  and  one  that  is  sufficiently  pliable 

mathematically  to  be  solvable  without  requiring  extraordinary  mathe- 

-2 

matlcal  procedures.  The  H2  molecular  ion  is  similar  electronically 
to  the  He,  molecule.  Since  Huzinga  (A9)  has  obtained  the  ground 
State  wavefunction  similar  to  his  was  chosen  for  the  present  calculation. 
The  total  wavefunction  is  given  by  the  Slater  determinant 

(1)  t(|)  (l)a(l) 

O 


^„(2)$(2) 

O 


.j>^(3)a(3)  4>^^(4)3(4)| 
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where  a(n)  Is  the  spin-up  and  3(n)  is  the  spin-down  spin  portions  of 

the  total  wavefunction  for  the  nth  electron,  and  the  (}i  and  (}>  are 

g u 

spatial  portions  of  the  total  wavefunction,  and  are  defined  below.  The 
use  of  the  Slater  determinant  takes  care  of  the  requirements  Imposed  on 
the  wavefunction  by  the  Pauli  exclusion  principle,  and  Insures  that 
the  proper  symmetry  of  the  wavefunction  is  maintained.  In  this 
(standard)  shorthand  notation,  the  elements  given  are  the  diagonal 
elements  of  the  determinant. 

The  space  part  of  the  wavefunction  of  the  doubly  negative  hydro- 
gen molecular  ion  is  chosen  as  a superposition  of  two  combinations  of 
atmolc  S orbitals,  A and  A . These  are  defined  by 


(2) 

and 


= N (x 

g '^a 


(3) 


({)  = N (x*  - X^) 
u u a b 


where  Xg(Xg)  and  Xjj(x^)  are  IS  atomic  wavef unctions  centered  about 
nuclei  a and  b respectively.  The  atomic  wavefunctlons  are  defined  by 


(4) 


and 

(5) 


e 


e 


If  a and  y are  different,  then  x x’  and  Xk  Xv* 

Si  Si  D D 

The  constants  N and  N are  normalization  constants,  so  thJt 
g u 


44 


n 


<|)  <J)  dV  dV,  = 1 

g g a b 


and 


^ ^CO 

I I 


dV  dV.  = 1 
a b 


Then 


(6)  N = 


and  N = 


® v^Cl+S)"^  “ 

where  S and  T are  the  overlap  Integrals 


(7)  S 


•CO 

n 

and 

T = 

r 


If  we  define  W = aR  and  x = yR,  N and  N may  be  rewritten  as 

g u 


N 

2~ 

, . -W/,  , . w 

1 + e (1  + w + 

and 

N 

2 ■ 

1 - e"^(l  + X + V) 

g /2 

3 

“ v/2 

3 _ 

Before  proceeding  further,  it  should  be  pointed  out  that  the 
terms  involving  electron  Interactions  with  the  nuclei  involve  the 
interactions  of  a single  particle  with  two  centers.  It  is  natural, 
then,  to  consider  carrying  out  the  necessary  calculations  in  a 
coordinate  system  in  which  there  are  naturally  two  centers  present. 
There  are  two  conic  sections  in  which  two  focal  points  are  present: 
the  ellipse  and  the  hyperbola.  If  we  choose  a coordinate  system  in 
which  confocai  ellipses  form  one  coordinate,  and  hyperbolae,  confocal 
to  the  ellipses,  the  other  coordinate,  then  Fig.  9 might  illustrate 
these  coordinates. 

By  definition,  the  ellipse  is  that  geometrical  figure  fomed  by 
the  loci  of  all  points  having  the  sum  of  the  distances  from  the  point 
to  two  foci  as  a constant.  Similarly,  an  hyperbola  is  the  loci  of 


f 

constant 


^ constant 


Fig.  9*  Confocal  coordinate  system  used  for  evaluating  two-center 
integrals.  The  third  coordinate  is  the  rotation  about  an  axis  joining 
the  two  foci,  which  are  indicated  by  dots  in  the  figure. 
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w = otR 


X = yR 


TABLE  1 


NOTATION  USED  FOR  CONFOCAL  INTEGRALS 


screening  parameter  for  inner  electrons 
screening  parameter  for  outer  elettrons 
spin-up  portion  of  the  wavefunction 
spin-down  portion  of  the  wavefunction 
internuclear  separation 


distance  of  electron  from  nth  nucleus 


Useful  confocal  coordinate  relationships 


X = r(X^-l)(l-V^l 

coord  2 _ 


= XyiR 

^coord  2 


^coord’  ^coord’  ^coord  ^ 
rectangular  coordinates 


y2  4 ^ 

d2/,2  2.  ^ 

R (X  -y  ) 


L]  + 3y  3y]  9^  [(^2_^^'J^_^2)  ay] 


, - Rg  + y) 

’^a  “ 2 


^ R(X  - y) 
b 2 
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all  points  such  that  the  difference  of  the  distances  to  the  foci  is 
a constant.  In  Fig.  9 this  is  expressed  by,  for  point  P, 


r 

a 


r 

a 


= e 


= h 


along  the  ellipse 
along  the  hyperbola 


These  constants  are  related  to  the  coordinate  system  by  the  separation 
of  the  foci.  Thus 


V = 


r + r, 
a b 


£ 

R 


and 


X = 


h 

R 


We  also  have  the  Inverse  transformation 


Ta  = (|)  R (X  + y)  r^  = (i)  R (X  - y) 

which  will  allow  the  conversion  of  the  Hamiltonian  into  confocal 
(elliptical)  coordinates.  The  remaining  coordinate,  i{)  , represents 
the  rotation  about  the  axis  of  the  foci. 

Before  proceeding  further,  it  is  also  necessary  to  express  the 
Laplacian  operator  in  confocal  coordinates.  In  rectangular  coordinates, 
this  operator  is  given  by 


By  expressing  X,  p,  and  (|)  in  terms  of  x,  y,  and  z,  and  using  the  chain 
rule  to  obtain  the  necessary  partial  derivatives,  the  result  may  be 
obtained  (51) 
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n2,^2  2. 

R (A  -u  ) 


(A^-1)  -^  + 2x  + a-/)  ^ 

3A  ° 8y 


A*--! 


2'  2 
1-y  3({> 


Transformation  to  this  coordinate  system  will  allow  the  greatest 

simplification  possible  in  the  individual  energy  terms. 

In  calculating  the  energy,  it  is  cumbersome  to  carry  along  a 

large  number  of  (unit-determined)  constants.  Thus  it  is  convenient 

to  adopt  the  so-called  atomic  system  of  units.  In  this  convention, 

masses  are  measured  in  terms  of  the  mass  of  an  electron,  m , and 

e 

charges  in  terms  of  the  charge  of  a proton  (equivalent  in  magnitude 
to  that  on  an  electron,  but  differing  in  sign).  Angular  momentum 
in  this  system  is  measured  in  units  ofii.  One  can  thus  express  the 
necessary  conversion  factors  from  this  system  of  units  to  MKS  units  as: 


-31 

1 at.  mass  unit  = m “ 9.1  x 10  kg 
e 

-19 

1 at.  charge  unit  = e = 1.67  x 10  Coulomb 

1 -34 

1 at.  ang.  mom.  unit  =11  = 6.67  x lO  J-sec. 

2-n 

We  still  have  at  our  disposal  one  unit  which  may  be  arbitrarily 
chosen  to  determine  all  other  units.  (In  the  MKS  or  SI  units,  the 
meter,  the  kilogram,  the  second,  and  either  the  Coulomb  or  the  Ampere 
(=  Coulomb/sec)  are  chosen  as  the  basic  undefined  units.  Four  units 
have  been  found  sufficient  to  determine  all  other  units  of  physical 
measurements.)  It  is  convenient  to  consider  for  the  last  unit  a 
length;  conventionally,  the  length  unit  chosen  is  the  radius  of  the 
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Innermost  electron  shell  of  hydrogen,  as  given  by  the  calculation  of 
Bohr.  This  unit,  designated  a^,  is  knovm  as  a Bohr. 

a = 1 Bohr  = 0.5292  A = 5.292  x lO"^^  m. 
o 

As  the  unit  of  energy,  two  choices  are  frequently  seen  in  the  literature. 

The  potential  energy  of  a single  electron  attracted  by  a single  proton 

2 

(i.e.,  a hydrogen  atom)  is  given  by  -e  /r,  where  r is  the  distance 
between  the  two  particles.  If  we  take  the  energy  at  a unit  distance 
in  atomic  units,  we  obtain 

2 

|-e  /a^l  = 1 Hartree  = 1 At.  energy  unit. 

This  energy  is  twice  the  energy  of  the  hydrogen  atomic  ground  state. 

Thus  the  energy  may  be  related  to  other  units  by 

1 Hartree  = 27.21  eV. 

The  other  choice  of  units  for  energy  frequently  made  is  the  Rydberg, 
where 

1 Rydberg  = (1/2)  Hartree  = 13.605  eV. 

The  latter  unit,  better  adapted  to  calculations  on  hydrogen,  will  be 
used  in  the  calculations  to  follow.  The  choice  of  these  units  will 
effect  a great  simplification  in  the  form  of  the  equations  without 
any  loss  of  physical  content. 
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The  energy  of  the  system  will  be  given  by 


fto  fCO  ^ fO> 


E = 


i|)*Hi|)dVj^dV2dV^dV^, 


.00— oo~oo~co 


where  H is  the  Hamiltonian  operator  appropriate  to  the  system.  For 
the  hydrogen  molecular  ion,  H is  given  by 


(8) 


2 _ _2_  _ _2 2 2 2_ 

’^al  ^a2  ’^aA  ^bl  ’^b2 


where  R is  the  internuclear  separation.  The 
of  the  ith  electron  from  the  jth  electron,  and  r^^^ 
of  the  ith  electron  from  nucleus  a (nucleus  b) . 
rewritten  in  a more  compact  form  as 


are  the  distances 
(r,  .)  are  the  distances 

Dl 

Equation  (8)  may  be 


(8.1) 


4 

H=  I 


i=l 


where  the  sum  over  i is  a sum  over  the  electrons  present  in  the  system. 
After  expanding  the  wavefunction,  grouping,  and  cancelling  terms  which 
are  zero  (see  Appendix  A) , the  energy  can  be  expressed  in  terms  of 
the  original  set  of  basis  functions  as 


(9) 


£jd*V 


£j,d‘v} 


51 


where 


(10)  - <t>g(l)<^y(2)<^^(3)^^(4)  -<  <i.g(l)<|.^^(2)4.g(3).J.^(4) 


+ <|.y(l)4.g(2).J.g(3)<t.^(4)  - <t.y(l)<i.g(2)<t^(3).J.g(4) 


is  the  spatial  portion  of  the  wavefunction  In  teinns  of  the  original 
basis  set,  and  f^^  is  the  complex  conjugate  of  f.  Equation  (10)  may 
be  further  reduced  by  noting  that  f^^*  = f^^,  and  that 


O ^ 

<j.g(l)<{.g<l)dV^  = J <j,^(i)4,^(l)dV^  = 1. 


It  is  then  easily  oeen  that  the  terms  of  the  energy  may  be  written 
in  a form  analogous  to  that  of  Huzlnaga  (49),  For  ease  of  comparison, 
his  notation  is  used  here. 


(11)  E=2H+2H+J  +J  +4J  -2K  +~ 

g u gg  uu  gu  gu  R 

where 

(12)  H = [ ♦ (1)(-  7^  - (i)dv 

g J s -L  g J- 

•■CO 

(13)  H^  = | +^(1,  (.  ,2  - ^ (IMVj 

al  bl 

mOO 

—00— CO 

o ^ 

J <l>^(l)<|>^(2)(^)4.^(l).j)^(2)dV3^dV2 


(14)  J 


(15)  J 


(16)  J 


gu 


= I”  <|>g(l)<^y(2)(-^)<j.g(l)4.^(2)dV^dV2 
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(17) 


K =r 

gu  J 
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4.g (1) <1.^ (2)  (^) ^^(1) (2) dV^dV^ 


These  integrals  are  most  conveniently  handled  by  separating  them 
further  into  the  atomic  wavefunctions  comprising  them.  Note  that 
the  Hg  and  integrals  refer  to  interactions  of  the  dectrons  with 
the  nuclei,  the  J integrals  (Jgg»  ^uu*  “^gu^  Interactions  be- 
tween "equivalent"  electrons,  and  the  Kg^  to  interactions  between 
"unlike"  electrons. 


B.  Evaluation  of  Terms 

Each  term  may  be  expressed  in  terms  of  the  atomic  wavefunctlon 
basis  set  and  then  separated  into  simpler  integrals  over  atomic 
states.  If  we  consider  Hg  first. 


(12) 


♦ (1)  (-  'i  - 

® al  bl  ® 


"g  = "U  '’‘a  ^ 


(12.1)  H = N {2  , 
S S J 


9 2 9 

Xa<-  'l  - 77  - 

al  bl 
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The  first  integral  is 


(18) 


X (-  ^ 

1 'al 


-^)x  dv, 

^bl  ^ ^ 
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2 
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and  the  second  is 
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(19) 


X (-  ? ^)X,dV,  = - ae"''r(2+2w)  (2-a)  + (l+w+^)a| 

^ ^ "al  "bl  ^ ^ ^ 


Substitution  of  Eqs.  (18)  and  (19)  into  Eq.  (12.1)  gives 


(20)  H = 2oV{a  - 2 (1+-)  (l-e“^'')  - e"'*  [(I+m)  (A-a)  + ^ 

g g w L J 


In  a like  manner,  H is  evaluated  as 
u 


(21)  (y  - 2(l+^)(l-e"^’^)  + e”’^|~(l+x)(4-Y) 

The  integrals  Involving  the  electron-electron  interaction  terras 
are  naturally  more  difficult.  Each  will  be  taken  in  turn. 


(14.1)  J = r { 
gg  g 


+ 2 


M I Xa(2)Xb(2)(— )dV^dV 


"12 


+ 4 


X3(l)X,(2)Xb(l)Xb(2)(^)dVidV2  + J J X3(1)X,(2)  (— )dVj^dV2 


2./--  \ . 2 / 2 
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+ 2 
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Xb(l)X3(2)Xb(2)(^)dV^<iV2+  X3(2)Xb(l)(7-)dVidV2 


2/-n\  2/,\  / 2 


12  - - ; / - - 12 
■*■  _0O— CO 

2/-\  2,-\  / 2 


X,tt)X,a)x^(2)(^)gVjdV2  + ] J X^(l)X^(«(^)dV^dVj} 


If  like  terms  are  equated,  J becomes 

gg 
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(^,aVidV2  + 2 I" 


Xa<W>!b<«>'a<«Xb(2><J^>JVidV2 


n X^l)X^2)(^^)■iVldV2} 


This  may  be  written 


(14.2)  J = 2N  {J  , + 4J  - + 2J  , + J ,} 
gg  g ggl  gg2  gg3  gg4' 


Equation  (14.2)  is  evaluated  term-by-term  as  follows. 


fCO 


(14.2.1)  J 


ggl 


Next, 
(14.2.2)  J 


gg2 


X^l)xJ(2)(^)dV^dV2  = f 


X^l)X,(2)Xb(2)  (^^)dV,dV2 


The  third  term  is 


(14.2.3)  - [”[”  x,(l)x,C2)x,(l)x,(2)(^)dV3dVj 


2a  r -2w,  25  . 23w  . 2 , w"^..  . 6 ro2,p.(>_v 

•^ggs  = T ® T^ir  w (C+£nw) 


+ (S')^E^(-4w)  - 2SS'E.c(-2w) 


> 


X^l)X,(2)XtC2) 


The  first  terra 


where 
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S = e"”(l  + w + y) 


o1  +W/1  J-  w X 

S - e (1  - w + -J-) 


C = 0.57722....  (Euler's  Constant) 

and  t.'!  (-x)  is  the  exponential  integral  function  (52) . The  last 
term  is 


(IA.2.4)  J 


88^  ~ J 


X^(l)X?(2)(-^)dV  dV 
a b 1 2 


88^ 


[2 

= a 

[w 


-2w.2  . 11  ^ 
e (—  + -7-  + 
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In  an  identical  manner,  may  be  dete' mined  , From  this  point  on, 
limits  will  be  omitted  as  integrals;  the  limits  of  integration  are 
the  ranges  of  the  variables  of  integration 


(15)  J 


uu 


Equation  (15)  may  be  written  in  terms  of  the  atomic  orbitals  as 


(15.1)  J = 2N^' { 
uu  u 


- 4 


I x;'(l>x;(2)x;(2)(^)dV^dV2  + 2 II  x;tt)x;(2)xi<l)x'(2)(^)dV,dV 


x;^i)x;^2)(^)dv^dv^) 


(See  Appendix  A for  details  of  the  expansion.) 
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Equation  (15.1)  may  be  written 


(15.2)  J = 2N^  {J  , - AJ  _ + 2J  , + J ,} 

uu  u uul  uu2  uu3  >m4 


where  J ,,  J o>  ^ o>  and  J , are  defined  and  evaluated  below, 
uul  uu2  uu3  uuA 
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(15.2.2)  J 
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(15.2.3)  J 
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where 


T = e ^(1  + X + y) 


T’  = e^(l  - X + ^) 
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Next,  consider  the  firs^t  of  the  interaction  terras  between  the 


two  types  of  electrons,  J . 

gu 
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+ II  Xjl)x;^(2)Xb(l)(^^>dV^dV2  - II  xJ(l)x;(2)xi(2)(^)dV^dV2 

- j X,(l)X,(l)x;(2)x;<2)(^)dV^dV2  ■(•  II  X^l)x;'(2)(^^)dV^dV2 

+ II  X^(l)Xb(»x;^2)(^)dVjdV2  - II  X,(l)Xb(l)x;(2)x;(2)(^)dV^dV. 

- ] X^l)x;<2)xi(2)  (^)dV^dV2  - II  xJ<l)x;(2)xi(2)(^)dV^dV2 

- ] x.(l)Xt(l)x;(2)xi(2)<:^)dV^dV2  + II  x^Uxi'(2)(^^)dV,^dV, 

+ II  X,(l)X^(l)x;^2)(^)dVjdV2) 


If  in  Eq.  (16.1)  is  expanded  in  terms  of  the  basis  set,  tnay 
be  written 


(16.2)  - 2»V  (||  x2(1)x;^(2)  (v^ldV^dV^ 

+ xf(l)x;^2)(-;:^)dV  dV  - 2 [f  X„a)Xt(l)x'<2)x;(2)(-^)dV  dV 

jj  a b ^ JJ  a b a D r^2  ^ 

- 2 |j|  X^l)x;(2)xi(2)(^)dV^dV^  - ||x;^(2)x^(l)x,(2)<^^)dV.,dV, 


Then,  using  a symbol  to  represent  each  integral,  J may  be  written 

O 
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gu 


2 2- 
rN 
g u 


(16.3)  J_.  = 2N^Nt  {J_.,  + - 2J_..,  - 2(j  ,-J 


gul  ‘'gu2 


gu3 


gu4  gu5 


]> 


These  terms  are  defined  and  evaluated  below. 


(16.3.1)  J 


gul 


xJ(i)x:^(2)(~)dV  dV  = 2y 
12  ^ 


2 2-1 
1 _ y - - 

(a+y)^  2(a47)^ 


(16.3.2)  = JhJ<">X;'(2)(^)dV^dV2 


j = r_  r Q _ -2x.  .IjoL-)  _ _L 

gu2  R ^^3^  3 

® Y Y 


(x  - 1)  + e ^^(x  + 1) 


(e  ‘ - e ) (1  + w + x)  e ^(1  - w 4-  x)  - e ^(1  + w - x) 

2 ^ 2 
(a  - y)  (a  + y)  (a  + y)  (ct  - y) 


-2x  -2w  ^ 

ale  -e  |(14-wd-x) 

2 

(a  - y)  (a  + y) 


(a+Y)  (a-Y)' 


t 


^ -2x,  .2„2 


(16.3.3) 


gu3 


(a-fY)^R^  -w-x 
p6w(l+w+x+  3 ^ 

(a  + Y)^R 


e~^(14x) 
y2  _ 


+ e ^~^(l+w+x)  ^ 

(a+Y)^ 


2(a-fY)  Vy 
we  ^ ^(l+w4x+ 3 

(a+Y)\ 


6S'TE-t(-2w) 

2 

Y X 


+ 


6S'T'E.t(-2(w+x)) 

(6S’e~^''"^) 

Y^(1+2w+2x)  (1-x) 

2 

Y X 

2 

_ Y X 

_ 2(a+Y)2 

(16.3.4) 


a (a  + y) 


(a  + y)' 


L e~^(2  - 2w  + w^)  - e~^^~^(2  + 2w  + vh  , 

2ot  (a  + y) 

V = 11  x;^2)X,(l)x,tt)(j^)dVjdV2 


I - — f2e  ^(1  + w)R  (e  ^ (1  + w + x) 

gu5  “ R ^ 2 ■ / ^ n2 

a Y(o  + y) 


-2x-w,.  , , -w. 

e (1  + x)  + e (x  - 1 


-w  -2x-v. , , , 

e -e  )(l  + w + x + 


Y (a  + y) 


(a  + y)' 


**w  >/%  n • 2n  ”2x*^w  / — to  I 

^ e (2  " 2x  + X ) - e (2  + 2x  + x ) 


2y  (a  + y) 


The  final  set  of  Interaction  terms  which  must  be  considered  are 

in  the  K term:  it  is  the  most  troublesome  of  all  the  terms  in  the 

gu 

evaluation  of  the  energy. 
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r 


(17)  V = 


4>g  (1)  <1.^  (2)  (^)  4-^  (1)  (t-g  (2)  dV^dV2 


It  is  easy  to  see  why  this  integral  is  difficult;  it  involves  the 
interchange  of  an  electron  in  one  orbital  with  that  in  another 
orbital.  The  expansion  of  in  terms  of  the  basis  set  is 

<17-»  V ° ‘II 

- If  X^»)x;(l)x;(2)Xt(2)(^)dVj^dV2  - 2 II  Xjl)x;(l>xj2)x;(2)  (^)dV^dV2 

- II  X^(l)x;<l)Xb(2)x;(2)(^)dV^dV2  + ||  X Jl)xi(l)Xj(2)x;(2)  (^)dVjdV2 


+ 2 


X,(l)x;(l)x;(2)Xi<2)(-;^)dV,dV2) 


12 


Equation 

(17.1)  ' 

(17.2) 

K = 2 
gu 

where 

(17.3.1) 

\ul 

(17.3.2) 

’^gu2 

(17.3.3) 

K , 
gu3 

(17.3.4) 

K / 

gu4 

fV  (K  - K - 2[k  . - K J - K ) 


'xJl)x;(l)X3(2)x;(2)(^^)dVidV2 


■12 


X,(l)x;(l)X,(2)x;(2)(^^)dV^dV2 


r. . 1i  mil  »*i^‘*^riiglil*il*i>iiiiffi<WT'l  iilMMl  ■miii'Mni  ~ irtiiirt1ir<iiiTi  * hfHiiii  iinliTi  ' iumi  i' 


61 


(17.3.5)  K 


ff 


gu5 


X^tt)x;(l)x,(2)x;(2)(^)dV^dV2 


x;  (1)  X^,  (1)  x;  (2)  Xj,  ( 2)  (7^)  dV^dV2 


The  term  will  have  to  be  handled  separately — It  is  the 
only  one  of  the  Integrals  Involved  which  cannot  be  solved  in  closed 
form.  (However  difficult  the  other  integrals  may  appear,  they  each 
possess  "closed"  form  solutions  in  terms  (at  worst)  of  transcendental 
functions.)  The  values  of  the  other  terms  are: 


3 3 


(17.3.1.1)  K . Majc 


33  3 3 

(17.3.2.1)  K , = { 16-e“”~^  fie+ll  (w+x)+3  (a+y)  - 1 > 


(a+Y)  R 

3.. 3 „ r -X/  „ . 2 ..2s 


(17.3.3.1)  K 


8a 


II-  f8re~^(-2Y+(a  -Y^)R)  +_2Ye~^ 

^ \ O O 


/ 2 2s2_ 

(a  - Y ) R 


, -X  -2w-Xs  on  /I  \2„2 

— 2^  (|  + I (w+x)  + ) 

a(a  + y) 


_ e~^^'^(wfl)  + e~^(w-l)  ^ e ^(w^-2w+2)  - (w^-f 2w+2)  ^ 

a^(a  + y)R  4a\ 


fT7  A K ■ 8a\^  ■ f8e'”(-2a+(Y^-a^)R)  + 16  e 

(17.3.4.1)  K - 3 t 2 2.2^ 

* (a+Y)  (Y  - a ) R 


-x 


***W  “W'“2X\  n A / 1 \ ^ O 

<?..  e.  ).  4 I ^ ^2) 

Y(a+Y)  ^ 


— c.*  tiaa-.urc.  . g 


-.'..-xSXiA  .K,„  W ' ai  V 


, i ' < 
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e ^ ^^(x+1)  + e ^(x-1) 


Y (a  + y)R 


e (x^-2x+2)  - e ^^(x^+2x4-2) 
3 

4y  R 


The  remaining  integrals  (designated  as  are  best  handled 

by  returning  to  a combined  form,  so  that  the  integrand  will  possess 
symmetry  properties  making  the  solution  a bit  easier  to  obtain.  In 
this  form, 


(17.3.5.1) 


K 


gu5 


(X^(1)x;(1)(7^))(X^(2)x;(2)  - x;(2)x^(2))dV^dV2 


The  evaluation  of  this  integral  will  be  treated  in  detail,  as  this 
integral  presents  the  greatest  difficulty  of  those  encountered. 
Consider  first  the  integral  over  the  separation,  namely 


(22) 

As  the  integral  (Eq.  (22))  involves  (•^)  , where  is  the 

12 

separation  of  the  two  electrons,  it  is  necessary  to  express  (- — ) in 

12 

terms  of  confocal  (elliptic)  coordinates.  This  expansion,  normally 
termed  the  Neumann  expansion,  is 


(23)  (z^)  = I I I (-l)"'(2k+l) 


'^12  ^ k=o  m=-k 


Rk- 

jn 

_(k+ 

m 

)L 

n 2 


X P, 


m 


X(a)  X(b)  P I 


m 


m 


im((}.^-(^2) 


where  (X^,  (}i.jj^)  are  the  coordinates  of  the  ith  electron,  and 

rx(a)l  is  the  lesser  and  [X(b)]  the  greater  of  X^^  and  X2.  The 


m 
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and  Qj^  are  the  associated  Legendre  functions  of  the  first  and 
second  kinds,  respectively  (60). 

The  integral  of  Kq.  (22)  over  <{)j^  and  412  may  be  carried  out 
Immediately  by  noting  that 


2tI  flTJ 


* fi  I 

J l^t=0  m=-k 


k im((t)j^-(}.2) 


] '*1' 


d(j),  = 0,  m 0 


since 


27T 


-im(|) 
e ^diji 


= i PI"  - ^ = » 


Thus  4 becomes  a nugatory  coordinate,  and  m is  zero  for  a non- 
vanishing integral.  For  m = 0, 


fZ-n 


d4, 


2Tr 


^*^2^^^!’  ^2’  ^1’  ^2^  * ^2*  ^2^ 


Since  m - 0,  there  is  no  azimuthal  dependence,  so  the  separation 
between  the  two  electrons  can  be  expanded  using  the  expansion 


(l)x'  (1)  = g-(a+Y)RX/2  ^-(a-Y)Rli/2 


'a'"'b 


The  integral  of  Eq.  (22)  can  be  written 


(22.1)  = (-St^)(^)(2^)(i)  I (2k+l) 

k=0 


(A^-u^)P,[x(a)] 
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in  confocal  coordinates,  where  the  integration  over  azimuthal 
dependence  has  already  been  carried  out.  This  may  be  rewritten  as 


k=0 


(22.2)  I (2k+l)  (x2-p2)Pj^[x(a)]Qj^[X(b)l 


-1  1 


X ^-(a-Y)RWi/2 


-(g-Y) Em 

It  is  the  presence  of  the  exponential  factor  e ^ which 
makes  this  integral  difficult:  this  factor  prevents  simplification 
by  use  of  orthogonality  conditions  of  the  Legendre  functions . 

To  proceed  further,  make  the  expansion 


(24) 


-(a-Y)Ru/2  _ 


= I 

j=0 


(-(a-Y)f)^(fi 


r)  = I 

j=0 


(a-ir)R 


j uis 
j! 


The  , however,  may  easily  be  replaced  in  this  series  by  Legendre 
polynomials,  which  will  be  truncated  after  only  a few  terms.  This 
allows  the  use  of  the  orthogonality  properties  of  the  Legendre 
polynomials . 

It  will  be  noted  at  this  point  that  a "straightforward" 
expansion  in  orthonormal  functions  (i.e.,  in  Legendre  polynomials) 
could  have  been  used  in  preference  to  the  treatment  adopted.  The 
drawback  to  this  procedure  lies  in  the  problem  of  representing  the 
exponential  adequately  with  only  a few  terms;  approximately  four 
times  the  number  of  terms  would  have  been  required  for  acceptable 
accuracy  with  the  "straightforward"  approach  as  were  actually  used 
with  the  approach  given  here. 
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The  original  series  will  thus  be  replaced  by  the  new  series 


E A.p,(u) 

j=0  J J 


where  the  coefficients  will  be  evaluated  by  direct  comparison 
of  the  power  series  expansion  with  the  expansion  in  Legendre 
polynomials  when  the  power  series  is  expanded  in  terms  of  Legendre 
polynomials. 

Substitution  of  the  expansion  of  Eq.  (24)  into  the  integral 

gives 


OO  f \ fCO 

(22.3)  I I (2k+l)  (xJ-upP^[x(a)]Q^[X(b)] 

X P^<bi)P^(W2)  I A P (p^dX^dp, 


j=0 


,2  2. 


The  dependence  of  (A^-Pj^)  involves  only  ^0(^2)  ^2^^1^*  Hence, 
the  integrand  will  be  expanded  to  obtain  the  dependence.  Then 
the  summation  over  k involves  only  a finite  number  of  terms.  Thus, 


2 2 

the  expansion  of 


(25)  (A^  - p^)  = (A^  - - I 


2A 

(26)  (xj  - p^)  I^A^Pj  (Pj)  - Ay  - i)P„(Pj)  - y 


Ik 
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TABLE  2 

FIRST  FEW  LEGENDRE  POLYNOMIALS 

n P (x) 

n 

1 X 

2 (\)  (3x^  - 1) 

3 (|)(5x^  - 3x) 

A (|) (35x^  - 30x^  + 3) 

5 (|) (63x^  - 70x^  + 15x) 


Q„(*) 

: i)  - 

1 

(^)  (3x^  - 

l)£n(- 

X 

- r '•2^ 

X 

(f)  (Sx^  - 

3x)£n(; 

X + 1 5 

X - 4 

(f) 

(^)  (35x^ 

- 30x^ 

+ 3)£n(^ 

(f)x3  4. 

( 0 A ) X 

<-T 

(63x^ 

- 70x^ 

+ 15x)£n(- 

X 

- (f)-' 

+ (f)x2 
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2A 


• * « • 


= I (^1  - - I I 

j “ - J " 

(For  convenience,  Table  2 contains  the  first  few  Legendre  polynomials.) 

Next,  the  products  of  Pj(p)P2(y)  are  expanded  and  the  series 
truncated  for  j = 3. 

^1^^1^^2^*^1^  = -5  + J 


These  results  reduce  the  Integral  of  Eq.  (22.3)  to  the  form 


00  /OO 

(22.4)  = >^5^  V?  I (2k+l)  J Pj^[X(a)]Qj^[x(b)]Pj^(y2) 


-1  1 


X e-<“^>‘‘Xr/2  j-„2  . - I I BjP^(„j)]  ax^d„^ 


where  the  coefficients  are  obtained  by  a direct  summation  of  the 
coefficients  of  P^  (pj^)  in  the  expansion.  The  orthogonality  of  the 
Legendre  polynomials  is  expressed  by  the  equation 


(27) 


P^(2>Pj(z)dt  . (^) 


-1 


This  allows  Eq.  (22.4)  to  be  written 
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(22.5)  I I <2M)(^)  6 

k,2. 

X I”  P5.[X(a)]Qj.[X(b)]Pj,(y2)e"^“‘^^^»'l/2  [(X^  - dX^ 

1 

Due  to  the  involved  nature  of  the  integrals  in  the  series 
expansion,  it  was  found  necessary  to  truncate  the  series  after 
three  terms.  The  entire  contribution  to  the  total  energy  amounts 
to  only  about  five  percent  of  the  total  energy  (the  term  is  positive, 
and  thus  is  a weakening  of  the  bonding) , so  that  its  truncation  did 
not  lead  to  a significant  error  in  the  total  energy.  For  an  inter- 
nuclear  separation  of  one  Bohr,  this  truncation  resulted  in  an  error 
of  about  eight  percent  in  the  value  of  the  term,  as  compared  with 
the  results  given  for  a similar  term  by  Huzinaga  (51)  (see  Appendix  F 
for  the  numerical  values  in  the  comparison).  This  accuracy  was  felt 
to  be  sufficient,  as  the  addition  of  a further  term  to  the  series 
would  have  greatly  lengthened  the  calculation  and  required  an  un- 
reasonably long  delay  in  the  completion  of  the  work. 

At  this  point,  it  is  convenient  to  give  the  coefficients 

resulting  from  the  expansion  of  the  exponential  and  the  product 

2 2 

of  this  expansion  with  the  (X^^-Pj^)  term.  Coefficients  are  given  in 
Table  3 for  an  expansion  of  the  exponential  through  the  fifth  power. 
Only  the  terras  with  j<3  were  used  in  the  calculations. 

After  the  integration  over  p^, 

(22.6)  = 2 I {A^  | P^Cx (a)]Qj^[x (b)] P^(P2) 

^ 1 
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X - i)dX^  - I” P^Cx(a)]Qj^[X(b)]P^(P2) 

1 


In  Eq.  (22.6),  X(a)  is  the  lesser  and  X(b)  the  greater  of  X^^,  X2  • 
Thus  the  integration  over  X^^  involves  integration  over  the  P^[x(a)^ 
for  values  of  X^^  less  than  X2»  with  Qj^(^2^  Independent  of  X^^,  and 
integration  over  Qj^[X(b)]]  for  values  of  X^^  greater  than  X2,  with 
^k^^2^  constant.  It  is  therefore  necessary  to  take  the  integral 
separately  over  the  range  X^^  < X2  and  X^  > X2» 

Equation  (22.6)  will  be  expanded  term  by  term.  Let  m = (a+Y)R 
for  convenience.  Let  us  consider  first  the  k=0  term: 


(28)  I23  = 


Po(w) 


X - ^)dX^} 


The  first  of  these  is  evaluated  in  Table  B-1  of  Appendix  B.  The 
second  integral  does  not  contain  the  value  Xj^=l,  so  there  are  no 
singularities  to  be  avoided. 

Using  the  values  tabulated  in  Table  4,  the  integral  of  Eq.  (28) 
may  be  written 


(28.1) 


^28  ^o^^^  ^2^'^^X-l^ 


(in/2)^ 


(1+Y  + 


2 

m 

8 


70 


4 


li 

3 

'-tn 


.2,2 


(1  + "I.  + - (1  - f + n(i+i)) 


+ (1  + I + f(X-l))  - f 


_2_ 

3m 


-mX/2p ; /X+1n 


- f(X+l))  + e-^^hU-  f(X-l)) 


Equation  (28.1)  may  be  rearranged  as 


(28.2)  ^28  " “3  ^2'^’^^^^  e"™^^(^+m+2)j  - + 2 

m “ 


+ + ? - T7)E-i(-  ?<>>-«>  - 


2 12 


The  comparable  B term  for  k=0  is 


(29)  Ij,  = Pg(yj) 


, . e""^^Ei(-  f(X+l))  + e"'^Ei(-  S(x-1») 

III  X-1  ^ 


2]-  ; / m / 


The  term  of  Eq.  (22.5)  is 


(30)  I 


30 


Pj^[X(a)JQ3^[x(b)JPj^(y2)e"’"^^^^(^l  " 


Pl(y2)  - •j)dXj^  - X. 


-mXi/2(x2  _ 


X c»  X,+l 
^ ' X 


fftnTi  * fr 


-^^-•'--'-tii/a«rMiifkiBri‘ir.finrTi-ai-t,^ 





71 


table  3 

EXPANSION  COEFFICIENTS 


A. 

3 

(a-Y)^R^  , (ct-Y)V 

>■  -24 ■*■  1940 

(a-Y)R  (g-Y)  V _ Ia-Y)V 

- 2 80  63(3840) 

(a"Y)^R^  + 

12  679 

(a-Y)V  (a-Y)V_ 

120  5040 

4 4 

2(g-Y)V 

3395 

5 5 
(g-Y)  R 

30240 


B. 

3 


15 


187A, 


4A, 


15 

4A 


15 


_A  + th.  + ^ 

21  21  3 


130A, 

99 

40A 


56A. 


231  35 

20A^  2OA2 

TiT 


5 


lOA^ 

33 

42A5 


0 


143 


Nils 


74 


16 


m 


.2,2 


- (l-B+t) 


(-f(X+l))  + (1+f +4)  !<’■-«>] 


96 


m 


ro 


2 3 

mX  , m^X^  , m^X^v  -mX/2p  A+ls  ^i_®4.n! — — 


(1+- 


8 48 


2 3 


24(32) 
5 


m 


j,,  , mX  , m^X^  , m^X^  l 


234  /o  „ 

. /-.  . n>  1 ni  j.*"  j-JS CX-l'^'i 


ij.L  uiaiiuipiLjii  ijiijjijRiimpi..iij|M«piiniiii^ 


m 
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Of  the  integrals  in  Eq.  (30),  the  first  two  are  given  in  Table  B-1, 
and  the  third  in  Table  4.  Then  the  integral  of  Eq.  (30)  becomes 


(30.1)  I3Q  = P^(y2)  {Q3(X2) 


(•-)^ 


' -ni/2  „ . m , , -inX,/2 . ”^2  ^2  .“*2, 

a a+2+T  + 48>  - ® ^ + 


o\  /2  « _ 

2X«e  ^ li'X  2 ^ - /V* 


. . 2.2  3.3 

IllA  r “^0  ™ '^9  ^ 

I n + — -+ — -+ — -) 

2 8 48  ^ 


m 


e 2 £^(_) 

- (l-|  + 4-f^)e“^^E^(-f  (X+D)  + (l+f+^+|^)e"”'/^EX(-|(X-l)) 


]4 

"4 

2 

-J  J • 

m 

- (l-|)e'"^^EX(-f  (X+D)  + (l+|)e"'"^^E^(-|(X-l))  } 


Since 


X„  X„+l 
Qj(X2)  = Y 


expression  (30.1)  may  be  evaluated  as 


(30.2)  ?^(M)  {|  [L(^)  +E.t(-|(X-l))]  |'(D3)  - i (Dl)j  - |^(D3) 

1*1  4.  A fii  n ^ ft 

M +2  4 <1-2+T'48>  - 7T<1-2  T 

-J  m 3m  J 


- (D3X)  - ~ (Dl)  + (DiX) 
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EU-  f(X+l))-  + 4mX)e""'^''^  - x[{ 

m ■- 


(D2X)  - 


“inX/Z  "1 

^-— ]} 

3(f)  J 


where  the  notation  of  Table  4 has  been  used. 

The  corresponding  B term  for  k = 1 in  Eq.  (22)  is 


(31)  I 


31 


= Pj^|>(a)]Q^[X(b)]e“"'^l/2 


(31.1)  13^  = P^(M)  {(Dl)|'£n(|^)  + Ec(-  |(X-1))  - 


— (1--) 
2 U 2 
m 


e”^^  Eli-  f(X+l))  + 


-mX/2  -1 

^ 

(t)  -* 


The  next  term  in  Eq.  (22.5)  is 


a" 


(32)  I 


32 


P2D(a)]Q2[X(b)]P2(U2)e''"^l^^  (X^  - 3)  dX^ 


In  the  notation  of  Table  4, 


(32.1)  I32  = P2(U2) 


'>2»2> 


r3[(D4)-(Da)]  - 2[(D2)-(D2X)] 


3m 


e-m/2  _ g-mX/2l  "| 


[e  - , j j 


+ P2^^2^ 


-mX/2  , 2,2  3,3 

/-  , inA  , in  A , m A v 

(1+—  +-  Q—  + --^,0  ) 


(-)^ 

''2'^ 


2 8 48 


I + I { (D4X)  In  (^) 


'x-r 


2 3 A 

- 24(^)(l-f  + ~-fg  + ^)e'"^^  Eli-  |(X+1))+  (D4)  Eli-  |(X-1)) 
m 


3 2 

-,,32,,m  ,2 , , 3ni  -mX/2-.  1 r/p.,,,  „ ,a+i^ 

- 24 (--5)  (^  X +-^  + -^)e  } - 2 {(D2X)  £n(x::3;) 


X+1, 


2 

+ Ei(-  f(A+l))  + <D2)  s.i(-  f(A-l))  - 

m m 


+ -g"  ~ f(^+l))  + Ei(-  f(A-l))}|} 


>] 


Now 

Q^a)  = 


so  terms  can  be  combined  in  the  above  expression  to  obtain  a simpler 
form  for  the  integral. 


(32.2)  PjCKj)  ((•^^)  in(p|)  r3|(D«)-(Da)l  - 21  (D2)-(D2X)  1 


+ A [e-"/ 


-m^/2  o -mX/2  2 2,2  3,3 

eA  9e  /..m.mX.mX 

2 ~r  ^^■*'T-*‘~ir+"4r 


) 


(y)  (f) 

+ |(D4A)  E-c(|— ) - 18(^)(l-|+^^-|^+-|^)e“''^  E^(- f(A+l)) 


+ |(D4)  El(-  f(X-l))  - 18(^)(S^  + ^ + .^)g-n.X/2 

m 


- |(D2A)  -en(|^)  + -^(l-|i+^)e”‘^^  &c(-  |(A+1)) 

m 

- -1(02)  E4(-  f(A-l))+ 

m 


ra/2  -m/2 

® EZ  (-  f ( A+1) ) + (-  f ( A-1) ) > } 


6m 


-( 

2 


Equation  (32.2)  may  be  simplified  as 


(32.2)  P2(b) 


6^  [^-<Fr)  !<'-«)] 
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2 3 A 

®_4...S[ 

2"  8 ■ 48  38A 


m , m in  , m ..  8,,  m m -..ll  m/  2 m, , . s 

r)(i--+— -7^+-^)  ^ ® E4(-2(a+i)) 


- [(isx^) 

I-  m 

e"”^''2  rfl-HnX/2)  9(l-hnX/2-hnV/8+m^X^/48)  144  ,m V , m^X  . .3m  , , ^ \ 

+ 2"  2 2 "T^"96‘+  8 +4  ^ 


(-) 


(-) 


The  B,  term  for  k = 2 is 
k 


(33)  I33  ■ P3(M3) 


(33.1)  - PjIUj) 


+ P2{U2) 


r-  . 2 -1 

/OO 

3X,-1 

1 

Xo 

4 

P,[>(a)lQ,ti(Wle  dl, 


3X^-1 

-Q2(X2)(-T-)e“”'^^^^  dXi  + P2(X2) 


£n(|^)  e dX^) 


r“  3X  I/O 
(-  dX 

X2 


(33.2)  1^3  = P2(y)t  (^^3~'^) 


3 (D2) -3 (D2X) - (D0)+(D0X)+3 (D2X) - (DOX) 


2 2 
+ [3(D2)-(D0)j  EZ(-f(X-l))  + 


48^,  m , m 2 
“(I-I-^T)  " m 
m 


m-  f(X+l))  - 


^(f)  + 6(4)  (1  + 

m m 


-mX/2-. 
e ; 


If  one  regroups  and  cancels  terms,  Eq.  (33.2)  becomes 

2 

(33.3)  I = P,(P)  (-^)(  13(D2)-(D0)){  fa(j^)  + Ei(-f(»-l))) 


- (l?+4-^)  Ei(-  ?(X+1))  - 4 (2+^)e-'’^'2) 


3 2 8 24 


m 


Proceeding  to  the  next  level  of  integrals  in  Eq.  (22.5),  k = 3, 


the  A term  is 
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(3A)  I 


34 


(34.1)  l3^  - PjCm^)  {-^ 


3^  - 

2 ■ 2 9 


X»  1- 


- /oe> 

7X?  X,1 

1 

- X 

_4 

6 ^ 4_ 

-m\j^/2 


5xJ  - Y + \ 


0 /Xjd.,  -raX-i/Z 
£n(-j^)  e ■‘-  dXj^ 


dX, 


■"'^1/2  dX, 


Again,  Tables  4 and  B-1  are  used  to  evaluate  the  integrals  of 
Eq.  (34.1).  There  are  at  this  point  nine  separate  and  distinct 
terms  in  tthe  evaluation  of  the  integral  of  Eq.  (22.5).  These  are 
collected  in  Eq.  (22.7)  and  are  further  simplified  in  Eq.  (22.8). 


r 

(22.7)  ^Kl  ° ^3 ^^2^  ^ 2 ' 


5(D5)-5(D5X)--y  (D3)+~  (D3X)+(D1)-(D1X) 


+ - 4 |5(D5X)  £n(gi)  + 5(D5)  Ei{- 

t/  2 3 4 5 

r/iort\/04\/^  in,n)  m , id  m \ ni/^  t-»/  i^nn 

- 5(120)(-g)(l- 2 + y-:^+384-  38-^>® 

m 

2 

- 5(120)  (^)  (2 (|)V  + 16(|)V  + (y+72)(|)^X  + (3m^+72)  (|))e""'^2/2l 

m J 

- P3(X2)  [|(D4X)  - |(D2X)  + |(|)e"‘"^2/2"| 

P„(X  ) r , X„+l  1 / nc  2 3 

- -3-(D3X)  ■^^(Y-ij;)  + -y(D3)  E.t(-  -(X^-D)  - 2’’'T“  48  ^ 

4 *-  2 m 

2. 

m X, 


m/2  ra,.  ..  14/96.. 2 , , . -mX,/2 

2 E-c(-  - —(—)(— 2"+  4m)e  2' 

m 


biim 
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r „ A r 

+ — ^—  (DlX)  ^nCy-zj)  + (Dl)  Eii-  fCX^-D)  - -^a-\)e 
- 2 m 


X2+I 


4 zi.m/2 


EZ(-f(X2+l))  } 


P,(X) 

(22.8)  = P3(P2){ 


{5(D5)  - y(D3)  + (Dl)}{£n(|^) 


+ E4(-f(X-l))} 


/ 1/0^  23  ,,  2345 

f ^ /I  ^ N 14/^0^  /,  in  j tu  m \ /.|  id  , n id  , in  m ^ 

- l-^(l-2)  “T^— ^°°^"6^^^'2’^T“48'*-m‘  3840  ^ 

mm  m 

e"*^^  Eli-  f(X+l))  - {700(^)(2^  + 2mV  + (S|-  + 72)  ^ 

m 


j-  /"i  2 , X^/'96v/mX  . / \ . 1 1 -mX/2 

+ (3m  + 72)  (-j))  - 3‘(“4)(~2"  + ® 

m 

- {240(^)(l+f  + a^  + ffi^+!^)  + 

m ra  -* 


- (• 


5X' 


-f)  [= 


5(D5)  - y(D3)+(D1)-5(D5X)  + y(D3X)-(D1X)J  } 


14, 


"1. 


The  B term  for  k = 3 in  Eq.  (22.5)  becomes 


f<a 


(35)  I 


35 


P3(M2)P3[^(a)  ]Q3[X(b)  dX^^ 


(35.1)  I35  = P3(h2) 


'2  dx^  + PjCXj)!"  w 
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(35.2)  = P3(M2)  ['{5(D3)-3(Dl)}{£n  (^)  + E^-^CX^-l))) 

2 3 

f A80  m m m . 12  /,  m.-,  m/2  c-,  m..  .t... 

- i~  “ ~~2  ® 

m m 

-^df+  _ io(„2i)  - i|  r"^2/2| 

m -' 

5X2 

- (-^  - -f)  |'5(D3)-5(D3X)-3(Dl)+3(DlX)j} 


At  this  point,  the  desired  level  of  approximation  for  the 

integral  has  been  reached.  It  is  now  necessary  to  take  the  terms 

developed  in  Eq.  (22)— (35)  and  include  these,  with  the  proper 

coefficients,  in  Eq.  (22.6),  for  I , and  carry  out  the  integration 

KX 

over  the  second  set  of  coordinates.  Thus 


(17.5.5.1)  = I I^i(X,(2)xi<2)  - x;(2)Xb(2))  dVj 


is  the  integral  that  must  be  solved. 

At  the  present  level  of  approximation,  the  p coordinate  appears 

til 

in  the  factor  only  through  P^(P),  the  i*"  Legendre  polynomial. 

One  property  of  these  Legendre  polynomials  is  that  each  contains  only 
even  or  odd  powers  of  the  argument,  as  i is  even  or  odd.  This  property 
allows  some  simplifications  for  the  integral. 

Note  that 


(36)  XJ2)X:(2)  - X:(2)XJ2) 


- /aV- 


-ar  -Yr,  -yr  -ar, 
a b a b 

e e - e e 


>] 


(35,1)  = ['g-(a+Y)RX/2|j  |■^-(a-Y)Rp/2  _ ^(a-Y)Rp/2j 
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The  second  wavefunction  is  an  odd  function.  Since  the  integral 
over  M is  carried  out  between  symmetric  limits, 


(37) 


r= 

^-1  L 


-(a-Y)RlJ/2  _ e(a-Y)RM/2 


P^(y)dM  = 0 


and 


(38)  I I e 
-1 


r[ 


-(a-Y)RM/2  (a-Y)RV!/ 


- e 


P^(Vi)m' 


dp  = 0 


for  even  values  of  i.  Thus  the  even  terms  in  the  expansion  for  I 

K1 

contribute  nothing  to  the  final  integral. 

The  lowest  term  remaining  which  contributes  to  the  integral, 

after  combining  terms  (see  Eqs.  (30.2)  and  (31.1)  for  individual  term 
values)  is 


„ r 48A. 

(39)  = 2/aV  R P^(P){' 


m"  . m^.  ,m . . m^  , m^. 


(1-  — + (!-—)(-) +—  + —) i (1+ 

4 72  72  2 8 48  2 

“ m m 


Xe 


-ml  2 


- ra 


^ (1-|)1  Xe”/2  e^(-2.(A+1))  + ^ k(4-!|)  - B (^)1 


-mX/2 


l-96Ai 


2 3 


_i  ^ 

L m m 


e-m/2} 


Equation  (39)  may  be  simplified  by  letting 


(40)  T 


11 


I-  48Ai  „2  J.  ™ n,2  3 2Bi 

-ru- 

■-  m ra  - 


8 3 


(41)  T 


12 


-^a---a--^)ij)+  8 48^  2 ^ 2i| 

L m m 


Equation  (39)  then  becomes 


(39.1)  135  = 2/a^Y^  R^P^(y)  ■*■ 


- 2T  e“'"/2} 
11 


To  complete  the  calculation  of  the  first  term  in  integration 

must  be  performed  over  the  remaining  electron  coordinates: 


-1  1 


Consider  the  y integral  first.  Pj^(y)  “ P>  so 


I-  ^2  2 r-(a-Y)RM/2  (a-Y)Ry/2' 

(43)  J (X  -y  ) [e  - e 


ydy  = 2(X  Bg^j  - Cg^j) 


where 
(44)  B 


gu 


■ (f)  r !-[» 

-1 


-(a-Y)Ry/2  (a-Y)Ry/2‘ 


- e 


2" 


dy 


(44.1)  B 


gu 


(cx-y)^R^/4 


and 
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(45)  C =4) 
gu  '2 


(45.1)  C = + 


-1 


-(a-Y)RM/2  ^(a-Y)RM/2' 


- e 


2“ 


dVi 


(a-Y)R/2 


gu  (a-Y)V/16 


6 - 3(a-Y)R  +7-(a-Y;  V 
4 


3 3 

(g-Y)  R 
8 


-(a-Y)R/2 


3„3 


— — — [s  + 3(a-Y)R  + -(a-Y)^R^  + 

(a-Y)VVl6  L 


8 


With  these  definitions,  the  integral  over  X can  be  obtained: 


(43.1)  Ix^B  -C  1 

L gu  guj 


e-“^/2  {X 

^ ir 


X ^ (-  f(X-l))] 


- T e 
12  ® 


,m/2 


Eli-  f(X+l))  + ^ 


-mX/2 


- 2T  dX 

11 


(43.2)  aW  {BgJ  [£h(|^)  + E4(-|(X-1))J  e""’^/^dX 


- B 


gu 


I T e"''^^X^  Ec(-  f(X+l))e"™^^^dX  - 2B  [ ' 

12  2 gu 


Tiie‘"’/2e-'"^/2^2^X 


+ B 


f (M)e-"'X/2 


gu  Jl 


m 


- 2 9 — 

A //  ® ^ Tj  A \2  “inX/2  .A 


- c 


gu 


Zni^)  +BL(-  f(X-l))] 


-inX/2  . 
e dX 


+ CgyJ  I’xie^^^Xe  ^ Eli-  ^(X+l))dX  .(.2c 


gu 


T ^“"*X/2 

e dX 


- C 


guJ 


(^)  e dX  } 
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The  majority  of  the  integrals  of  Eq.  (43.4)  appear  in  Tables 
B-1  and  4.  The  remaining  integrals  will  now  be  evaluated  in  steps 
(44)-(52). 


(44) 


'00 

j^3e-mX/2 

1 ^ 


To  integrate  by  parts,  let 


dv  = X^e“"^^^dX 


„ . . 1|  ,-"X/2  ^ ^ (|)2  ^ ^ 


U = El(-  |(X-1)) 

-m/2  (X-1) 

du  = 

(X-1) 


Integral  I may  then  be  replaced  by 


(44.1)  I = [uv]”  - f”vdu 

E 1 n 


Noting  that  (-  ^(X-1))  cannot  be  evaluated  at  the  lower  limit 


and  replacing 

this 

by 

tim 

Et(- 

?(B-1)).  Eq. 

(44.1)  be 

6-^1 

L 

(44.2) 

I = 
E 

tim 

6-»-l 

m 

g-m/2 

^48 

2 

mm. 

X)  El  (■ 

- f(6-l)) 

^00 

-mX/2 

fm^X 

3 

<—  X 

2^2 

IB,  L.  4. 

mX  , 

-|  -m/2  (X-1) 

1 ® 

dX} 

^41® 

L^8 

8 ^ 

2 + 

1 • 

J (X-1) 
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Consider  now  the  remaining  integral  of  Eq.  (44.2).  This 
integral  can  be  evaulated  by  making  a change  of  variables,  X-  1 = z. 


»5)  I45  = 


+ ^ (z-1)  + 1|  dz 


m 


1 


(00  -mz  I-  3 
e ^ 

5-1  2 m 


2 

(z^-3z^+3z-l)  + ^ (z^-2z+l) 
o 


(45.1)  I, 


The  Insertion  of  the  result  of  Eq.  (45.1)  into  Eq.  (44.2)  yields 


(44.3, 

m 


2 3 

o * /96  _*"m/2\  /-i  X ® .1.  ^ X 

tun  {-r  e ^ T T 15^ 
6-^1  m^  2 8 48 


[ec  (-  f (6- 


D)  -El  (-ni(6 


-1))] 


Step  (44.3)  has  to  be  substituted  back  into  Eq.  (43.2).  In 
Eq.  (43.2)  the  integral  over£n(-|^)  must  be  evaluated  before  pro- 
ceeding to  the  limit.  The  value  for  the-t*^  integral  is 


(46)  {e““/^(l  + I + ^ + |i)  tun  \ In  {2)  -In  (5-1)  + EH-  5(6-1) )1 

m S-*"!  ^ -* 


,,  m . m^  m^.,  -m/2 

- (1  - 2 + T " 48^  ® Et  ( 


-m)  - i[4  + (f  “ <?)] 


.-m/2 


Thus  the  first  term  in  the  K integral  (Eq. 

gu5 

including  coefficients,  is 


(17.3.5.1), 
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2 0 

(47)  B T (^)  e"™  f ~ 4 52  + 1 _ e“™  (1  - 2 + SL  _ SL) 

gy  1C  1 <i~o  /.o' 


16  16  3 


2 8 48' 


E^(-m) 


+ e'^d  + 2 + m + m ) T £jj(2)  - £n(6-l)  + 2 EX(-  -(6-1)) 

641  ^ ^ 


- EK 


-m(6-l))]} 


The  last  term  in  Eq,  (47)  has  already  been  evaluated  in  Appendix 
B,  and  has  the  limiting  value 


C + £n(^)  - Zn(2)  + £n(2) 


giving  a final  value  for  Eq.  (47)  of 


(47.1)  B T (^) 
gu  11  ^4 


r -m  .m^ 

L®  % 


. 5m  . 2v 

+ IE  + 3^ 


- 2 + T - !S) 


+ e-'^d  + 5 + ^ + |^)(C  +^n(f)) 


The  corresponding  terra  in  Eq.  (43)  with  coefficient  C may  be 

g^ 

solved  in  an  identical  manner,  with  a final  result  of 


r -m 

(48)  C T (—)'  ^ “ 


4 \ I _ /I  _ m'^  m^,  , m^  „ ^m 


gu  11  n,2  L 2 


- (1  - |)  EX(-m)  + e-'d  + j)(C  +ln  (f))] 


The  integrals  of  the  Et(-  ^(X+1))  terms  of  Eq.  (43.2)  will  be 
considered  next.  The  approach,  as  above,  is  by  an  integration  by  parts 


/OO 

(49)  El{-  f(X+l))dA 
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dv  = 


V = (1  + "X  + m^X,^  + !n^)(ii)p-in^/2 
''■^  2 8 48  ' 


u = Eli-  |(X+1)) 


-m/2(X+l) 

du  = ^ dX 

(X+1) 


This  gives 


(49 . 1)  I . (2|)  + !^  + 1 + 1)  Ei(-m) 

m 


+ [”(96.  -mX/2  5sV  + slA!  + + i)  i 

j,  4^®  '■48  8 2 

1 m 


^ -m/2  (X+1) 
(X-1) 


The  remaining  term  in  Eq.  (49.1)  is  integrated  by  making  the 
substitution  X + 1 = z. 


(50)  ^ f 

m"* 


9|  g+m/2  e U (^3  - 3z^  + 3z  - 1)  + ^^  (^2  _ ^ 

J 2 z IJf8  o 


+ ^ (2  - 1)  + 1 dz 


= 26  iV”"  p(^)  + z(^)(l  - |)  + |(1  - S + 


+ (1  - f + 1 


2 3 


-53)<7>_ 
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96  m/2 
- e 

m 


3 Oo"2m 

(=^)  V 

. A8 


2 -2m 

(l+2m4-2m2)  + (^){1  - H)  §— (1  + 2m) 


2'  2 
m 


2 -2m 

fB')  n _ B + B_^  Q 

(2HI  Y + 8 ^ m 


m/2 


96e‘"'"  m m^  cv  Om'i 

“ " "2  ■*■  T " E^(-2m) 


m 


Combination  of  like  terms  and  multiplication  of  this  integral 

by  its  coefficient  in  the  K term  yields  the  final  result 

gu5 


2 3 

+m/2  f96  1"  -m/2  ,,  ,m,m  .mv^./^v 

"12  V 7 b 5 T ^ 

Ei(-2»)1  + T (t  I + <> 

— ‘ m 


m , m^  m^. 

- " 2 + T ■ ® 


e-3tn/2} 


Identical  techniques  give  the  solution  of  the  equivalent  C 
Integral  of  Eq.  (A3. 4) 


gu 


(52) 


(52.1) 


'^12^gu  J 


= T.nC 

12  gu 


^m/2^e“'"^/^  £c(-  •|(X+l))dX 

-ra 

{— 2"  (1  + Et(-ni)  + ^^2 
m m 


/ ™ 

^ (1  - ^)  E4(-2m)} 


The  remaining  integrals  of  Eq.  (43.2)  are  tabulated  in  Tables 
B-1  and  4. 

Substitution  of  these  terms  into  Eq.  (43.2)  gives  the  first 

non -zero  contribution  to  the  K . integral.  As  the  integrals  were 

gu5 

to  be  used  in  a computer  program  for  finding  the  numerical  value  of 
the  energy  as  a function  of  internuclear  separation,  no  additional 
simplification  was  attempted  with  this  integral. 


90 

The  next  term  which  yields  a non-zero  contribution  is  the  k = 3 


term  in  of  Eq.  (17.3.5.1).  The  integral  multiplying  the 
coefficient  in  Eq.  (22.6)  is 


0 r i,, 

(53)  I53  = r2p3(p)  |{5(D3)  - 3(D1)}  itni^)  + (- f (A-1)) } 


_ (i_^+  !L  _ ^m/2  £^(_^x+l)) 

m4  2 8 48  2 2 


- {Mlil  ^m_X  ^ g-mX/2  _ 3^o(d2X)  - 

m4  2 3m  _ 


5A. 

2 


-f] 


2 2 2 2 3 3 

5(D3)  - 3(D1)  - {(1  - -SL-)  + (1  - 7!r:)(’V)  + 

4 120  120  2 8 48 


m 


Proceed  by  again  considering  the  p integral  in  Eq.  (17.5.5.1) 
ueing  Eq.  (53)  and  (32.6).  Now 


(54)  I j (X^  - [;-(a-Y)I>M/2  . ^<a-v)I!p/^ 


dX  = X B + C , 
gu3  gu3 


where  B „ and  C „ are  determined  in  Eqs.  (55)  and  (56),  respectively. 
gu3  guj 

(55)  . (-1)  - 3P) 


(55.1)  B =5  I 

gu3  IS 


.3„3 


2_2 


(a-Y)''R''  3(a-Y)"R”  I 


+ 3(a-Y)R  + 6 


3 


J (a"Y)^R^/16 


_ 5 + S (q-Y)-^-  + 3(a-Y)R  + 6 

8 4 


-(a-Y)R/2 

e 

(a-Y)^R^/16 
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- ■) 


[<2:^-1] 


(a-Y)U/2 

e 

21 — 

R'‘/4 


+ 3 


TT 


-)  e 


-(a-Y)R/2 


The  integrals  required  for  evaluating  B are  listed  in  Table  B-1. 

gu3 

Similarly,  C ^ is 
gu3 


(56)  = i [3u='  - 5m^ 

(56.1)  C , = . 3(2^  3.  3 . ,-]  ^(c-y)!>/2 

(a-Y)^R^  L 8 4 J 

- toiV  + n ,-(c.-y)R/2, 

. _W0  . 5<sa^  + 5(oa)V  . I5(a-Y)'R' 

(a_Y)6R6  L 32  16  2 


+ 60(a-Y)R  - 120]  ^ — 

J L 32  16  2 

22  , ^ I -(a-Y>R/2, 

+ 15(a-Y)  R + 60(a-Y)R  - 120  e } 


With  these  definitions,  the  integral  of  Eq.  (54)  can  now  be 
evaluated  by  integrating  over  the  final  coordinate.  The  integral 
of  Eq.  (54)  is 

(54.1)  ^ f [x^B  + C 1 [(5(03)  - 3(D1)} 

2 l_  gu  gu3J  4 |_ 

{ ^n(M)  + ei(-  1))}  - {2110^(1  _ S + ^ _ ]:|(1  - S)} 

A-l  ^4  2 8 Z 

X fiii-  2(X+1))  - {3QU8)  ^ g-my/2  _ io(D2X)  - — 

^ ^ ^ .jin 
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,5X2  2.  r 


- - f)  [5(03)  - 3(D1)  - ((1  - 2^)  + (1  - gj)(f) 


. m^X^  . m^X^  , -raX/2  "1 , 


For  convenience,  let 
(57)  = 5(D3)  - 3(D1) 


(58) 


T - r / I H 4.  ssi  N jL2  , ra.  ~|  m/2 

The  integral  of  Eq.  (54.1),  on  substituting  for  P3(X),  becomes 


(54.2)  (x\„3  + C^„3)  {(^)’-i3  + eU-  f(X-l))] 


- (^V^)Ti4  EK-  ^(X+l))  + (^  - 1){ 


5(D3X)-3(D1X) 


-[i 


- 3« 


. (5|?  . 2,,^^  . (5^,  20,5^  ^ ,„,^-nX/2, 


m 


e-"'^/2dx 

Each  term  in  the  integral  of  Eq.  (54.2)  will  be  evaluated 
separately.  Taking  these  in  order, 


^59  = "nij 


[£>i(|^-|)  +Eli-  ^(X-1))]  e "^^^dx 


5T  2 3 

(59.1)  I = Zm—~  {(D3)[£n(2)  -.£ia(6-l)]  - f 

59  6-)-i  ^ m 


•*■  (23)  Eli-  ®(6-l))  - “4  + ^m)e-"'/2} 


m 


8 
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+ I EU-  ?(X-l))dX  - tim-P  {(Dl)t£n(2)  -Zn(S-l)] 

8 Ji  2 8 

- 4 (1-?)  E^(-m)  + (Dl)  Ei(-  ?(6-l)) 

m2  ^ 2 

3T 

- X eu-  ?(x-i))dx} 

8 2 


The  integrals  of  Eq.  (59.1)  can  be  integrated  by  parts. 
In  the  first  integral,  let 


. >3  -mX/2 

dv  = A e 


dX 


96  -mX/2.,  4.  mX  , m^X^  , m^X^s 
2 8 48  ^ 


u = E4(-  f(X-l)) 
-m/2(X-l) 


du 


(x~i) 


With  this,  the  E-4(-  |(X-1))X  integral  in  Eq.  (59.1)  becomes 


2 3 


(60)  I = ton  {2^  + 

59  6->l  2 8 48  2 


+ 26 
m^ 


« „3i3  2,2  , , -m/2(X-l) 

1 8 2 


In  the  last  integral,  let  z = X - 1. 
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/ /■“  r 3 2 2 "I 

/r -1  \ 96  I I m 2 j.  in  /it.  ni\  i ni/-i.i.  tn  t.  ^ \ I 

m ' 1 - 


e""'^dz 


+ i^  e 


-m/2 


i:  t-" 


™ + IIL  + 
8 


3-  , 

m /_1\ 

48 J 


The  first  of  these  integrals  in  Eq.  (61)  is  equal  to 


(62)  2|  e-m/2  ji  + i (1+|)  + i(l  + I + I|-)1 


and  the  second  is  equal  to 


(63)  - ^ e““^^  (1  + ? + V + ^)  E^(-m(6-l)) 

2 o Mo 


By  similar  methods,  the  other  integral  over  E^C" 
Eq.  (59.1)  is  evaluated. 


(64)  tan 
6-»-l 


A p-n>/2 
m 


2e 


-in/2 


^ e-^'^  (1+B)  E^(-  f(5-l))  + „ 

2 2 2 2 


4e' 


-m/2 


(1+f)  E4(-  >n  (6-1))] 


Noting  that  the  coefficient  of  those  terms  which  are  ill- 
behaved  at  the  lower  limit  in  Eqs.  (59.1),  (63)  and  (64)  is  the 
same  for  each  pair  of  integrals,  the  discontinuity  can  be  reduced 
to  evaluating 


tun  -'En(6-1)  + 2Zn(-  ^(6-1))  - £n(-m(6-l)) 
6->-l  L ^ 


(65) 
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This  is  the  same  term  which  is  evaluated  in  Appendix  B;  its 
value  is 


(65.1)  -in(2)  +£k(|)  + C 


Using  this  result  in  each  of  the  two  pair  of  ill-behaved  terms 
of  Eqs.  (60),  (63),  and  (6A), 


Pi  2 

3m  ^ 

+ 4 _ _ 5^}  - 2 

_ 8 

8 2_  ® 

Wl3 


, m/2  -m/2 

4^  (l-f)Ei(-m)  +2^  ) 


Taking  the  integral  over  X in  Eq.  (54.2)  next,  the  same 
sequence  of  steps  is  followed. 


(66)  I 


66 


dX 


(66.1)  I 


66 


"^X^  £m(~)  dX  + 

1 


X^  E£(-  |(X-1))  dX 


The  first  integral  on  the  right-hand  side  of  Eq.  (66.1)  may  be 
evaluated,  using  Table  B-1,  It  is  found  to  have  the  value 


(67)  tur!{(D5)[£n(2)  - £n(6~i)]  + (D5)  E£(-  S(6-l))  - 120^ 


6->l 


m 


r“ 


2 n,3  4 

n “ 4-  SL  EL  4.  J3L_ 

2 8 48  384 
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5 

ji_ 

3840 


4 2 2 

+ CV  + 2m^+  + 72)  (5^)  + (3m^  + 72)(^)) 

o 2 4 ^ 


The  other  integral  in  Eq.  (66)  can  then  be  integrated  by  parts, 


with 


dv  = X5 


dX 


V = 


-I20e 


-raX/2 


, 2.2  3.3  4.4  5.5 

n 4-  4-  ^ I ^ ^ 4-  ^ 4.  5LJl_') 

m.6  ^ 2 8 48  384  3840^ 


(-) 

'■2'' 


u = EU-  |(X-1)) 


du  = e 


-m/2(X-l)  dX 
(X-1) 


It  is  found  that 


(68)  I 


68 


x5e-mX/2  |(X-1))  dX 


-m/2 


(l  + ? + 'v  + ^ + 


6-»-l 

120(64) 


m 


2 8 48  384  3840 


) Eli-  f(6-l))} 


m 


\ In  \ 2.2  3.3  4.4  3.5  “’^(X“l)/2 

-mA/2  . mX  , ra  X , m X , ra  X , m X ^ e ^ 

- U + o + Q + + 3g^  + 3g40‘'  - 


(A-l) 


dX 


The  terms  in  Eq.  (68)  are  then  evaluated,  yielding 


"HI  f 2 I"  0 

(68.1)  [jfe  + i (1+f)  + id+l  + 'l') 


- Ml  + s + ^ ti(-  f(s-i)) 

6*^1 


in'- 


Combination  of  Eq.  (68.2)  and  Eq.  (60.1)  gives  the  value  of 


Ijg  in  Eq.  (59.1)  as 


^eul'^is  [(DS)  (£n(f)+c)  - “ 2 ■*■  T " fs 


gu3  13 
5 


m 


384  3840 


m 


-m/2  4 _ 3 

r ' 4 2 


) E^(-m)  I-  (-  T “ ^ - 36m  + 3 + ^(1+^) 


2 2 3 4 - 

+ |(2-hn+I^)  + |(6  + 3m  + ^ + !|-)  + |(24  + 12m  + 3m^  + ^ + ^)) 


4 8 

3 4 


2 3 


- I ['(D3)(Mf)+C)  - (1  . m ^ !L.  _ |_)  U{.-m) 


I6e 


-m/2 


m 


^ (4  + i|£5  + ^ - 4m)|} 


The  next  of  the  major  terms  in  Eq.  (54.2)  are  presented  in 
Eqs.  (69) — (74.2).  This  term  is 


(69)  Ui-  ?(>.+l))  dA 


gu3  ^gu3^ 


The  method  of  solution  will  be  similar  to  that  used  on  the  ill- 


behaved  integrals  just  considered,  i.e.,  expansion  by  parts  follow -id 
by  a substitution  of  variables. 
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>00 

(69.1)  1,0  = (5X3-3X)  E.c(-  5(X+1))  dX 

oy  j 1 ^ 


= 5 


.00  /CO 

X3  m(x+i))  dX  - 3 X (_  |5(X+i))  dX 

1 a “ 


-m/2 


(69.2)  Igg  = 5 [(DS)  Ec(-m)  + ^ ^ 2 “ f 


2 „3  nr 

+ ^ - 7^)E^(-2n.)  - 3 (Dl)  EX(-m)  + 

8 48  L 


2e 


-3m/ 2 0^/2 


2e 


m 


(1+^)  E^(-2m)'| 


Similarly, 


(70)  5 r X^  EU-  f(X+l))  dX  - 3 j^X^ 

' 1 

(70.1)  = 5 [(D5)  EX(-m)  + 


-3m/2  4 o o 2 , o 

64e  /SL  + „3  . 2m_  + 3m.^  3 

6 ''2  2 2 4 


+ |(1-|)  (m3  + ^ + ^ + |)  + 5(2  - m + 5^)(m3  + m + |) 


2 4^ 


+ 5(i+2m)(6  - 3m  + - 5L)  + |(24  - 12m  + 3m3  - !|-  + ^)) 


2 3 A 5 ” 

64(120)  /I  m L m m jn _ m „m/2  p/c_9m'\ 

- -6  " 2 ■*■  T “ 45  + 384  3840^  ® 


m 


- 3 


IP  -3m/2  2 qa„"’/2 

(D3)  EZ(-m)  + (^  + i + ^)  - 


m ~2 


m 


8 8 


n _ •_  + !!i_ 
U 2 8 


- a(-2m)] 


The  next  term  in  Eq.  (54.2)  is 
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(71)  I (X2Bg^3  + Cg^,3)(^  - i)  [5(D3X)  - 3(DlX)| 


-mX/2 


dX 


(71-1)  = 5 B (S2) 

m 


, - 2,2  3,3  , 

X^6  + 3mX  + ^ e""'^  dX 

4 o 


2^gu3 

(i^) 

m2^  . 

2,2 

3, 

1 3m  X 

+ slA 

4 

8 

X^l  + 5^)  dX  + (|Cg^^3 


- i\u3>  C ^ 


dX 


3 '^gu3  ' 


,00 

m 1 


(6  + 3«X  f 32^  + S^)  e-"'  dX 

4 o 


3 ^gu3 
® m 


(1  + y)  e dX 


These  integrals  are  tabulated  in  Table  4.  Using  the  terms 
defined  there  and  rearranging  the  result  algebraically  yields  the  form 


(71.2)  {25  ||(C7)  + i|(C6)  + ^(C5)  + ^(C4) 

^ L m"^  m - 


- 15 


20 

3 


+ -(C5)l 
m J 

_ 2Q 
3 

[-(05) 

Ira 

+ y(C4)] 
in  " 

r^(c3) 

+ 2|(C2)] 

+ 4 

Wi.01) 

-tn 

Tn  - 

-fn 

+ (25 


-(C5)  + i§(C4)  + ^(C3)  + 

m^  ra"*  m 


5(C2)] 


15 


■y(C2)  + ^(C3) 


4n 


>1  - f f- 

J 3 Iji 


(C3)  + y(C2)  + ^(Cl)  + ^(C0)1 
m'^  m 


+ 4 


-^(CO)  + -(Cl)|} 


>1 
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The  next  term  t.'  be  considered  in  Eq.  (54.2)  also  has  its 
component  integrals  given  in  Table  4: 


™ l“fV3 


. ,^^gu3"^\u3.  r ^2  r40w,  . mX 

+ ^ T " 5^1 

2,2 


+ '-2-5^-) 


m_l!x  _ Al  e-'"^  dX 
3m 


-i^Su3 


e"*"^  dX 


If  Eq.  (72)  is  evaluated  by  use  of  the  tabulated  integrals 
and  combining  terms,  its  value  is 

(72.1)  B , r^(C7)  + ^(C6)  + (1^  - |5)(C5)  - ^(C4) 

8u3  L2m  n,2  ^3  6m  ^2 

- - 2)(C3)]  + C^„3  ^(C5)  + ^(C4)  + - g)(C3) 

- ^(C2)  - - ^,)(C1)] 

m^  m'^ 

The  next  term  from  Eq.  (54.2)  is 


^3'"l3 


' (5f  - H Cg„3)  e-"*/?  dX 


= (^)  r ^ ("x^  dX 

J-j  Z Z -5  j 


^ T,3  f e-^/2 
3 ) I 


dX 
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Equation  (73)  is  evaluated  as 


(73.1) 


73 


5B  5C 

— ^ Tj^3(D4)  + T]_3 


2B 


Ru3 


2C 


)(D2)  - 


^ Tno(DO) 


^13' 


The  final  term  in  Eq.  (54. 2)  is 


/■"“  ...  3 


(74) 


-mX/2 


dX 


Referring  to  Table  4 for  the  solutions  to  these  integrals, 
and  rearranging,  the  final  result  becomes 


(74.1)  B 

L2 


(C6)  - i|(C4)  + i^(C5)  - i^(C3)| 

m"^  -• 


m 


m 


+ C 


gu3 


Z5(c4)  - 15 (C2)  + 1^(C3)  - i^(Cl)] 
m^  m^  m^  m"^ 


Substitution  of  Eqs.  (59.2),  (67),  (68.1),  (69.3),  (70.1), 

(71.2),  (72.1),  (73.1),  and  (74.1)  into  (54.2)  yields  the  next  terra 

in  the  series  for  Eq.  (17.3.5.1).  As  the  calculation  of  the  energy 

was  to  be  done  numerically  on  a digital  computer,  each  term  was  used 

in  the  form  presented  above  without  explicit  back-substitution. 

From  comparison  of  the  final  value  of  the  K term  with  values 

gu 

at  T = 1 given  by  Huzinaga  (49) , and  for  large  R,  it  was  determined 
that  this  was  a sufficient  number  of  terms  to  give  an  overall 
accuracy  which  was  adequate  for  the  calculation  of  the  ground  state. 


varying  at  worst  by  only  about  five  percent  from  the  best  numerical 
solutions  available.  The  work  involved  in  adding  additional  terms 
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to  improve  the  accuracy  of  this  integral  (the  only  one  where  any 
approximation  has  been  necessary)  did  not  seem  commensurate  with 
the  improvement  in  accuracy  attainable  with  the  additional  terms. 

For  the  smallest  values  of  R considered,  the  total  contribution 
of  this  term  was  less  than  one-tenth  of  the  total  energy,  dropping 
to  less  than  five  percent  in  the  range  of  internuclear  separations 
of  interest,  so  the  error  involved  in  this  truncation  was  small. 

The  final  determination  of  the  ground  state  energy  curve  was 
accomplished  by  using  a computer  program  to  evaluate  the  terms 
given  above  in  Eq.  (11)  for  values  of  R,  the  internuclear  separa- 
tion, from  0.2  to  15.0  Bohr.  The  screening  constants  a and  y 
were  chosen  at  each  value  of  R to  produce  a minimum  in  the  energy 
curve.  The  basic  scheme  used  by  the  program  was  to  first  select  an 
R value,  and  then  to  compute,  term-by-term,  the  value  of  the  various 
interactions  for  values  of  a ranging  from  /2  to  2,  the  range 
appropriate  to  the  screening  constant  for  hydrogen.  From  a considera- 
tion of  the  physical  meaning  of  the  screening  constant,  it  is 
evident  that  the  values  of  y appropriate  to  the  problem  are  between 
0 and  /I,  as  the  value  a = y = v2  gives  a minimum  in  the  energy 
curve  for  large  R.  This  corresponds  physically  to  the  equivalence 
of  the  two  electrons  in  a hydride  ion  (H~) . The  program  then 
recalculates  the  energy  for  another  pair  of  a and  y values,  finally 
storing  the  energy  minimum  and  the  values  of  a and  y corresponding 
to  this  minimum.  Since  this  portion  of  the  calculation  was  to  be 
achieved  by  use  of  the  computer,  no  attempt  was  made  to  further  com- 
bine terms  beyond  the  level  presented  earlier.  In  hindsight,  it  is 
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evldeiu  Llwit  there  would  have  been  some  virtue  in  having  combined 
terms  further:  for  small  values  of  R,  round-off  errors  began  to  be 
a considerable  factor  in  the  calculations  on  the  computer.  This 
was  particularly  evident  in  the  evaluation  of  the  Bi(-x)  terms, 
where  double  precision  calculations  were  necessary,  in  the  evalua- 
tion of  £n(S)  and  C,  and  in  the  highest  of  the  terms  in  the 
integral,  where  the  problem  of  round  off  was  never  completely 
eliminated.  This  problem  was  an  additional  inducement  to  truncate 
the  series  at  the  point  chosen,  as  it  was  evident  that  the  problem 
would  continue  to  be  a troublesome  one  for  the  even  more  involved 
higher  power  terms.  Had  it  been  possible,  extended  precision  (21 
or  more  decimals)  would  have  avoided  the  worst  of  the  problems 
encountered  with  round  off,  but  the  computing  facilities  currently 
available  do  not  support  this  arithmetic. 

One  additional  word  concerning  the  programs  Is  in  order.  Some 
effort  was  made  to  reduce  the  time  required  per  run  of  the  program, 
which  would  have  taken  some  two  to  four  hours  of  computer  time  for 
evaluating  the  ground  state  energy  for  R values  spaced  0.05  Bohr 
apart.  This  was  achieved  partially  by  first  determining  the  rough 
curve  shape,  and  then  by  selecting  R values  which  were  widely  spaced 
in  regions  where  the  energy  curve  was  smooth,  using  closely  spaced 
values  only  in  regions  of  large  curvature.  This  allowed  a great 
reduction  in  machine  usage,  but  necessitated  a somewhat  more  Involved 
program  for  plotting  the  results. 
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C.  RESULTS  OF  MOLECULAR  ION  CALCUUTIONS 

The  ground  state  energy  of  the  doubly  negative  hydrogen  molecular 
Ion  can  be  calculated  using  the  functions  developed  In  Section  IV  B. 

The  details  of  the  results  are  given  here.  Five  cases  will  be  considered 
separately;  these  Illustrate  effectively  the  results  obtained  and  Indi- 
cate the  levels  of  accuracy  In  the  calculations.  Of  these  cases,  two 
are  chosen  prlmarally  to  determine  the  accuracy  of  the  numerical  results 
by  comparison  with  previously  known  values,  and  a third  Is  chosen  for 
the  Insight  which  It  sheds  on  the  results  obtained  for  the  doubly 
negative  Ion  ground  state.  These  will  be  considered  In  the  following 
order:  neutral  hydrogen  ground  state  with  fixed  and  then  with  variable 
screening,  the  doubly  negative  Ion  with  and  without  variations  In  the 
screening,  and  then  the  He^  ground  state,  for  comparison  with  the 

results  for  the  K Interactions, 

gu 

Case  I.  Neutral  Hydrogen  with  Constant  Screening 

The  ground  state  energy  of  the  neutral  hydrogen  molecule  was 
determined  using  the  previously  developed  terms  for  comparison  with 
known  results,  to  verify  the  accuracy  of  these  terms.  The  calculation 
Is  derivable  from  that  of  the  doubly  negative  Ion  by  the  presence  of 
only  two  electrons  in  the  molecular  orbitals.  Tne  value  obtained  for 
this  energy  as  a function  of  internuclear  separation  is  given  as  Fig.  10. 
Note  that  the  energy  approached  asymptotically  is  not  the  anticipated 
value  of  -2  Rydbergs,  but  is  (-2.0  + ■!■  ) Rydbergs.  To  obtain  the  value 
of  -2  Rydbergs  requires  that  the  molecule  separate  into  two  neutral 
atoms.  However,  in  the  MO  approach,  the  molecule  may  also  separate  into 
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+ H . Such  a system  differs  in  internal  energy  from  the 
system  by  the  calculated  (^)  Rydberg.  This  point  is  discussed  in  Slater (51) 
in  great  detail,  and  the  result  of  the  present  calculations  are  in 
agreement  with  the  values  given  by  Slater  for  this  approach.  The 
screening  constant  a has  the  value  /z  . 

For  finite  internuclear  separations,  the  values  found  by  this 
calculation  are  close  to  those  obtained  by  the  best  available  wavefunc- 
tion  choice  (52) , and  for  small  separation,  will  be  about  as  good  as  can 
be  obtained  for  a simple  wavefunction  choice,  wherein  the  screening 
constant  is  held  fixed. 

Case  II.  Neutral  Hydrogen  with  Variable  Screening 

In  order  to  obtain  more  accurate  energy  values,  the  screening 
parameter  may  be  varied  with  separation  by  choosing  values  for  the 
parameter  to  minimize  the  calculated  energy  at  each  separation  R.  As 
seen  in  Fig.  11,  the  behaviour  at  infinity  is  improved,  although  it 
still  shows  the  behaviour  characteristic  of  MO  calculations.  Figure  12 
shows  the  variation  of  the  screening  parameter  with  internuclear 
separation  for  the  calculations  of  Fig.  11.  This  variation  of  the 
screening  parameter  is  characteristic  of  this  system  (51);  as  the  elec- 
tron "orbit"  will  be  that  predicted  by  the  Bohr  model  for  infinite 
internuclear  separation,  and  will  gradually  lessen  its  distance  from  the 
center  of  the  nuclei  as  the  two  nuclei  approach  each  other  closely. 

This  leads  to  an  increase  in  the  screening  parameter  with  decreasing 
distance  (as  the  screening  parameter  is  inversely  proportional  to  the 
separation  of  the  electron  from  the  nucleus) . The  energy  minimum 
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obtained  by  this  approach,  -2.27  Rydbergs,  is  consistent  with  the  better 
values  given  by  MO  calculations  cited  by  Slater  (51) . Values  obtained 
by  Wang  (53)  by  a similar  method  give  an  energy  minimum  of  -2.278  Ryd- 
bergs, In  each  case,  the  position  of  the  energy  minimum  is  at  about 
1.4  Bohr  intemuclear  separation. 

The  above  comparisons  with  previous  work  on  the  neutral  hydrogen 

molecule  were  performed  to  test  the  validity  of  the  present  computation 

_2 

scheme.  The  transition  from  the  four-electron  calculation  of  the 
molecular  ion  to  the  two-electron  molecule  involves  contraction  of 
the  state  function,  but  individual  terms  remain  unaltered. 

-2 

Case  III.  Ground  State  Energy  with  Fixed  Screening 

If  the  electrons  in  the  two  types  of  orbitals  are  equivalent 

(i.e.,  have  the  same  screening  constant),  the  curve  of  the  energy  as  a 

function  of  intemuclear  separation,  shown  in  Fig.  13,  does  not  possess 

a minimum  for  finite  values  of  the  intemuclear  separation.  The 

screening  constants  used  for  calculating  this  curve  (a=Y=*^)  were  chosen 

to  give  the  lowest  energy  for  infinite  separation  of  the  two  nuclei. 

At  small  separations,  the  energy  is  everywhere  above  that  of  the  neutral 

_2 

hydrogen  molecule.  Indeed,  the  energy  of  the  H2  is  greater  than  that 
of  the  separate  H ions  removed  to  infinity.  Thus  the  curve  of  Fig.  13 
actually  lies  in  the  range  of  free  states  of  the  neutral  molecule  (see 

_2 

Fig.  11).  As  the  separation  tends  toward  infinity,  the  energy  of  the 
drops  below  that  of  the  neutral  molecule;  this  is  in  accord  with  the 
known  stability  of  the  Hydride  (H  ) ion.  The  energy  obtained  for  very 
large  R values  agrees  acceptably  with  the  known  binding  energy  of  the 


Kig.  L3.  (Iround  8 


molecular  ion  with 


lii 

-2 

excess  electron  of  the  Hydridve  ion.  The  asymptotic  energy  of  the 
system  should,  on  the  basis  of  the  single  hydride  ion  calculation  of 
Chandrasekhar  (55),  be  about  2.05  Rydbergs.  The  present  value  of  1.85 
Rydbergs  is  in  accord  with  the  asymptotic  results  using  similar  tech- 
niques and  wavef unctions  obtained  by  Fischer-Hjalmars  (48)  for  the 
ground  state.  The  values  reported  there  for  the  equivalent  method 
(her  Method  II)  give  the  results  of  1.85  Rydbergs  for  the  asymptotic 
value  of  the  energy. 

Case  IV.  Doubly  Negative  Hydrogen  Molecular  lOn  with  Variable 
Screening 

The  next  set  of  curves.  Figs.  14,  15,  and  16,  illustrate  the 

results  obtained  for  the  hydrogen  molecular  ion  when  the  screening 

parameters  are  allowed  to  vary  with  R.  Figure  14  gives  the  ground  state 

energy  for  the  molecular  ion  as  a function  of  internuclear  separation. 

Figure  15  gives  the  screening  parameters  corresponding  to  the  energy 

curve  of  Fig.  14.  Of  more  importance,  however,  is  the  behaviour  of  the 

curve  in  the  region  about  7 Bohr.  Figure  16,  which  is  a composite  of 

Figs.  11,  13,  and  14,  shows  that  the  energy  curve  of  Fig.  16  is  actually 

a composite  of  the  two  curves  of  Figs.  11  and  13,  with  a transition 

occurring  at  the  point  where  the  two  curves  cross. 

The  behaviour  of  the  screening  constants  alpha  and  gamma  also 

gives  insight  into  what  is  happening.  It  will  be  noted  that  the  shape 

of  the  variation  of  alpha  with  separation  in  Fig.  15  is  much  the  same 

as  for  the  parameter  if  Fig.  12,  with  alpha  increasing  smoothly  with 

-2 

decreasing  separation.  In  Case  III  (H^  ),  the  two  constants  were  held 
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constant.  In  Case  I (H2) , only  one  screening  parameter  was  required 

_2 

and  it  was  made  a function  of  R.  In  case  IV  (H^  ) , two  screening 

parameters  were  used,  each  a function  of  R.  The  second  screening  cons- 
tant, gamma,  in  Case  IV  is  the  same  as  alpha  for  R greater  than  7 Bohr, 

but  gamma  drops  abruptly  to  zero,  (or,  rather,  to  its  minimum  allowed 
value)  for  smaller  values  of  the  internuclear  separation.  This  implies 
that  the  second  electron  is  no  longer  bound  to  the  system,  as  the 
screening  constant  may  be  interpreted  as  (1/r^),  where  r^  is  the  separ- 
ation between  the  electron  and  the  nucleus  to  which  it  is  attracted  at 
which  the  radial  portion  of  the  molecular  orbital  wavefunction  has  fallen 
to  (1/e)  of  its  maximum  value.  Thus  the  minimum  of  Fig.  14  is  actually 
the  minimum  for  the  neutral  hydrogen,  and  is  not  the  value  which  the 
doubly  negative  molecular  ion  would  have  at  an  equivalent  separation 
(that  is  best  approximated  by  the  curve  of  Fig.  13).  Thus  what  the 
energy  function  is  actually  doing  is  properly  portrayed  by  Fig.  16,  where 
a transition  has  occurred  between  the  doubly  negative  ion  and  the  neutral 
hydrogen  at  an  internuclear  separation  of  about  7 Bohr.  This  transition 
occurs  by  the  ejection  of  two  electrons  from  the  system,  resulting  in 
a new  system  consisting  of  the  neutral  molecule  plus  two  free  electrons. 
The  very  considerable  energy  difference  between  the  two  curves  in  the 
region  of  the  minimum  of  the  neutral  curve  indicates  the  lack  of  sta- 
bility of  the  molecular  ion.  Although  the  approximations  inherent  in 
the  choice  of  the  wavefunction  do  not  allow  the  exact  position  of  this 
transition  to  be  determined,  its  existence  and  approximate  location 
must  be  close  to  the  indicated  7 Bohr.  The  magnitude  of  the  difference 
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between  the  two  curves  at  the  position  of  the  neutral  hydrogen  minimum 
is  believed  to  be  correct. 

The  above  may  be  interpreted  in  terms  of  hydride  ions  which 
approach  each  other  from  infinity.  When  their  internuclear  repulsion 
becomes  sufficient  to  remove  the  extra  electrons,  allowing  the  formation 
of  the  neutral  molecule  from  the  molecular  ion.  In  constructing  Fig.  16, 
the  curve  of  the  energy  using  variable  screening  has  been  shifted 
slightly  with  respect  to  the  axes  to  make  it  stand  out  more  clearly  from 
the  other  two  courve,  so  that  it  does  not  directly  overlay  either  of 
these  lines.  It  is  readily  seen  that  this  curve  first  follows  the  curve 
of  Fig.  13  and  then  crosses  over,  following  the  curve  of  Fig.  11  at 
about  7 Bohr.  There  appears  to  be  a small  transition  region  near  the 
crossover  in  which  the  energy  determined  for  the  curve  allowing  vari- 
ation in  the  screening  parameters  is  somewhat  below  either  of  the  other 
two  curves. 

This  calculation,  on  the  basis  of  the  initial  wavef unction 
choice,  does  not  treat  the  additional  electrons  on  the  hydride 
separately.  It  is  therefore  impossible  in  this  calculational  scheme 
for  only  a single  electron  to  be  ejected  from  the  system.  However,  the 
values  obtained  by  Fischer-Hjalmars  (^iS)  for  the  ground  state  energy  of 
the  system  would  seem  to  indicate  that  the  singly  negative  ion 
ground  state  energy  is  near  that  of  the  neutral  hydrogen,  but  probably 
somewhat  greater  (i.e.,  less  strongly  bound).  It  is  thus  probable  that 
the  results  indicated  above  would  have  occurred  had  the  calculation 
allowed  the  possibility  of  the  ejection  of  electrons  singly.  However, 
in  that  case  there  might  have  been  an  additional  transition  region. 
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Case  V.  Neutral  He^  Molecular  Calculations 

As  an  additional  check  on  the  accuracy  of  the  terms  used  in  this 
calculation,  a calculation  was  also  made  of  the  ground  state  energy  of 
the  He^  molecule  (or,  if  one  prefers,  the  interaction  between  two  He 
atoms) . The  calculation  is  identical  in  principle  to  the  one  presented 
by  Huzinaga  (49) , so  that  the  values  of  the  terms  obtained  here  may  be 
checked  directly  with  those  given  in  his  paper.  A comparison  of  the 
values  obtained  for  this  calculation  and  by  Huzinaga  are  given  in 
Appendix  D.  The  total  energy  for  the  He  system  is  shown  in  Fig.  18, 
while  the  electronic  terras  only  are  shown  in  Fig.  17.  The  values  of  the 
screening  constants  alpha  and  gamma  are  shown  as  a function  of  internu- 
clear  separation  in  Fig.  19.  As  the  values  used  by  Huzinaga  are  tabu- 
lated values,  his  results  are  limited  in  the  scope  of  variations  of 
the  screening  constants,  so  that  the  results  obtained  here  are  somewhat 
more  accurate  for  these  constants.  The  difference  of  about  2 % between 
the  values  obtained  by  Huzinaga  and  those  obtained  by  this  calculation 
for  the  energy  at  an  internuclear  separation  of  1 Bohr  are  almost 
entirely  due  to  differences  in  the  value  for  the  term  wiiere  a 

series  of  approximations  are  necessary  in  order  to  obtain  a "closed 
form"  function  of  the  screening  constants  suitable  for  performing  the 
variation  of  parameters.  To  obtain  a more  accurate  energy  curve  for 
this  system,  it  would  be  possible  to  correct  this  term  at  selected  points 
by  a numerical  integration  using  the  values  of  the  screening  constant 
obtained  from  the  approach  used  here.  The  additional  accuracy  would 
be  significant  only  for  small  internuclear  separations,  as  this  term 
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drops  off  rapidly  with  increasing  separation.  Figures  17,  18,  and  19 
serve  to  demonstrate  the  essential  accuracy  of  the  individual  teirms 
used  in  the  evaluation  of  the  molecular  ion  ground  state  energy,  and 
I to  give  limits  on  the  accuracy  of  the  calculation. 
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Y.  DISCUSSION 


If  space  charge  in  an  insulator  plays  an  Important  role  in  dielec- 
tric breakdown  by  causing  the  growth  of  a highly  dense,  highly  energetic 
gas  in  thi  solid,  there  remains  the  necessity  of  explaining  the  manner 
in  which  the  energy  of  the  space  charge  is  concentrated  into  a suffic- 
iently small  volume  of  the  dielectric  to  allow  the  start  of  channel 
formation.  Two  points  must  be  considered:  firstly,  the  space  charge 
distribution  prior  to  formation  of  a gas,  and,  secondly,  the  space 
charge  distribution  during  the  actual  channel  formation,  after  the 
process  has  been  Initiated.  Since  defects  serve  as  trapping  sites,  the 
generation  and  movement  of  defects  requires  further  analysis. 

A.  Effects  of  Charge  Imbalance  on  Breakdown 
1.  Initial  Charge  Distribution 

There  are  several  lines  of  Investigation  which  suggest  the  impor- 
tance of  an  initial  imbalance  of  the  charge  distribution.  The  effect 
noted  by  Budenstein  and  Hayes  (3,17)  in  investigations  of  breakdown  sites 
by  visual  observation  with  an  electron  microscope  is  perhaps  the  most 
directly  suggestive  of  this  importance.  Wlien  the  capacitor  under  inves- 
tigation was  illuminated  in  the  region  of  a "likely"  site  for  the  occur- 
rence of  breakdown,  the  breakdown  event  would  always  occur  at  this  site. 
If  some  other  region  of  the  capacitor,  not  containing  any  such  site,  was 
Illuminated,  the  breakdown  never  occurred  in  the  illuminated  region.  The 
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beam  of  the  electron  microscope  was  therefore  insufficient  to  cause 
breakdown,  but  the  extra  imbalance  in  the  charge  was  sufficient  to  deter- 
mine the  breakdown  site.  Studies  of  dielectric  materials  with  much  higher 
Intensity  high  energy  electron  beams  (23,24)  show  the  patterns  of 
breakdown  characteristic  of  electrical  breakdown  in  the  region  of  pene- 
tration of  the  beam.  In  the  beam  experiments,  an  additional  effect  is 

noted  in  polycrystalline  samples the  breakdown  will  frequently  occur 

along  or  in  the  vicinity  of  grain  boundaries. 

Also  of  relevance  in  this  connection  are  high  energy  beam  studies 
in  plastic  dielectrics.  In  such  cases,  when  a beam  (several  MeV)  illum- 
inates the  plane  surface  of  a plastic  disk-shaped  specimen,  breakdown 
ultimately  occurs  with  the  breakdown  channel  originating  on  the  illumi- 
nated surface.  The  observed  channel  goes  directly  to  the  space  charge 
layer , then  runs , with  branching , along  this  layer  and  parallel  to  the 
surface  (11) . Localization  of  the  breakdown  channel  in  a plane,  perpen- 
dicular to  the  electric  field,  that  contains  a high  space  charge  density 
is  strongly  indicative  of  the  Importance  of  space  charge  in  the  breakdown 
process. 

The  mechanism  for  space  charge  build-up  under  application  of  an 
external  electric  field  is  frequently  assumed  to  involve  carrier  injec- 
tion in  regions  where  there  is  i geometrical  enhancement  of  the  applied 
field  due  to  asperities.  Thus  for  example  Smith  (10)  assumes  a hemispher- 
ical electrode  projection  causes  a field  strength  of  three  times  that  of 
the  region  near  the  smooth  portion  of  the  electrode.  Ridley  (45),  using 
a more  complex  geometry,  shows  that  Smith’s  results  yield  a lower  limit 
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to  the  enhancement  of  the  field,  a more  pointed  projection  giving  values 
substantially  higher.  If  the  injection  mechanism  is  Fowler-Nordheim 
emission,  then  the  (local)  injected  current  density  is  a very  strong 
function  of  the  electric  field  in  the  vicinity  of  the  injecting  point: 

J - J C exp(-k/E)  ) 

where  only  the  dependence  of  the  current  density  on  the  electric  field  is 
indicated.  Thus  only  a small  increase  in  the  electric  field  can  give 
rise  to  a very  large  change  in  the  injected  electron  density.  If  the 
injected  charges  are  trapped,  then  the  cathode  field  is  weakened. 

However,  if  the  field  continues  to  rise,  a point  will  be  reached  when 
traps  are  filled.  Subsequent  injected  electrons  will  then  be  accelerated 
for  long  distances  and  be  able  to  acquire  the  energy  necessary  for  cre- 
ating ionization  and  defect  centers,  such  as  the  Vj^  center. 

Point  defects  have  mobility  and  will  move  in  response  to  a concen- 
tration gradient  by  diffusion.  Mechanical  stress,  the  presence  of  dislo- 
cations, grain  boundaries,  and  other  point  defects  will  influence  the 
rate  of  diffusion.  If  the  defect  contains  an  unbalanced  charge  or  a 
dipole  moment,  then  it  will  be  influenced  by  the  presence  of  the  electric 
field  or  its  gradient,  respectively. 

Dislocations  and  grain  boundaries  are  regions  of  major  misfit. 

Such  regions  tend  to  trap  defects,  thus  immobilizing  . am.  Cooper  (19) 
has  shown  that  carefully  annealed  single  crystals  of  alkali  halide 
e'.hlbit  about  one  third  of  the  scatter  in  breakdown  field  strengths 
compared  to  similar  unannealed  specimens.  Annealing  does  not  appreciably 
alter  the  centroid  of  the  distribution.  If  unannealed  crystals  tend  to 
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concentrate  charged  defects^  then  a lower  breakdown  strength  Is  antici- 
pated. If,  on  the  other  hand,  the  built-in  imperfections  strongly  Inhibit 
motion  of  defects,  then  the  defect  mobility  is  decreased  and  it  is  more 
difficult  to  reach  the  critical  charge  density  requisite  to  start  the 
growth  of  a gaseous  channel.  This  explains  how  defects  in  unannealed 
crystals  broadened  the  range  of  breakdown  voltages  without  appreciably 
altering  the  mean  breakdown  strength. 

Cooper  (20)  has  also  shown  that  the  breakdown  channels  of  carefully 
annealed  single  crystals  tend  to  follow  certain  crystallographic  direc- 
tions, even  when  the  field  is  45  degrees  from  these  directions.  In 
unannealed  crystals,  the  channels  follow  the  field  directions.  If 
mobile  charged  defect  centers,  such  as  centers,  are  generated  contin- 
uously at  a site  and  if  these  migrate  more  easily  along  the  favored 
crystallographic  directions  (as  is  known  for  H centers) , then  a "spike" 
of  Vj^  centers  is  produced  along  this  direction.  If  the  density  of  charge 
reached  the  critical  value,  then  conversion  of  the  dielectric  to  the  gas 
would  occur  along  the  spike  and  the  resulting  channel  would  be  along  the 
favored  crystallographic  direction. 

2.  Space  Charge  Distribution  After  Channel  Initiation 

Following  the  initiation  of  a gas  pocket,  the:  charge  distribution 
will  alter  considerably  from  that  prior  to  the  formation  of  the  first 
gas  pocket.  If  the  assumption  is  made  that  the  channel  will  form  in  the 
region  of  highest  initial  charge  density  imbalance,  then  one  expects  that 
th<*  region  containing  the  pointed  electrode  projection  will  be  the  point 
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of  initiation  of  the  channel,  Previously,  the  dominant  process  for 
removing  the  defect  centers  from  the  region  has  been  diffusion,  Once  a 
gas  pocket  has  forroed,  the  effect  of  the  Internal  field  inhomogenletles 
within  the  dielectric  will  begin  to  have  a more  marked  effect  on  the 
distribution.  As  was  shown  in  Chapter  III,  after  a gas  pocket  has 
formed,  there  will  be  forces  on  the  centers  tending  to  concentrate  them 
in  the  region  of  the  gas  pocket.  This  will  lead  to  a lessening  of  the 
net  diffusion  from  the  region,  and  may  lead  to  a reversal  of  diffusion, 
as  the  defects  are  now  incorporated  within  the  gas  pocket.  Further,  the 
gas  pocket  will  offer  a lower  resistance  to  the  passage  of  current, 
thereby  leading  to  further  enhancement  of  injection  from  the  electrodes 
in  this  region.  The  greater  density  of  current  in  this  region,  and  the 
immediately  surrounding  solid  dielectric,  will  lead  to  an  enhancement 
in  the  production  of  centers.  The  additional  centers  will  be  placed  in 
the  dielectric  in  a manner  leading  to  their  rapid  incorporation  within 
the  gas  pocket,  leading  to  its  growth.  In  addition,  the  greater  current 
injection  entering  the  pocket  will  traverse  the  gas,  and  will  perforce 
have  to  leave  the  pocket  through  largely  undisturbed  dielectric  material. 
This  will  lead  to  the  creation  of  further  centers  within  the  solid  along 
a line  determined  largely  by  the  field  orientation  and  by  the  orientation 
of  the  gas  pocket  relative  to  the  material.  The  combination  of  these 
two  effects,  then,  serve  to  establish  the  defect  density  in  the  region 
v;here  it  is  most  likely  to  lead  to  further  growth  of  the  gas  pocket. 

If  the  gas  pocket  leads  to  a higher  rate  of  creation  of  centers 
in  this  fashion  by  acting  as  a "new"  cathode  (as  is  suggested  by 
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Lopez  C56)  ) , then  one  might  expect  there  to  be  another  such  pocket  . 
formed  further  into  the  dielectric,  spaced,  relative  to  the  first  gas 
pocket,  much  as  the  first  gas  pocket  is  spaced  from  the  cathode.  In  such 
a situation,  both  pockets  would  then  repeat  the  process,  etc.  Throughout 
this  time,  each  gas  pocket  would  be  growing  as  it  is  formed,  increasing 
its  voliane  both  through  the  direct  process  of  incorporating  defect 
centers  directly  into  the  gas,  and  by,  it  may  be  assumed,  a process  of 
heating  at  the  walls  making  use  of  the  energy  liberated  on  the  incorpo- 
ration of  the  defect  centers.  Such  a process  would  be  likely  to  give 
rise  to  a channel  having  a "beaded"  appearance,  such  as  is  actually  obs- 
erved in  some  breakdowns. 

The  presence  of  extensive  defects,  such  as  slip  planes,  where  def- 
ects may  have  been  concentrated  during  the  initial  charge  build-up 
stages  of  the  process,  will  also  have  a considerable  effect  on  the  form 
and  direction  of  the  channel  within  the  dielectric.  Such  planes,  in 
single  crystals,  are  normally  oriented  along  certain  directions  of  the 
dielectric.  This  may  lead  to  a preferred  direction  of  the  breakdown 
channel  in  such  an  orientation  as  to  incorporate  the  maximum  available 
charge  within  the  channel.  The  exact  direction,  however,  has  not  as  yet 
been  determined  to  give  rise  to  this  condition.  The  existence  of  such 
extended  defect  trapping  sites  in  polycrystalline  samples  would  thus 
depend  on  their  orientation  with  respect  to  the  field. 
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B.  Relevance  of  Molecular  Results  to  Breakdown 

The  gaseous  breakdown  model  requires  a solid  to  become  energetically 
unstable  at  sufficiently  high  charge  densities.  External  excitation  must 
produce  intermediate  states  within  the  solid  which  can  decay  into  a mole- 
cular  gas.  If  the  H2,  H2  » and  H2''"^  molecules  are  taken  as  prototypes 
for  molecular  bonds  with  different  numbers  of  extra  electrons,  some  gen- 
eral conculsions  can  be  drawn  from  the  calculations  on  the  ground  state 
energies  of  the  hydrogen  negative  ion  sequence. 

The  bonding  energy  in  the  ground  state  of  the  neutral  hydrogen 
moleucle  at  its  equilibrium  internuclear  separation  of  1.4  Bohr  is 
0.347  Rydbergs  (51),  or  about  4.7  eV.  Fisclier-Hjalmars  (48)  has 
calculated  the  ground  state  energy  for  the  H2"^  ion,  obtaining  a dis- 
sociation energy  of  2.54  eV  at  its  equilibrium  separation  of  1.65  Bohr. 
The  minimum  of  the  ^2”^  ion  curve  lies  above  that  of  the  H2  by  1.9  eV, 
with  the  remainder  of  the  difference  in  bonding  energies  being  due  to 
the  stability  of  the  hydride  (H  ) ion.  The  H2“’^,  then,  represents  a 
possible  intermediate  state,  analogous  to  the  CI2”  state  discussed  in 
Chapter  III. 

_2 

The  results  for  H2  show,  for  the  case  in  which  the  screening 
parameters  are  held  fixed,  an  energy  at  the  position  of  the  neutral 
minimum  of  about  -0.75  Rydbergs,  or  about  21  eV  above  the  neutral  mini- 
mum. The  state  is  highly  repulsive,  and  could  not  form  during  breakdown 
except  as  a short-lived  resonant  state,  as  it  lies  far  into  the  continuum 
range  of  states  for  neutral  hydrogen.  The  repulsive  character  of  the 
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state  is  sufficiently  pronounced  that  the  interaction  of  an  additional 
electron  with  the  112"^  ion  may  cause  detachment  of  the  excess  electron 
in  the  negative  ion,  yielding  a neutral  molecule  plus  a free  electron, 

If  the  additional  electron  Interacts  elastically,  then  the  freed  electron 
would  have  an  energy  up  to  1.9  eV,  The  action  of  the  calculated  ground 
state  energy  of  H2  with  a variation  of  screening  parameters  also 
supports  this  conclusion.  The  calculations  indicate  the  doubly  negative 
ion  changes  into  an  molecule  plus  free  electrons  at  an  internuclear 
separation  of  7 Bohr.  Thus  the  H2''^  state  becomes  energetically  unfav- 
orable at  large  internuclear  separations.  The  value  of  the  internuclear 
separation  for  the  conversion  gives  an  approximate  scale  of  interaction 


ranges  for  the  interaction  of  a free  electron  with  the  singly  negative 
hydrogen  ion,  as  the  extra  electron  would  have  to  be  accomodated  within 
the  bonding  system  of  the  H2~^  ion.  This  implies  that  the  interaction 
would  create  an  intermediate  state  during  the  interaction  of  112“^,  and 
so  gives  a range  for  the  interaction.  It  is  not  possible  from  the  cal- 
culations presented  here  to  obtain  significant  energy  information  from 
the  state  energy  of  the  H2”  ion  obtained  by  the  variation  of  parameters, 
due  to  the  ionization  of  the  extra  electrons  from  the  ion  and  consequent 
conversion  into  neutral  hydrogen,  as  discussed  in  Chapter  IV. 

From  the  progression  of  ground  state  energies,  it  can  be  seen  that 
a sufficient  number  of  excess  electrons  introduced  into  the  bonding 
system  will  lead  to  the  disruption  of  molecular  bonding.  This  is  shown 
in  Table  5,  where  the  ground  state  energies  are  given  relative  to  a 
zero  energy  for  separated  nuclei  and  electrons  at  infinity.  The  intro- 
duction of  two  electrons  into  the  bonding  system  in  addition  to  those 
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TABLE  5 

GROUND  STATE  ENERGIES  AND  EQUILIBRIUM  INTERNUCLEAR 
SEPARATIONS  FOR  THE  HYDROGEN  NEUTRAL  AND  NEGATIVE  ION 
MOLECULAR  ION  SEQUENCE 


SPECIES 

GROUND  STATE  ENERGY 
(Rydbergs) 

EQUILIBRIUM  SEPARATION 
(Bohr) 

-2.73 

1.40 

-2.21 

1.65 

V' 

-0.75 

none 
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normally  present  in  the  neutral  is  sufficient  to  disrupt  totally  the 
bonding.  In  the  case  of  excited  states,  Taylor  et  al  (57)  have  shovm 
that  the  singly  negative  ion  possesses  no  bound  excited  states;  this  is 
also  obviously  true  with  regards  to  the  doubly  negative  system. 

The  energies  associated  with  bonds  in  condensed  matter  vary  from 
a few  tenths  of  electron  volts  to  several  electron  volts.  If  the 
binding  energies  of  the  112°,  systems  are  taken  as 

representative  of  other  molecular  charge-variable  systems,  a perspective 
is  provided  for  the  proposed  gaseous  breakdown  mechanism.  The  addition 
of  an  extra  electron  to  a single  bonding  pair  decreases  the  bond  stabil- 
ity by  about  2 eV  with  only  a small  relaxation  of  the  Internuclear 
separation.  An  additional  electron  makes  the  system  strongly  unstable. 
Instability  arises,  then,  for  an  excess  charge  of  between  one  and  two 
extra  electrons.  The  case  considered  in  the  detailed  analysis  was  that 
of  an  isolated  molecule.  In  condensed  matter,  molecular  bonds  between 
two  adjacent  atoms  are  not  completely  Independent  of  those  of  surround- 
ing atoms.  If  excess  electrons  are  added  in  low  concentrations  in  a 
solid  Insulator,  they  generally  become  trapped.  In  so  doing,  they  modify 
the  local  energy  level  structure.  Each  such  electron  will  contribute 
to  the  electrical  potential  energy  of  a-1  atoms  containing  trapped 
electrons.  If  the  trap  density  becomes  high  enough,  then  the  perturbing 
potential  would  reach  the  condition  where  nonbonding  orbitals  are  pro- 
duced. (In  the  sense  used  here,  the  term  nonbonding  includes  both 
nonbonding  and  antibonding  orbitals.)  This  condition  is  exemplified  by 
the  lack  of  stability  of  the  molecular  ion. 
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An  estimate  pf  the  charge  density  necessary  for  the  production  of 
nonbpnding  orbitals  may  be  obtained  from  breakdovm  experiments  on 
electron-djeam  irradiated  plastic  specimens,  If  a thick  slab  is  irradi-^ 
ated  by  a beam  having  an  energy  of  several  MeV,  the  breakdovm  occurs 
when  the  field  due  to  the  trapped  charges  is  10^-10^  V/cm,  and  the 
discharge  path  goes  from  the  upper  surface  dovm  to  the  charge  layer  and 
through  this  layer  (parallel  to  the  surface  and  almost  perpendicular  to 
the  field  direction)  to  the  lateral  surfaces,  A field  of  10^-10^  V/cm 
Implies  a charge  layer  of  lO"''  or  lO"*"  electrons/cm^.  If  this  charge 
is  concentrated  in  a layer  1 micron  thick,  the  average  distance  between 
defect  centers  is  about  200  A. 

Tang  ^ (58)  have  studied  the  increase  in  electrical  conductivity 
with  applied  field  in  silver  azide  (AgN^)  and  explained  their  results 
in  terms  of  formation  of  an  intermediate  N^”*  molecular  state.  Their 
experimental  arrangement  involves  the  study  of  a single  crystal  of 
silver  azide  with  applied  fields  of  up  to  3.5  KV/cm  under  high  vacuum 
conditions.  Electrode  contact~related  phenomena  were  also  studied  In 
their  experiments,  and  mass-spectrographic  analysis  was  carried  out  on 
the  reaction  products  and  evolved  gas  during  the  reactions.  Decompo- 
sition of  the  crystal  was  detected  with  applied  fields  as  low  as 
0.3  KV/cm,  although  the  decomposition  rate  varied  widely  between  samples. 
For  sinusoidally  time  varying  fields,  decomposition  was  not  observed  to 
occur  for  frequencies  above  15  Hertz  and  fields  up  to  1.7  KV/cm  (peak 
to  peak) . When  only  one  contact  was  made  between  the  electrodes  and 
the  crystal,  decomposition  occurred  only  when  that  electrode  was  the 
cathode.  The  decomposition  product  monitored  was  molecular  nitrogen; 
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In  some  samples,  "bursts"  of  nitrogen  were  detected  superimposed  on  a 
general  background  level.  The  decomposition  reaction  proceeds  through 
the  formation  of  an  Intermediate  N^",  having  an  energy  of  about  2 eV 
above  the  ground  state  of  the  system.  In  all  the  azides,  the  breakdown 
will  proceed  slowly  at  lower  field  strengths  than  those  necessary  for 
channel  formation  and  other  features  normally  associated  with  breakdown. 

Capacitors  subjected  to  high  field  strengths  over  long  periods  of 
time  show  aging  effects.  Including  reduced  breakdown  field  strengths. 
Aging  may  be  due  to  Intermediate  states  formed  over  a period  of  time  at 
field  strengths  less  than  those  required  for  breakdown.  Deterioration 
would  be  associated  with  formation  of  a gas  molecule  rather  than  a 
return  to  the  original  stable  lattice.  This  effect  would  thus  lead  to 
regions  in  which  the  dielectric  is  less  stable  than  the  original  mater- 
ial; such  centers  might  serve  also  as  nucleation  centers  for  additional 
intermediate  states  once  they  are  formed.  This  would  thus  link  the  phe- 
nomenon of  aging  in  dielectrics  directly  to  the  breakdown,  with  the 
primary  difference  being  largely  one  of  time  scale. 

C . Summary 

Dielectric  breakdown,  according  to  the  model  proposed  by  Budensteln, 
involves  the  growth  of  a localized  conducting  gaseous  channel  through 
the  dielectric  under  the  action  of  an  external  electric  field.  In  this 
study,  two  questions  relevant  to  the  early  stages  of  the  solld-to-gas 
conversion  have  been  addressed  and  placed  in  improved  perspective; 
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Cl)  H(?w  Is  energy  added  tg  the  solid  so  that  the  solid-'to-gas  conver- 
sion can  occur?  C2)  What  is  the  mechanism  for  the  solid-to-gas 
conversion? 

The  first  question  has  been  explained  by  showing  that  point  defects 
can  be  Introduced  that  have  excess  energies.  The  energy  of  formation  of 
these  defects  comes  from  the  external  electric  field;  electrons  in  the 
conduction  band  of  the  insulator  are  assumed  to  be  accelerated  to 
sufficiently  high  energies  (several  eV)  that  they  can  produce  collisional 
ioziatlon.  and,  thereby,  point  defects.  The  number  of  defect  centers 
required  for  creation  of  a gaseous  breakdown  channel  is  about  the  same 
as  the  number  of  atoms  of  the  solid  that  are  converted  to  gas  in  this 
channel. 

The  second  question,  on  the  nature  of  the  mechanism  of  solid-to-gas 
conversion,  has  been  answered  by  hypothesizing  the  creation  of 
nonbonding  orbitals  when  the  local  excess  charge  density  is  increased 
sufficiently.  To  obtain  quantitative  information  on  the  role  of  excess 
charge  density  on  molecular  bonding,  the  bonding  characteristics  of  the 
molecular  sequence  H2°,  H2~  , and  was  studied.  This  required 

ab  initio  calculation  of  the  molecular  ion.  Both  H2°  and  H2  ^ have 

stable  ground  state  configurations,  while  the  H2  is  strongly  repulsive 
for  all  internuclear  separations.  The  results  of  the  molecular  calcu- 
lations were  applied  qualitatively  to  bonding  within  a solid.  By  taking 
the  field  due  to  space  charge  to  be  of  the  same  order  of  magnitude  as 
the  breakdown  field,  it  was  estimated  that  an  excess  of  one  electron 
exists  for  about  every  10^  atoms  in  the  breakdown  location,  For  the 
rapid  conversion  of  solid  to  gas,  an  energy  of  several  electron  volts 
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per  atom  Is  required.  Defect  centers  have  energies  of  formation  in  this 
range.  Trapped  electrons  in  the  vicinity  of  defect  centers  will  alter 
the  energy  level  structure  at  the  defect.  Nonbonding  orbitals  are 
assumed  to  arise  when  the  trapped  electrons  interact  stpongly  enough 
with  the  defect  center.  This  is  the  sought  after  conversion  of  solid 
to  gas. 

Features  of  breakdown  in  solids  that  have  been  placed  in  Improved 
perspective  by  the  above  considerations  Include  the  following. 

1.  The  formative  period  after  the  Initiation  of  voltage  stress  is 
here  associated  with  the  need  to  store  energy  in  the  solid  so  that 
the  solld-to-gas  transformation  can  occur. 

2.  The  formative  period  is  here  related  to  the  generation  of  point 
defects  due  to  collisonal  ionization. 

3.  Point  defects  should  form  more  readily  in  regions  already 
partially  disordered  and  in  regions  of  high  field.  The  electric 
field  may  be  locally  enhanced  through  geometric  irregularities  and 
space  charge  effects. 

4.  The  evolution  of  nonbonding  orbitals  has  been  followed  for  a 
particular  molecular  system;  H2  with  zero,  one,  and  two  extra 
electrons.  The  latter  case  has  been  evaluated  for  different  inter- 
nuclear  separations,  including  those  of  the  H2°  and  H2~^  molecules. 
The  ideas  of  this  evolution  have  been  applied  to  the  creation  of 
nonbonding  orbitals  in  solids  due  to  trapped  charges  about  defect 
centers.  This  is  taken  as  the  mechanism  for  converting  atoms  of 
the  solid  to  atoms  of  a gas. 
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5,  PQlnt  defects  can  gerye  ag  a source  pf  prebreakdown  recombin- 
ation radiation, 

6,  The  radius  of  a trapped  electron  is  greater  if  the  dielectric 
has  a high  dielectric  constant.  Hence,  a lower  density  of  electrons 
is  needed  to  alter  significantly  the  energy  levels  at  defect  sites. 
Thus,  nonbonding  energy  levels  will  arise  at  a lower  applied  field 
and  the  solld-to-gas  conversion  will  be  initiated  at  a lower  field. 
Dielectric  breakdown  will  occur  at  lower  fields  for  materials  with 
higher  dielectric  constants. 

7.  A material  that  is  disordered  has  a wider  spectrum  of  trapping 
levels  than  a material  that  is  more  ordered.  Deep  traps  will  yield 
more  localization  of  the  trapped  charges  than  shallow  traps.  The 
detailed  trap  distribution  will  play  an  Important  role  on  the 
breakdown  threshold.  This  distribution  will  be  Influenced  by 
a-nealing . 

8.  Defects  can  form  at  lower  than  breakdown  fields,  as  can  space 
charge  due  to  electrons.  Thus  molecular  bonds  can  be  broken  in  a 
limited  region.  This  constitutes  aging  of  the  dielectric.  The 
rate-limiting  process  in  the  aging  is  probably  the  creation  of  the 
defect  centers.  Thus  chemical  effects,  absorption  of  moisture, 
diffusion  of  active  gases  and  similar  perturbations  can  provide 
the  energy  for  the  defect  centers. 
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D.  Suggestions  for  Future  Work 

The  present  work  has  increased  the  credibility  of  the  gaseous  model 
of  electric  breakdown  in  solids,  but  much  work  reitiaina  before  the  model 
can  be  accepted  without  reservation.  Tasks  to  be  done  include  the 
following. 

1.  Quantitative  investigation  of  point  defect  formation  through 
collisional  ionizaiton  and  other  mechanisms.  Analysis  of  stability 
and  mobility  of  defects  in  high  uniform  and  nonuniform  electric 
fields 

2.  Quantitative  investigation  of  the  role  of  space  charge  on  the 
energy  level  structure  of  defect  centers  and  the  creation  of 
nonbonding  states. 

3.  Dynamics  of  growth  of  the  gaseous  channel  prior  to  the 
completion  of  the  channel. 

4.  Growth  dynamics  of  the  gaseous  channel  after  completion  of  the 
gaseous  bridge  between  electrodes. 

5.  Application  of  the  ideas  of  defect  formation  and  bond  breakage 
to  low  fields  and  long  times  to  explain  aging  effects. 
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APPENDIX  A 


i 


EXPANSION  OF  THE  SLATER  DETERMINANT  AND  GROUPING  OF  THE  WAVEFUNCTION 

As  given  in  Eq.  (7),  the  ground  state  wavefunction  for  the 
double  negative  hydrogen  molecular  ion  is  given  by  the  (expanded) 

Slater  determinant 


(l)Qi(l) 

O 

<{.  (2)a(2) 

cS 

<}>g(3)a(3) 

i}>  (4)a(4) 
g 

1 

<l>g(l)B(l) 

(2)8(2) 

O 

<|>  (3)8(3) 

O 

<1.  (4)8(4) 
& 

""  >P}A 

<|.^(l)a(l) 

<!>  (2)a(2) 
u 

<l>y(3)o(3) 

(j)  (4)a(4) 
u 

«1>^(2)8(2) 

<}>y(3)8(3) 

(|.^(4)8(4) 

This  determinant  is-  expanded,  giving  the  next  form  for  the 
wavefunction;  it  will  be  noted  that  here  the  space  and  spin  com- 
ponents of  the  wavefunction  are  now  grouped  separately  rather  than 
kept  together.  The  expanded  wavefunction  is 


(A-2)  tj) 


{<!)  (!)<{•  (2)<t.  (3)(J)  (A)a(l)e(2)a(3)6(4) 
8 8 « “ 


+ 


(l)(|>g(3)(|.y(2)<{i^(4)a(l)e(2)8(3)a(4)  + 4.^(l)<^^(4)(|.^(2)(J.^(3)a(l)a 


(2) 
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a(l)a(2)6(3)B(4)  - <J>„(1)4.„(4)<{>  (2)6  (3)a(l)e(2)a(3)8(4)  + 6„(l)<i>„(2) 

6^(3)6  (4)e(l)a(2)B(3)a(4)  + 6„(l)6„(2)6„(3)6^(A)a(l)a(2)B(3)B(4)  + 

6^(l)6g(2)6y(3)6g(4)B(l)a(2)a(3)B(4)  - 6g(l)6g(2)6y(3)6^^(4) 6(l)a(2) 

a(3)B(4)  - 6y(l)6g(2)6g(3)6^(4)B{l)a(2)B(3)a(4)  - 6y(l)6g(2)6^j(3)6g(4) 

a(l)a(2)B(3)B(4)  + 6g(l)6y(2)6g(3)6^(4)B(l)a(2)a(3)B(4)  + 6y(l)6g(2) 

6^(3)6^j(4)B(l)B(2)a(3)a(4)  + 6^(1)6,  (2)6_(3)6_(4)a(l)B(2)a(3)B(4)  - 

6 (1)6,  (2)6„(3)6^(A)B(1)B^'"  :(3)a(4)  - 6„(l)6_(2)6-(3)6„(4)o(l)B(2) 

a(3)B(4)  - 6y(l)6y(2)6g(3)6g(4)B(l)a(2)a(3)B(4)  + 6g(l)6y(2)6y(3)6g W 

6(l)B(2)a(3)a(4)  + 6,,(l)6g(2)6y(3)6g(4)a(l)B(2)B(3)a(4)  + 't>y(l)<l>y(2) 

6 (3)6„(4)B(l)a(2)B(3)a(4)  - 6„(1)6  ,(2)6  .(3)6„(A)a(l)B(2)B(3)a(4)  - 
g 8 g u u g 

6 (l)6„(2)6„(3)6,,(4)B(l)B(2)a(3)a(4)  - 6„(1)6  (2)6.  (3)6p(4)a(l)B(2) 
118^8  ^^88 


6(3)a(4)} 
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The  space-portions  of  the  wavefunction, 

nth  electron  are  defined  by  Eqs.  (1)  and  (3)  of  Chapter  IV,  and  the 
spin-portion  of  the  wavefunction  is  denoted  by  a(n)  for  spin-up  and 
g(n)  for  spin-down. 

By  regrouping  and  factoring  into  the  (six)  spin  components, 
one  obtains  the  form 

(A-3)  i|<  = — r{4'„(l)<{'„(4)<|.,  (2)4>  (3)  - i (2)6  (4)  + <{.  (1) 

|_gguu  gguu  u 

^ (2)<fr  (3)6  (4)  - 6 (1)6_(2)6  (3)6„(A)1  a(l)a(2)B(3)6(4)  + (1)6  (2) 

6u(3)<tu(^)  - <i>g(l)<{'g(A)6^j(2)6y(3)  + 6y(l)<t>y(2)6g(3)6g(A)  - 6^(l)6g(2) 

6 (4)6  (4)1  a(l)6(2)a(3)6(4)  + r6„(l)6„(3)6„(2)6  (4)  - 6„(1)6„(2) 

guj  I_gguu  gg 

6 (3)6  (4)  + 6 (1)6_(2)6„(3)6_(4)  - 6 (1)6^(2)6„(3)6_(4)]  o(l)B(2)B(3) 

a(4)  + U (1)6^ (2)6,  (3)6^ (4)  - 6p(l)6„(2)6,,(3)6,  (4)  + 6 (1)6,  (2)6 J3) 

6 (4)  - 6 (1)6  (2)6^(3)6^(4)]  B(l)a(2)a(3)6(4)  + [6  (1)6.(2)6„(3)6„(4) 

- 6„(l)<frg(2)6g(3)6„(4)  + 6,,(l)6„(2)6g(3)6g(4)  - 6g(l)6„(2)6^(3)6g(4)] 

g(l)a(2)6(3)a(4)  + r6„(l)6„(2)6„(3)6„(4)  - 6„(1)6„(2)6„(3)6„(4)  + 6_(1) 

jugjju  gugu 

6 (2)6„(3)6„(4)  - 6,,(1)62(2)6„(3)6„(4)1  B(l)6(2)a(3)a(4)} 
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Denote  each  space  portion  bv  and  each  spin  portion  by  s^. 
Then  the  wavefunction  takes  the  form 


1 

(A-4)  \l)  = + foS-  + f„s-  + £,s,  + f_s_  + ffS-} 

2 2 3 3 4 4 5 5 6 6 

where  s^  = a(l)a(2)B(3)6(4) 
s^  = a(l)B(2)a(3)6(4) 
s^  = a(l)6(2)6(3)a(4) 
s^  = 6(l)a(2)a(3)6(4) 

Sj  = 6(l)a(2)6(3)a(4) 

Sg  = B(l)6(2)a(3)a(4) 

The  may  easily  be  obtained  from  direct  comparison  with  the  previous 

form,  Eq.  (A-3) . Thus,  in  a more  compact  notation, 

1 6 

(A-4.1)  t = — I 

24  i=l  ^ 


Since  'll*  = for  this  vravefunction,  the  ground  state  energy  is 
given  by 


(A-5)  E = ]i(;*H>j;dV  = 


6 6 

( I )H(  I f s )dV 
i=l  j=l  J 1 


Also, 
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(A-5.1)  E = 


6 6 


( I I f Hf  s s )dV 
i=l  j=l  ^ J ^ J 


since  the  Hamiltonian  function  H is  spin  independent.  But  s^s^  = 6^^ , 
so  that  the  double  summation  reduces  to  a single  sum: 


(A-5.2)  E=-^ 


(66  . f 6 

( I I = ^ ( I f .Hf  )dV 

i=l  j=l  ^ J ■'  i=l  ^ ^ 


Within  the  integral,  the  electron  numbers  may  be  permuted,  as 
they  are  only  dummy  variables.  If  the  order  of  integration  and 
summation  are  interchanged,  the  result  is  obtained  that 


(A-6) 


f^Hf^dV 


Returning  to  the  f^,  electron  numbers  may  be  permuted  and  nuclear 
labels  a and  b interchanged  (as  they  indicate  dummy  variables  in 
integrals) . Thus 


(A-7)  E = ^ I 


i=l 


f^Hf^dV 


or 


This  is  the  result  desired,  and  is  given  in  the  main  body  of  the  text 
as  Eq.  (9),  where  the  Hamiltonian  H has  been  expanded. 
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GROUPING  OF  THE  ENERGY  INTEGRAL 


From  Eq.  (9), 


(A-8)  E = i {f  I f dV 

^1=1  ai  ^ 


} ^ i>j=l  ’^Ij  ^ ^ 


f3L*fidV} 


Consider  first  the  last  term: 


(A-9) 


fi*fldV 


[i.g(l)(!.^(2)4.^(3)({.g(4)  - <|.g(l)4.^(2)4.g(3)4.^(4) 
+ <t'^(l)<(>g(2)<^g(3).^^(4)  - ({>^^(l)'^g(2)(J.^(3)(|.g(4)j  dV^dV^dV^dV^ 


For  simplicity,  denote  these  terms  as  ^3^3 » ^nd 


Then 


f • 


(A=9.1) 


fl^dV  - 2 - 2 


2 

f 

2 

2 

flldV  f 

fl2dV  + 

fl3dV  -h 

/ 

f,2fi3dV  - 2 


fl4dV  - 2 , 


fl3fl^dV  + 2 I f^^f^3dV  + 


^12^4^" 


The  Integral  of  a diagonal  terra  is  definitely  non  zero:  explicitly 


(A-10) 


<{.g(l)<^^(2)(j.^(3)(|.g(4)  dV 
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^g(l)<l)g(l)dV^  ^^(.2)^J2)dV^  (f.^(3)4.^(3)dV2  'l'g(4)(}.g(4)dV^ 


= (1)(1)(1)(1)  = 1 


and  similarly  for  the  other  diagonal  teirms. 


Of  the  off-diagonal  terms,  (choosing  the  first  term  as  represen- 


tative) 


(A-11)  I = |[4>g(l)4>^(2)<l.^(3)4.g(4)j|<^g(l)*^(2)4.g(3)4.^(4)JdV 

= I <{.g(l)<}.g(l)dV^  I <^^(2)<|.^(2)dV2  I (j.y(3)4.g(3)dV3  I <^g(4).l.^(4)dV^ 


■»-(3>*g(3>dV3  - Y„  j K - \ 


* 

= N N 
g u 


■yr,  -yr,  -ar„  -ar,  rT~T 

(e  ^-e  ^(e  ^ + e ^ ^ dV 

IT 


-NU  3!25I  ' 

g u TI 


W I -(«+Y)r^  -ar^  -Yrj,  -Tr^  - 


- e e “ + e e ‘')dV 


vOTJ  f -(a-hr)r 
N N { e ^dV  - 

g U IT 


-yr  -or, 
a b ..o 


J ® ^dV  - 


-ar  -yr 
e ®e  °dV  + 


e e dV/ 


In  the  second  and  fourth  integrals,  interchange  the  labelling  of 
atomic  sites.  Then  Eq.  (A-11)  becomes 


t. 


(A-12)  N N -2_1_  { 

g U IT 


U9 

‘dV  - I e 


dV  - 


-ar  -yr. 

8 b,,, 

e e dV 


-yr,  -or 
e e dV 


and  it  is  Immediately  seen  that  the  first  and  second  sets  of  integrals 
add  to  zero,  pairwise. 

Thus 


(A-13) 


fllfl2dV  = 0 


and  similarly  for  the  other  off-diagonal  terms. 
This  gives  the  result  that 


(A-IA) 


f j^f  j^dV  = 4 


Next  consider  the  electron-nucleus  interaction  terms.  Denote 


2 2 2 

H = - 7^  ~ — — 

en  i r . r.  , 

ai  hi 


These  terms  are  then 


4 

1=  1 

1=1 


f*H  f . dV. 
i en  i i 


But  since  the  electron  number  does  not  affect  the  result,  this  is  equal  to 


f,  H , f,  dV, 
1 enl  1 1 


(A-15) 


1 = 4 
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^ ^ I ^11  \nl^^ir^l2‘*'^13"^l4^‘*^l 


- 4 


^12  ”enl^^ir^l2’^^13"^14^^'^l 


+ 4 


^13  ^enl^^ir^l2'*‘^13“^14^‘^'^l 


" ^ ^ "14  ^enl^^ir^l2'^^13“^14^‘^'^l 


In  this  integral,  only  one  electron  is  affected  by  the  operator. 

In  the  cases  of  off-diagonal  terms,  such  as 


'll  “.„1  ^13  “’l 


the  components  in  each  case  break  down  into  integrals  of  the  form 


4 rr 


(A-16)  I f 
i=l  ^ 


9^(2)(^g(2)dV2  J <i>^j(3)(^g(3)dV^  j 


^g(4)(^^(4)dV^ 


<j.g(l)<^^(l)dV^ 


^<2)\„l^g(2)dV2  j 


'^u(3)<^g(3)dV3  <).g(4)<j.^(4)dV^ 


^g(l)<^u(l)dVl  j 


(2)dV2  ♦„(3)H^^^.K3)dV3 


♦ (4)*  (4)dV, 

^ Cl 


i (1).^  (l)dV  , 
g U 1 J 


^,(2)^g(2)dV2  J 


4.^(3)<^g(3)dV3 


In  each  such  circumstance,  there  will  always  be  a pair  of  component 


integrals  of  the  form 
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<{.  (2)6  (2)dV, 
u § ^ 


One  of  the  two  integrals  will  always  be  zero.  Thus,  for  the  same 
reasons  as  in  the  first  case  considered,  the  off-diagonal  terms  must 
be  zero,  and 


(A-16.2)  T = 4 


f,,H  ,f,,dV-  + 4 
11  enl  11  1 


^12^enl^l2‘^'^l 


+ 4 


f.H  ,f,,dV,  +4 
13  enl  13  1 


f,,H  Tf,,dV, 
14  enl  14  1 


This  can  be  further  expanded,  noting  that  the  integrals  over  all 
electrons  except  the  first  are  now  unity,  as 


I = 4{  I ^g(l)\„i6g(l)dV^  + j 


g enl  g 1 


I « 4{2 


(A-17)  I = 4{2H  + 2H  } 
g u 


Turning  finally  to  the  last  of  the  terms  in  Eq,  (A-8) , 


(A-18)  J = { 


fl*(  I -^)f,dV} 

i>j=l  ""ij  ^ 
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In  this  integral,  unlike  the  previous  cases,  there  are  terms  remaining 
from  many  integrals  which  were  previously  zero.  Equation  (A-18)  may 
be  expanded  as 


(A-18.1)  J = 


^^ir^l2'^^13“^14^ 


^^12  ^^13  ’^14  ^23  ’^24  ^^34 


^ ('ir'i2-^'i3-'i4^‘'^i 


J = 


11 


2,2,2,222 

^12  ^^13  ^14  ^^23  ^34  ^^24 


(fir^2-^^3-'i4>‘’''i 


-1 


^12 


_2_4._2_4.^4.^4.  1-  + ..^ 

"12  "13  "14  "23  "24  "34 


(fir'i2-^fi3-'i4>‘'" 


■13 


"l2  "l3  "l4  "23  "24  "34 


^^ir^l2‘*‘^13"^14^'*^ 


14 


"12  "l3  "l4  "23  "24  "34 


^'ir'i2-^'i3-fi4^‘'^ 


This  expands  into  a total  of  ninety-six  terms  which  must  be  regrouped. 


It  is  readily  apparent  that  the  four  diagonal  terms  are  combinable, 

4 

I 

as  are  terms  of  the  form 


f-^(  ^ )f  dV.  After  one  combines 

i>3=l  13 


terms,  Eq.  (A-18.i)  may  be  written  as 


(A-19)  J = 4 


‘11 


"12  "l3 


"l4  "23  "24  "34 


flldV 


+ 2{(-) 


11 


2 , 2 , 2 , 2 , 2 , 2 

"12  "l3  "l4  "23  "24  "34 


'■12^" 


\ 

U-i 

■-1-  + 

2 2 2 2 2 “ 

J 11 

_"l2 

"l3  "l4  "23  "24  "34 

fl3dV 
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"12  "13  ’^U  ’^23  "24  "34 


fl4dV 


■12 


"l2  "l3  "l4  "23  "24  "34 


^12  [r; 


2 + ^ + . _L. 

12  "l3  "l4  "23  "24  "34 


f,,dV 


fl4dV 


■13 


"12  "l3  "l4  "23  "24  "34 


fl4^V} 


Proceed  by  explicitly  expanding  the  diagonal  term  first: 


11 


^g(l)<(.^(2)4.^(3)<^g(4) 


"12  "l3  "l4  "23  "24  "34 


<(.g(l)<{.^(2)<J.^(3)<}.g(4) 


dV 


■11 


<{.  (1)41  (2)(-^)4  (1)4.  (2)dV,dV_ 
^g  'rj^2  8 “ 1 2 


+ If  ♦ (!)♦  (3)(— )♦  (1)«  <3)dV,dV,  , 


♦„(2>'»„(2)dV2 


* (4)t  (4)dV 
6 6 ^ 

<^^(4)4>^(4)dV, 

o o 


4.^  (1)  (4)  (-^) 4-„  (1)  4'„  (4) dv  dv, 

o a S 5 i H 


4.^(2)4>^(2)dV2 


4.^(3)4.^(3)dV3 


^u^2)^,(3)(7^.)  <^^(2)^^(3)dV2dV3 


4)g(l)4>g(l)dV^ 


<fg(4)4>g(4)dV4 


4>,.(2)4.„(4)  (~)4^j(2)4>g(4)dV2dV^ 


u g 


■24 


4>_(l)4>_(l)dV 

60 


4u(3)^^(3)dV3 


4>^(3)4'g(4)(— )4^(3)4'g(4)dV3dV4  j 


4-g(l)4>g(l)dV^ 


4.^(2)4.^(2)dV2 
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Here,  all  integrals  not  involving  a (-^)  factor  are  unity. 

’^ij 

others  may  then  be  combined  into 


/A_on\  T = f 
x-v/  '■ 


A (i)A  (4)(-^)a  (1)^  (4)av,dv. 

-S  "g  1:14  -g-  -g-  • 14 


♦u<«*u«>'7^>*u«H„(3)dV2dV3 


II 


+ 4 II  4.„(1)<J.,  (2)(-^)(^^(l)<j.,  (2)dV  dV, 

g u ri2  g u 12 


Expanding  the  next  term  in  Eq.  (A-19)  gives 


12 


'll 


_1_  + _2_  + -2-  + -2_  _2_  _L. 

^12  ^13  ^^14  ^23  ^ ^24  ^ ^34 


fl2dV 


j .<frg(l)^^(2)<^^(3)^g(4) 


I (—  + - 
^^12 


2 2 2 2 2 , 
13  ’^14  ’^23  ^24  ^^34 


<l>„(l)<f>„(2)6  (3)f  (4)  dV 

_ S u & u ^ 


12 


(3)<^„(4)  (-^)<}>„(3)6  (4)dV-dV 

« 8 ^24  8 ^ ^ H 


where  the  terms  involving  factors  of  the  form 


(A-13.1) 


<l>u(3)<l'g(3)dV3  = 0 


The 


have  been  dropped.  Similarly,  the  third  term 


13 
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f + _2_  + _2_  + ^ + _2_  + _2_)f  dV  = 0 

^12  ’^13  ’^lA  ’^23  ^24  ^34 


since  there  will  always  be  terms  of  the  form  of  Eq.  (A-13.1)  remaininR. 
The  fourth  terra  of  Eq.  (A-19)  is 


14 


f (_2_  + _l_  + _2_  + _l-+-i_  + dV 

^12  ^^13  ^14  ’^23  ^24  ^^34 


14 


<i>„(l)<t>.,(2)  (-^)f  (l)<!.^(2)dV  dV 
g u u g X 4, 


12 


after  cancellation  of  terms. 

The  next  term  in  Eq.  (A-19)  is 


23 


+ — + — + -^  + -^  + — )f  ,dV 
^12  ^13  ^14  ^23  ^24  ’^34 


23 


-}-  (l)<fru(2)  (:^)<{.^a)<i>  (2)dV^dV2 


'12 


Continuing  with  the  next  terra  of  Eq.  (A-19), 


24 


'12  'l3  'U  '23  '24  '34  I'' 


124  -O 


since  each  term  has  a factor  of  the  form  of  Eq.  (A-13.1) 
The  final  term  in  Eq.  (A-19)  is 


34  J 13  rj_2 


^ + — + — )f.,dV 
23  ^^24  ^^34 


rombiniuR  all  the  non  zero  terms  into  (A-19)  gives 


(A-21)  J '^11  ~ 20i2~-^13'^'^14‘^'^23"'^24'*'‘^34] 
j = 4 J,,  - 2[4  J,.] 

" ^Pll  " ^ 

The  term  is  given,  in  the  notation  of  the  main  text,  as 

(A-20.1)  '^gg  "^uu  ^ ^ ^gn 

and 

(A- 22)  3^2  ~ ^ ^gu 

The  combination  of  each  of  these  terms  into  the  energy,  given 
as  Eq.  (A-8),  gives  Eq.  (9)  in  the  main  body  of  the  text. 
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/U:riL11UJ.A  D 

EVALUATION  OF  INTEGRALS  ENCOUNTERED  IN  THE  GROUND  STATE  ENERGY 

CALCULATION 


Much  of  the  material  here  is  also  presented  in  Slater  (51) , his 
Appendix  Six  in  particular  containing  much  background  information  on 
these  calculations.  One  term,  '^gg3»  will  be  presented  in  detail  as 
representative  of  the  method  of  solution,  and  to  present  the  method  of 
evaluation  of  one  type  of  problem  Integral. 

CONSTANTS  AND  PRIMARY  INTEGRALS 

The  wavef unctions  $ and  'I'  are  normalized  wavef unctions,  that  is 
8 u 


$ $ dV  = 1 
8 8 


and 


$ $ dV  = 1 . 
u u 


The  N and  N are  the  normalization  constants.  They  are  given  by 
8 ^ 


r * 

$ $ dV  = 
8 8 


+X,)N  (X  +X.)  dV 
g a b g a b 


or  Ng  (X/Xi,)(x^+X^)  dV  = 1 . 


Thus 


XaXa  dV  +2j  x^X^  dV  +J  Xj^X^  dV  } = 1 . 


But  the  atomic  wavef unctions  are  already  normalized,  so  that 


XaXa  dV  = X^X^  dV  = 1 
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This  gives  the  relation  that 


N {1  + 2 
S 


XaXb  dV  +1  } = 1. 


So 


N { 1 + 
g 


X^Xb  dV  } = 1/2  . 


(B-1) 


Denoting  /x  Xu  = S,  then 

& D 

Ng  = (-^)(l/2)  , or 

N = ( — - — 

g ^ 1+S  • 


Similarly, 


N { 1 + 
u 


x^xi  <1V  1 - 1/2  . 


These  atomic  wavefunctions  are  also  normalized.  Then 

1 


(B-2) 

where 


N = 


u ■ ^2  ( 


T = 


Xa^b 


It  still  remains  to  determine  S and  T. 


S = 1 X^Xb  dV  = - 


^-ara^-arb 


dV 


(B-3) 


TI 


/'•i2  2\  “«B.X  ,,  , ^irR^, 

e dX  dy  (-^)  , 


-1  1 

where  a conversion  to  confocal  coordinates  has  been  made  to  facilitate 


the  evaluation  of  this  integral.  Here,  it  is  easily  seen  that 


~ 3p  3 

S = V-  ' 


,1 


dy 


-1 


-oRX, , 
e dX 


- 

-1 


-oRX 
e dX 
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Each  of  these  integrals  is  easily  evaluated;  the  final  value  for  S is 

S = ( 1 + „R  + ) . 

Or,  letting  w=aR  , 


(B-4)  S = ( 1 + w + 1^)  e"'' 


In  an  identical  manner. 


(B-5) 


T = 


^a'^b  Ti 


b dV 


It  is  immediately  seen  by  comparison  with  Eq.  (B-3)  that  the  two 
integrals  are  identical  except  for  the  interchange  of  one  constant  for 
another  (i.e.,  y a )•  The  result  for  Eq.  (B-5)  is  thus  given  by 

T = ( 1 + yR  + . 

Letting  x=yR>  this  may  be  rewritten 

(B-6)  T = ( 1 + X + ) e"^  . 


EVALUATION  OF  J „ INTEGRAL 
ggS 

The  J „ integral  will  be  worked  out  in  detail,  to  illustrate  the 

gg3 

technique  involved  in  the  solution  of  the  terms  given  in  Chapter  IV, 
and  to  demonstrate  the  method  of  evaluation  of  the  ill-behaved  integrals 
encountered  in  this  and  other  terms. 

The  problem  that  arises  with  the  similar  integrals  is 

that  the  integration  must  be  carried  out  in  confocal  coordinates,  where 
the  inter electronic  separation  ( — is  a complicated  expression. 

Its  value  is  given  as  Eq.  (23)  in  Chapter  IV. 


The  terms  it  will  be  necessary  to  deal  with  using  this  expansion 
are  of  the  form 


jfd)  [|  £(2)(^)  dvj  dv^ 


SO  one  must  first  evaluate  the  Integral 
(B-7)  = f(2)(-^)  dV^  . 

^19 


■)  dV^  . 


Only  the  case  of 

f (2)  = Xa(2)Xb(2)  = e 

need  be  considered,  as 

f(2)  = Xa(2)Xt(2)  = ^ e"Y^’'a2+^b2) 

is  adequately  treated  by  the  interchange  of  y “•  Using  the 


expansion  of  Eq.  (23),  Eq.  (B-7)  becomes 

(B-7.1)  I [ f UX2  -y2)(|)  ^ [ (-l)“(2k+l) 

88  ^ J„  i,  J,  2 2 R Z 


0 '11 


(k- 

)! 

2 

(k+ 

m 

)!_ 

pi"'  [Kail  qW  [A(bI|pl"lui)pi"lp2)a-^«2-'*l>e-“‘‘^2 


dX2dy2d<i>2  ^ 

The  integration  over  ({>2  will  yield  zero  except  when  m=0,  since 
2Tr 

a-^'“d*  - (^Xe"  - e-2^’)  . 


Physically,  the  vanishing  of  Eq.  (B-7.1)  except  for  npO  arises  from  the 
fact  that  an  S-type  wavef unction  has  no  dependence  on  azimuthal  angle. 
Thereby,  Eq.  (B-7.1)  simplifies  to 
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(B 


-7.2)  I = a3R2  I (2k+l)f  f (Xj  -uj  ) e“«^^2p  [x (a)_| Q [X (bH 

SS  k=0  il  M ^ ^ 


T»  \ T>  /%.  \ jy  j.. 


Now  the  factors  involved  in  the  y integration  are  expressable 

2 

in  terms  of  Legendre  polynomials  by  use  of  tabulated  values: 

(X2  -p2  ) = ( x2-i)  P (p  ) -|p,(y  ) 

22  2-^  02^22 

Making  use  of  the  orthogonality  properties  of  the  Legendre  polynomials, 


I becomes 
gg 


(B-7.3) 


I = a3R2{ 
gg 


(X2  - ^)e  “^^2p^[X(a)"|QQ[x(b)I|PQ(pi)  dX^ 


4 

3 


P [X(a)]Q2[X(b)J  P (ui)e'“^^2  ^^2  } 


1 


Substitution  of  tabulated  values  fo  P^,  P^,  Qq,  and  Q^,  and  reexpression 
of  the  integral  in  terms  of  the  actual  functions  Xj  and  X2  in  place  of 
the  dummy  functions  X(a)  and  X(b)  yields 

(B-7.4)  I = a3R2  ^ £n(Y^)  f \x2  - ■^)e"“^^2tjx2 
gg  J , 2 3 

00  ^ 

+ f (X2  - h e'“^^2d;^ 

Jxi  2 3 X2-I 

- f (y2-  -j)  ’ {(Xf  - j)  -en(^)  -2Xi}|^1"“^^2(x2  _ i)  dX2 

00  1 

- f o*'-  ¥ - ¥ I J '‘^2  > 


where  the  subscript  on  y has  been  supressed  as  superfluous. 
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I 


i 


i 


The  process  of  solution  which  will  be  followed  Is  one  of  breaking 
these  integrals  up  for  separate  handling.  The  first  of  these  is 

(B-8)  J ^^^2  ~ 3 ^ ’ 

This  integral  may  be  evaluated  easily  as 

<B-8.U  = (^)  { ( 1 -foR  + 2^)  - e"“'“  (1+t.RX 

+ x2.i)  ) , 

Continuing, 

00 

(B-9)  I-  = f (X^  - h e"“^^2  dX 

2 J 2 3 X2-1‘'  2 

From  the  symmetry  of  the  logarithmic  term,  it  is  obvious  that  there  are 
two  integrals  of  identical  form  to  be  considered,  each  depending  only 
on  the  sign  of  the  constant  in  the  argument.  It  is  therefore  possible 
to  replace  the  constants  +1  and  -1  in  this  term  by  the  single  constant 
-a,  and  by  later  setting  a=H-l  or  a=-l  obtain  the  specific  form  desired. 

Integration  by  parts,  followed  by  a change  of  variables  in  Eq.  (B-9) 
gives 

(B-9.1)  flt(Xi-a)e‘““‘[l-k.RXl+2^U5  -i)] 

+ e'“®(l-aaE  + i)  Ei(-aR(Xi-a))  + SS-S^  ^-«R(Xi-a) 

+ ^ <1  + «RX,  + aaR)  ) ! 

I 

1 

I 

Evaluating  Eq.  (B-9.1)  for  the  two  values  of  a allows  the  result  to  be  ■ 

! 

obtained  for  I„  that  | 

I 

I 

i 

I 


5 
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~2ij2  -1 

(B-9.2)  I-  = {1  +ctRX  + t)>  iln(\+l)-  MX-1)} 


2 
-aR 


+aR 


+ (1  +aR  + Et(-aR(A-l))  - ||p-  (1  -oR 


2e 


-aRX 


+ ^)  Ei(-aR(X+l)) 


Substitution  of  the  values  of  the  integrals  and  into 
Eq.  (B-7.4)  for  I is  now  possible.  Note  that  the  coefficient  of  the 
£t(-aR(X-l))  term  is  just  the  overlap  integral  S. 

Define  a new  fu  iction 
+aR 


S'  = e 


( 1 -aR  + 


a^R^ 


■) 


as  the  coefficient  of  the  Et(-aR(X+l))  term.  Substitution  of  values  for 
and  I2  then  gives 

(B-7.5)  I = f 'CnC-—)  - S'  E^(-aR(X+l)) 


+ S £c(-aR(X-l))  - aRe"“^^}  + - (p2-  i)  {XS  + X + 


Equation  (B-7.5)  is  the  initial  integral  for  a good  many  of  the 
integrals  involved  in  the  ground  state  energy  calculation. 

Next,  consider 


(B-10) 


'gg3  J^a'-'^b'^'^gg  '*■1 


JXa(l)x.a)I„  dV 

) 


J = — |e 
ggJ  ir 


I dV^ 
gg  1 


Transformation  to  confocal  coordinates  gives  the  form 

4 

{(»2  -i)  - (p2  -i)) 


(B-10.1) 


a^R^ 


gg3 


-1  •'1 


dX  dp 
gg 
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vrtiere  has  been  added  and  subtracted  in  the  (X^  - term  in  order  to 
obtain  the  factor  (y^  - -j),  which  is  directly  proportional  to  P2(y)* 
This  allows  the  orthogonality  properties  of  the  Legendre  polynois LmIs  to 
be  fully  exploited  in  simplifying  the  result.  Using  the  relation  that 


/-i  =.  |P2<U) 

a nonzero  contribution  to  the  integral  over  y arises  only  from  the 
products  between  terms  not  involving  y with  like  terms  and  from  products 
between  terms  proportional  to  y^  . Then,  noting 


(B-10.2) 


8_ 

45 


gg3 


60%^ 

5 

2a3R3| 


(X2-  i)  4(X)  dX 

I ^ 

7 -aRX,  ^ -2aRX-  aR 
I e X S - e ( j- 


1 


where 


4(X)  = S -en(^)  -S’  E4(-aR(X+l))  +S  Et(-oR(X-l))-aRe 


-aRX 


It  is  desirable  to  investigate  the  Improper  integrals  first.  Consider 

the  integrals  over  the  &c(x~l)  functions.  It  is  necessary  to  proceed 

by  an  integration  by  parts,  with  the  lower  limit  replaced  by  6 , since 

each  integral  is  improper  at  the  lower  limit.  Again  let  -a  be  used  in 

place  of  the  +1  and  the  -1,  and  investigate  the  integral 
/ 00 

(B-11)  J e"“^^(x2-|)  &t(-aR(X-a))  dX  . 

Integration  by  parts  gives 
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(B-11.1)  4^  E-C  (-aR(6-a))e  “^^{l+a6R  + i)  } 


a^R 


a^R^i  -a 


“ E^(-2aR(6-a)){l  + aaR  + -^}  e 


. 1 -aaR,  „ . 9„?.  -2oR(6-a) 

+ e ( aR  »•  aa^'R^)  e 

1 -aaR  -2aR(6-a) 


(1  + 2aR  -2aaR) 


A+l\ 


Combination  of  Eq.  (B-11.1)  with  the  integral  over  Eq.  (B-10.2) 

(noting  the  difficulty  with  the  lower  limit  in  both  cases) , and  taking 
the  where  this  limit  is  well-defined,  gives 


(B-10.3) 


^ ^ {s[s  tyi(2)  -S  ^(6-1)  -S'  Ei(-2aR)  -aRe 
+S  E-c(-aR(6-l))]  - S'[S  BC(-2aR)  -S'  U(-^aR)  + | e“^“^ 

+S[s  &c(-aR(6-l))  -S  E<.(-2aR(6-l))  +(|  + |^)e"“^  ] } 

f " -2aRX,,,  U ,,  2Sa3R2  f “ -aR\  ,, 

•j  e (X^-  j)  dX ^ — j e dX 


6a‘'R 


4p3| 


2a3R2 


1. 


” aR  , , . aRX^  ^ -2aRX  ,, 
( -^  + X + — ) e dX 

ton 


Since  the  terms  for  which  „ , are  not  well-defined  all  have  the  same 

o-*"! 

2 

coefficient,  S , they  may  be  grouped  together,  as 


(B-12) 


A = {-  Zn(6-1)  +2  Et(-aR(6-l))  - Et(-2aR(6-l) ) } 


By  expressing  the  EL  function  in  terms  of  the  integral  which  serves  to 
define  it,  and  rearranging,  it  may  be  demonstrated  that  this  term 
behaves  in  an  acceptable  fashion  in  the  limit  6->-  1. 


(B-12.1)  A = -p"  { - ^h(6-1)  -2 

o-> 


“X  J 

e dx 


/CO  _x 

e dx 


aR(6-l) 


} 


2aR(6-l) 
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Transform  the  integrals  in  Eq.  (B-12.1)  by  adding  and  subtracting 
the  same  quantity  to  the  equation: 


(B-12.2)  A = ^n(6-l)  -2 


+2 


5->l 

dx 

X 


e dx 


aR(6-l) 
dx  , 


4 ^ 

•'aR(6-l) 


dx 

X 


dx 


aR(6-l) 


2aR(6-l) 


2aR(6-l) 


2aR(6-l) 


Next,  combine  the  terms  into  single  integrals,  giving 


(B-12.3) 


-2 


A = ^n(6-l)  +2 


dx 


e dx 


, 1 


aR(6-l) 


(l-e~^) 


dx  + 

X J X 

' 1 2aR(5-l) 

+2  tn(6-l)  - ln(2)  - £n(aR)  - In(S-l) 


-X 


— +2  ^(aR) 

X 


Since  the  terms 
(B-12.4)  A = 


£n(5-l)  add  out,  then 


,1 


■(i^dx  - 

X 


0 


W -X  , 

e dx 


+ ^(aR)  - Inez) 


The  first  two  integrals  serve  as  a definition  of  Euler's  constant,  C. 
Ttius  A = C + in(w)  - ln(Z)  . 

llie  terms  not  yet  evaluated  are  elementary;  the  final  result  for 


(B-10.4) 


J o = { S2  (c  + Zn(w))  + (S')2  Et(-4aR)  -2SS'  £4(-2oR)  } 

gg3  5R 

-Zw,  5 23w  6w2  2w3  v 


APPENDIX  C 


COMPARISON  OF  RESULTS  OBTAINED  WITH  THOSE  OF  HUZINAGA  FOR  He-He 

In  the  paper  of  Huzinaga  (51) , detailed  values  are  given  for 
the  individual  integrals  only  for  the  case  of  R = 1 Bohr,  and  are 
given  for  three  sets  of  alpha  and  gamma  values.  These  have  been 
tested  explicitly  for  agreement,  as  have  other  values  given  by 
Huzinaga  in  the  form  of  the  energy  of  the  state  for  a given  R. value. 
The  following  table  presents  the  comparison  of  the  results  obtained 
by  both  methods,  for  the  value  R = 1 Bohr.  The  results  of  Huzinaga 
have  been  converted  to  Rydbergs  from  Hartrees  (each  equal  to  two 
Rydbergs) . 
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TERM 

HUZINAGA 

THIS  CALCULATION 

Alpha  = Gamma  = 1.75 

H 

88 

-8.080478 

-8.08132  (alpha=1.75C 

^uu 

-5.3238 

-5.3238 

“^88 

1.800708 

1.90114 

•'uu 

1.970024 

1.970024 

1.871924 

1.872140 

^8U 

-0.513378 

-0.513482 

^total 

-8.476912 

-8.47750 

Alpha  = 2.25;  Gamma  = 1.25 

^88 

-8.63269 

-8.632688 

«uu 

-4.97937 

-4.97937 

*^88 

2.29911 

2.299109 

•^uu 

1.49188 

1.491879 

1.72451 

1.724513 

^8u 

0.438298 

0.38693 

^total 

-9.41168 

-9.37126 
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TEM 


HUZINAGA 


THIS  CALCULATION 


Alpha  = 2.5;  Gamma  = 1.0 


H 

gg 

-8.69146 

-8.6914597 

«uu 

-4.55023 

-4.550229 

■^gg 

2.47925 

2.479253 

1.22889 

1.228894 

•'gu 

1.53580 

1.535799 

^gu 

0.33758 

0.21041 

^total 

-9.29956 

-9.18196 

All  values  given  are  for  an  intemuclear  separation  R of  R=1.00  Bohr 


APPENDIX  D 

CHARACTERISTICS  AND  PRODUCTION  OF  V , H,  AND  F COLOUR  CENTERS 

iC 

There  exists  extensive  literature  dealing  with  colour  centers, 
principally  the  F center  (39-41),  but  also  dealing  with  the  and 
H centers  (38,  59-65).  These  are  not  the  only  types  of  colour  cen- 
ters observed  experimentally,  even  in  alkali  halides,  but  these  are 
the  centers  of  immediate  interest  in  relation  to  the  breakdown  problem. 
These  centers  are  the  prototype  colour  centers;  many  of  the  other 
identified  types  are  related  directly  to  these  centers,  and  are  clus- 
ters of  two  or  more  simple  centers.  Ivey  (66)  has  written  a monograph 
on  electroluminescence  and  related  effects  which  is  also  of  interest 
in  this  connection.  In  this  section,  the  interest  will  be  focussed 
principally  on  the  characteristics  of  the  Vj^,  H,  and  F centers  and 
their  production  mechanisms  which  are  most  directly  related  to  the 
breakdown  problem.  For  a more  complete  treatment,  the  work  of  Town- 
send and  Kelly  (41)  and  that  of  Schulman  and  Compton  (40)  make 
excellent  starting  points  for  the  recent  literature,  while  the  paper 
of  Dienes  (38)  is  an  excellent  short  discussion  of  the  production 
mechanism  of  both  and  the  paired  H and  F centers.  The  discussion 
which  follows  will  follow,  in  large  measure,  the  works  of  both 
Dienes  (38)  and  that  of  Townsend  and  Kelly  (41) . Additional  information 


171 


172 


concerning  the  dynamics  of  the  interactions  of  these  centers  with 
the  host  lattice  is  given  by  Stoneham  (68). 

The  F center  consists  of  a negative  ion  vacancy,  in  which  an 
electron  is  trapped  to  neutralize  the  charge  of  the  region.  The 
center  is  normally  produced  as  an  F (also  known  as  an  a center) , and 
then  traps  an  electron.  It  is  optically  the  least  energetic  of  the 
centers  mentioned,  and  has  the  longest  lifetime  at  room  temperatures. 
The  F center  is  the  dominant  colour  center  in  the  alkali  halides. 

At  temperatures  above  about  150  K,  it  is  possible  in  some  materials 
for  an  additional  electron  to  be  trapped  at  the  F center,  giving 
rise  to  an  excess  local  negative  charge  within  the  region  of  the 
center. 

The  V,  center  is  of  much  more  interest  in  the  discussion  of 
k 

breakdown.  This  center  consists  of  a hole  trapped  by  a pair  of 
negative  ions.  After  trapping  the  hole,  the  combination  of  atoms 
relaxes  toward  a common  center,  in  effect  forming  a molecular  ion; 
this  molecular  ion  is  perturbed  by  the  presence  of  the  lattice.  In 
the  conditions  encountered  before  breakdovm,  the  Vj^  center  can  be 
formed  directly  in  the  perfect  lattice  by  impact  ionization.  To 
quote  Dienes  (38) : 

Ionizing  radiation  knocks  off  an  electron  from  a Cl 
ion  in  (KCl)  with  the  electron  pushed  off  into  the  conduction 
band.  The  original  Cl”  becomes  a neutral  chlorine  atom  (Cl'^) . 
This  configuration  is  unstable,  however,  with  respect  to  a CI2" 
molecule.  _This^  the  Cl®  grabs  part  of  the  electron  of  the 
nearest  Cl"  and  forms  a Cl  “ molecule  [as  illustrated  in  Fig.  5]. 
Thus,  the  center  is  a CI^”  ion  sharing  two  normal  lattice 
sites . 

As  shown  in  Fig.  5,  the  molecular  ion  will  orient  itself  along  [110] 
directions  in  the  normal  lattice.  Under  the  influence  of  a strong 
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electric  field,  however,  this  orientation  will  be  perturbed  by  the 
electric  field. 

The  H center  is  the  complement  of  the  F center;  it  consists  of  a 


Cl^  molecular  ion  occupying  a single  Cl  site.  This  center,  illustrated 


in  Fig.  21,  is  oriented  similarly  to  the  V center,  along  the  [110] 

iC 


direction.  Quantum  mechanical  calculations  by  Dienes  (57)  give  the 
right  configuration  of  the  center  but  the  wrong  orientation.  The  orien- 
tation of  the  H center  is  apparently  stabilized  by  interaction  between 
the  center  and  the  halide  ions  lying  at  either  end  of  the  H-center's 
axis.  Dienes  (38)  indicates  that  the  H center  can  move  easily  along  the 
[110]  direction  with  an  activation  energy  in  the  range  of  0.1-0. 2 eV. 
This  motion  may  be  significant  in  the  breakdown  process.  The  analagous 


motion  of  the  center  does  not  seem  to  have  been  studied,  but  it  might 


be  expected  to  be  similar  to  that  of  the  H center,  but  probably  having 
a somewhat  higher  energy  of  activation,  as  the  interaction  of  the  H 
center  with  neighboring  ions  would  tend  to  reduce  that  activation  energy 
for  motion  in  an  axial  direction. 

A formation  mechanism  for  the  paired  F and  H centers  has  been 
suggested  by  Pooley  (43)  whereby  the  energy  available  from  electron-Vj^ 
recombination  is  used  to  produce  the  pair  of  centers.  Dienes  (38) 
describes  the  process  as  follows. 

The  energy  available  for  the  formation  of  the  pair  of  centers  due 


to  the  recombination  of  the  Vj^  with  an  electron  is  about  9 eV  (for  KCl) . 


This  energy  is  the  sum  of  the  direct  electron-hole  recombination  energy 
of  about  6.4  eV  and  the  potential  energy  of  the  two  ions  once  the  recom- 
bination has  taken  place  (2.6  eV) . If  the  energy  is  distributed 


J 


174 


Fig.  20.  The  F center  in  an  alkali  halide  crystal.  This  center 
consists  of  a halide  vacancy,  occupied  by  a free  electron. 


(After  Townsend  and  Kelly) 


Fig.  21.  The  H center  in  an  alkali  halide  crystal.  Note  that 
four  halide  atoms  are  occupying  the  lattice  spaces  normally  occupied 
by  three  halide  ions.  (After  Townsend  and  Kelly) 
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asymmetrically  between  the  two  atoms,  it  is  possible  to  initiate  a 
replacement  chain  of  collisions,  moving  the  Cl  with  the  greater  kinetic 
energy  into  the  next  cell,  thus  forming  a closely  linked  F and  H pair. 

As  indicated  earlier,  the  energy  of  activation  for  motion  of  the  H center 
once  created  is  small,  so  the  motion  will  continue  along  several  cells, 
leading  to  a displacement  of  the  H center  by  several  unit  cells  from  the 
F center.  It  is  necessary  that  the  F center  retain  the  excess  electron, 
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as  the  energy  of  formation  of  the  pair  of  H and  F centers  requires 
additional  energy.  The  replacement  sequence  will  occur  along  a [110] 
direction  within  the  crystal.  It  appears  likely  that  the  mechanism 
leading  to  the  formation  of  the  H and  F centers  will  dominate  other 
defect  production  mechanisms  for  higher  temperatures  (above  ]50  °K) . 

The  effect  upon  the  H center  of  the  application  of  an  electric  field 
of  the  magnitudes  involved  in  breakdown  should  be  similar  to  that 
described  above  for  the  V center,  but  H center  mobility  is  probably 
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higher  than  that  of  the  center.  Motion  of  an  H center  entails  only 
a slight  displacement  of  the  added  Cl°  atom  to  an  adjoining  Cl  site. 

Thus  mobility  along  crystallographic  directions  close-packed  with  Cl 
ions  is  favored  over  other  directions  in  highly  perfect  crystals.  In 
crystals  with  high  elastic  stress  or  high  dislocation  densities,  other 
paths  might  be  preferred.  This  provides  perspective  on  the  directional 
effects  seen  by  Davisson,  Cooper,  and  others  (mentioned  earlier)  in 
breakdown  studies. 

There  exists  one  other  polarization  effect  which  calls  for  mention 

in  connectio.j  with  these  centers.  Both  the  V,  and  the  H centers 
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